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In the Office Action several documents have been cited that will be considered 
irrelevant to the average skilled person. 

Since it is believed that this might stem from the mentioning of terms used in the US 
patent application, like "seed laser", "modulation", "phase control" and "frequency-shifted 
feedback-cavity", certain aspects of the physics involved will be briefly explained to show 
that the principles and objectives of the present invention are entirely different. In addition, a 
list of peered publications analyzing the principles in detail is given. 

Only the work using an FSF-laser (Shattil, Phillips and Takeuchi) seems to warrant 

more detailed comments in the official reply as the average skilled person surely and easily 

realizes that the teachings of the other patents are completely unrelated to our scheme. 

L THE PRESENT INVENTION AND THE 
UNDERLYING PHYSICAL PRINCIPLES 

1. Major Elements of the current invention are ... 

a) ... for the laser source: 

a cavity, a gain medium, a frequency shifter (usually an acousto-optical 

modulator), 

a narrow band seed laser and a modulator in the path of the seed laser (with the 
Vattei being the key of the invention, see below), 
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b) . . . for the measurement set-up: 

a beam splitter, which separates the laser beam into two parts, one sent along a 
reference path, the other one sent to the object, the distance of which is to be 
measured, light collecting elements which collect light reflected (or scattered) 
back from the reference surface and the surface of the object, means to direct 
the reflected radiation towards a photo-detector followed by suitable 
electronics. The beam splitting and recombination is typically, but not 
necessarily, done in a set-up similar to a Michelson interferometer. 
The combination of detector and electronics records the light power. 
2. Comments on the major elements: 

a) The cavity and the gain medium 

The combination of cavity and gain medium exhibits net effective gain over a 
certain bandwidth. The larger the bandwidth, the higher the resolution and the signal 
enhancement due to the phase modulation (see below). Because radiation circulating 
in the cavity will always suffer from losses, "net effective gain" means that the gain 
medium itself must allow for an amplification factor larger than unity. 

The frequency shifter imposes a frequency shift A S hinper roundtrip on radiation 
circulating in the cavity. The shift leads to a comb of spectral components (order of 
10 4 components are typical), separated by A sh m, which extends over a certain 
bandwidth away from the seed frequency. Depending on the sign of A S hsn the comb 
can either extend to higher or to lower frequencies. The value of A S hift is chosen 
different from the free spectral range A FS r of the cavity. For A. sh if t = Afsr mode locking 
will occur, what is to be avoided. 

Normally, the effective gain is positive over a certain spectral range and 
negative outside this range. Typically, the seed laser frequency is set to a region with 
small positive (or negative) gain, but near one of the frequencies for which the 
effective gain is zero. The frequency shifter then shifts the frequency towards the 
other frequency, for which the gain is zero. As the frequency is shifted, the spectral 
components are amplified and reach their maximum power at the frequency (away 
from the frequency of the seed laser) for which the effective gain is zero. The comb 
will extend beyond that frequency, but the components will be successively attenuated 
Technical considerations lead to the conclusion that it is best, but not required 
for the method to work, that the net gain (i.e. "gain" minus "losses") is positive near 
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and around the seed wavelength and zero over an as large as possible bandwidth 
beyond the seed wavelength in the region covered by the spectral comb induced by 
the frequency-shifter. The power of the spectral components remains unchanged over 
the spectral range with "zero net gain". 

Such a cavity with a gain medium and a cavity-internal frequency shifter is 
known as a "frequency shifted feedback (FSF)" cavity. When the device shows 
positive effective gain, lasing will occur and the system is called a FSF-laser. 

This FSF-laser set-up is a central element in our ranging scheme, but is not 
claimed per se . The FSF-laser concept has previously been described by Kowalski 
Furthermore, the use of an FSF-laser (although without external phase-modulated 
seed) for ranging has been published by Nakamura et.al.. 

b) The measurement set-up: 

The Michelson-type setup is also central to our ranging scheme but - by itself - 
is not part of the patent claims, as it is a well-known text-book type concept. 

c) The effect of the seed laser and its modulation: 

Injection of a seed laser into the FSF-cavity will modify the laser operation as 
the output spectrum of the FSF-laser will be dominated by the frequency comb 
starting at the frequency of the injected photons. Without seeding, the spectrum is a 
combination of many frequency combs, each one started by a spontaneously emitted 
photon within the gain medium. Seeding the FSF-cavity without phase-modulating 
the seed will not improve the performance of the system for ranging. 

The seed laser is referred to as "narrow band*', but no specific limit to the 
bandwidth is set. It is the phase-modulation of the seed laser, in combination with the 
FSF-laser cavity and the Michelson-interferometer type setup, which leads to an 
orciers-of-magnitude improvement of the performance in distance measurement. 

Understanding the physics requires a detailed analysis of the radiation field 
reaching the detector. The spectral characteristics of the field (frequency, amplitude 
and phase) are determined by the features of the FSF-laser, the seed laser including its 
phase-modulation and the relevant propagation distances. The field at the detector is a 
superposition of the two fields delayed in time by the difference in propagation time 
along the reference path and the path to the object. 

The signal at the detector can be represented as a sum of terms [see eq. 22 in 
Opt. Commun. 282, 2212 (2009).]. Each term, labelled by an index n, has an 
amplitude and a time dependent phase. The detailed analysis of the latter term shows 
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that for a specific modulation frequency Q = Q* of the phase of the seed laser - and 
only for this modulation frequency- all terms have the same phase. 

Theory furthermore shows the relation of Q* to the desired path difference 
AL. The experimental signature for finding Q* is the maximum of the signal as Q is 
varied. When Q is different from Q*, then the phase of the (order of) 10 4 terms 
changes from term to term. 

The situation can be visualized by interpreting each one of the term (of eq. 22, 
see above) as a vector with amplitude and phase (the phase being the angle in a 
complex plane). Adding many such vectors with varying phase will lead to a sum- 
vector with is small. However, when the phases of all terms are equal, the vectors all 
line up parallel and the sum pf the vectors will increase by orders of magnitude. 

In conclusion, it is the phase modulation of the seed laser radiation and the 
option to vary the frequency of it that allows the finding of the frequency Q*, which 
leads to a dramatically enhanced signal and thus leads to the determination of AL. 

The phase modulation of the seed laser, which is the key element of this 
invention, works towards a dramatic improvement of the ranging scheme when 
combined with a FSF-laser and a Michelson-type optical arrangement. This scheme is 
entirely different from any other scheme which was referred to by the US patent 
officer, even when terms such as FSF-laser or seed laser or phase modulation are used 
in the cited documents. 
II NON-OBVIOUSNESS OF THE PRESENT INVENTION 

The reasons for and the consequences of imposing phase modulation on the seed are 
far from obvious or trivial. 

As a senior research professor at a technical university I am well familiar with the 
knowledge and experience of the average skilled person such as a PH.D. student or a 
postdoctorate. According to my conviction, the average skilled person would not have been 
able to find and implement the scheme proposed, not even in a trial and error process, based 
on the knowledge of an FSF-laser and the ranging concept, as demonstrated by e.g. 
Nakamura and any of the cited documents. 

The physics behind this particular ranging scheme is subtle, but well understood by 
now. Even to a trained person with relevant experience (not to speak of the average skilled 
person) the reason for the dramatic improvement of the signal for ranging with an FSF-laser, 
when phase modulated narrow band radiation is injected into a FSF-laser cavity, is not at all 
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obvious. Most importantly^ it can not be derived from any of the teachings of the cited 
documents, as the underlying physics is, in each case, entirely different. It is the complexity 
of the physics involved which classifies the phase-modulation scheme as a true invention. 

The suggested scheme was found only after several years of work on FSF-Lasers and 
their use for ranging applications at the University of Kaiserslautem and after extensive 
discussions among the inventors in Kaiserslautem including a careful analysis of the work by 
Nakamura et.al., who showed that an FSF-laser (without seed) can be used for ranging, 
although with rather poor performance. 

The present scheme was first identified by rather involved analytical and numerical 
calculations. It was only after the availability of those results that the intuitively 
understandable physics, underlying the concept was identified. Based on those discussions 
and results, the physics behind the success of bur ranging concept with orders-of-magnitude 
improvement over ranging with only an FSF-laser is now well understood. 

Only detailed expert-analysis of how the phase-variation, imposed on the seed laser, 
propagates through the components of the set-up and how it shows up in the sum of the two 
field at the detector, reveals the underlying physics. 

The detailed scientific analysis extends far beyond the scope of the present remark 
and has been subject of a number of scientific papers, all published in the mean time in 
reknowned journals after peer-reviewing 

Among the coauthors are the following inventors which all are prominent scientists : 

Yatsenko: Ukrainian State Prize Winner 1998 for his work in laser-based 
metrology 

Shore: Humboldt-Research Award winner 1997/1998, 

25-year employee of Lawrence Livermore National Laboratories , 

specializes in optics and laser physics 

single author of the "bible" of coherent excitation: 

"The Theory of Coherent Atomic Excitation", John Wiley 
(1990), 2 Volumes, Total of 1 700 pages. 
Bergmann Senior research Professor of the Technical University of 
Kaiserslautem, 

Chair to the "Laserzentrum Kaiserslautem" between 1997 - 2002 
Chair to the special doctoral training program("Graduiertenkolleg") 
"Nonlinear Optics and Ultrafast Physics" between 2002-2006, 
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Chair to the center for "Optical Technologies and Laser- controlled 
Processes" between 2002 - 2008 

has coordinated the Marie-Curie Trainingsite : Coherent Control in 
Matter Waves and Quantum Optics 

1972 Godecke Forschungspreis for the thesis work (topic: one of the 
first applications of lasers to collision dynamics), awarded by the 
University of Freiburg and sponsored by local industry (Godecke). 
1981 Physik-Preis (Physics - Prize) of the Deutsche Physikalische 
Gesellschaft (German Physical Society) "for his ground breaking work 
in laser state selection and the application of the method to scattering 
experiments using crossed molecular beams" 

1992 Fellow of the American Physical Society (fellow status awarded 

"for his development of innovative laser systems and their application 

to molecular scattering processes") 

2000-2004 Coordinator of European Research Training 

Networks: "Coherent Control of Molecular Systems and Processes - 

COCOMO", 

2006-2010 Coordinator of European Research Training Networks: 
"Engineering of Quantum States of Matter and Light - EMALI",. 
2003 Robert-Wichard-Pohl-Award of the Deutsche Physikalische 
Gesellschaft (German Physical Society) "in recognition of his 
outstanding scientific achievements, especially in molecular physics, 
which had an impact on a broad range of fields in physics" 
2003 Honorary Doctorate of the Latvian University of Riga 
2003 Max-Planck-Researchprize for International Collaboration 
awarded jointly by the Alexander von Humboldt- Foundation and the 
Max-Planck-Society "for outstanding research achievements" 
2005 Honored by the Ministerpresident of Rheinland-Palatinate Kurt 
Beck on December 8th with the "Verdienstorden des Landes 
Rheinland-Pfaiz" (order of merit). The "Verdienstorden" is the highest 
honor awarded by the state of Rheinland-Palatinate. It is given for 
outstanding service for the state and its citizens. 
The following attached publications in peered journals explain the theory of the 
present invention: 
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1. L.P. Yatsenko, B.W. Shore and K. Bergmann 
"Theory of a Frequency Shifted Feedback Laser" 
Opt. Commun. 236, 183 - 202 (2004) 

A new and thorough theoretical treatment of FSF-laser, however not yet 
including the competition between coherent radiation evolving from 
spontaneous emission and injected photons. In particular: showing the 
equivalence of the two pictures "stationary comb" and "moving comb". 

2. L.P. Yatsenko, B.W. Shore, and K. Bergmann 

"Ranging and Interferometry with a Frequency Shifted Feedback Laser" 
Opt. Commun., 242, 581 - 598 (2004) 

First discussion of the principles of laser ranging based on a FSF-laser, 
emphasizing the relevance of injection seeding by frequency- (or phase-) 
modulated narrow band radiation. 

3. V.V. Ogurtsov, L.P. Yatensko, V.M. Khodakovskyy, B.W. Shore, G. Bonnet, 
and K. Bergmann 

"High Accuracy Ranging with Yb 3+ -doped Fiber-Ring Frequency-shifted 
Feedback Laser with Phase-modulated Seed " 
Opt. Communications 266, 266-273 (2006) 

First discussion of experimental ranging results based on an FSF-laser with 
injection by phase-modulated narrow band radiation, confirming the theory 
presented in paper 2. 

4. V.V. Ogurtsov, L.P. Yatensko, V.M. Khodakovskyy, B.W. Shore, G. Bonnet, 
and K. Bergmann 

"Experimental Characterization of an Yb 3 * -doped Fiber Ring Laser with 

Frequency-shifted Feedback" 

Opt. Commun., 266, 627-637 (2006) 

Experimental characterization of an FSF-laser with injection seeding, 
confirming theoretical results of paper 1. 

5. V.V. Ogurtsov, V.M.Khodakovsky, L.P.Yatsenko, G.Bonnet, B.W.Shore, and 
K. Bergmann 

An all-fiber frequency-shifted feedback laser for optical ranging; signal 

variation with distance 

Optics Commun. 281, 1679 - 1685 (2008) 
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Experimental verification of a specific feature of the ranging method using a 
FSF-laser with injection of phase-modulated narrow band radiation, namely 
the variation of the signal with distance. 

6. L.P.Yatsenko, B.W.Shore, and K.Bergmann 

Coherence in the output spectrum of frequency shifted feedback lasers 
Optics Commun. 282, 300-309 (2009) 

The full ab initio theory of FSF-laser operation with treatment of growth 
starting from spontaneous emission and from injection seeding on equal 
footing. 

7. L.P.Yatsenko, B. W.Shore, and K. Bergmann 

An intuitive picture ^optical ranging using frequency shifted feedback lasers 
seeded by a phase modulated field 
Optics Commun. 282, 2212-2216 (2009) 

Explaining the physics behind the dramatic enhancement of the signal by 
injection seeding of phase-modulated radiation on the basis of an intuitively 
accessible picture. 

These documents explain the theory behind the invention should the explanation given 
above should not be considered sufficient. 

Ill COMMENTS ON PATENTS CITED BY THE U.S. PATENT OFFICER: 

As indicated in the introduction to this remark, only the work using an FSF-laser 
(Shattil, Phillips and Takeuchi) seems to warrant more detailed comments as the average 
skilled person surely and easily realizes that the teachings of the other patents are completely 
unrelated to our scheme. 

a) Documents unrelated to the principles of the present invention 
Palese: 6 570 704 B2 "High average power chirped pulse fiber amplifier array" 
This patent is about a high power fiber amplifier allowing the production of 
nearly arbitrary pulse shapes. The term "seed" is used but it means the radiation 
which is to be amplified after spectral dispersion by an array of fiber amplifiers. The 
terms "phase" and "phase modulator" are used but in an entirely different context. 
Like in a liquid crystal pulse shaper (also called modulator) the phase of the di fferent 
spectral components are controlled (NOT modulated) by phase modulators in order to 
achieve a certain pulse shape. 

Goldberg: 5 745 284 "Solid-state laser source of tunable narrow-bandwidth 
ultraviolet radiation" 
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I have no idea why this patent is cited in the given context. The patent is about 
frequency conversion starting from a diode laser, which is usually considered too 
weak for efficient conversion into the UV. Talk is about "seed laser" (3/7) and the 
"amplitude modulator" (3/15). The seed laser is very standard in the context of 
frequency conversion, as one typically starts from a source with low power and 
narrow bandwidth, which is then amplified and converted. The amplitude modulation 
serves to increase the peak power without enhancing the average power of the diode 
laser radiation in order to increase the efficiency in the non-linear conversion step. 
Thus, the teaching of this text is not even remotely related to our FSF-scheme for 
ranging. 

Gabl: 5 592 327 "Regenerative amplifier incorporating a spectral filter within 
the resonant cavity" 

Like the Goldberg text, the teaching of this invention is absolutely unrelated to 
our approach. The invention is about pulse amplification of radiation, called "seed". It 
describes means to allow for high amplification factors without reaching destructive 
intensity levels. This is done by stretching the pulse of the radiation to be amplified in 
time and recompress it after amplification. Also, the inclusion of a spectral filter is 
proposed to prevent or reduce amplification of spontaneous emission, which would 
reduce the gain in the wanted spectral region. Thus, the teaching of this text is not 
even very remotely related to our FSF-scheme for ranging. 

Mocker: 5 394 238 "Look ahead windshear detector by filtered Rayleigh and/or 
aerosol scattered light" 

Windshear detection relies on the measurement of the frequency shift of 
reflected (incoming) light with respect to the outgoing light. A simple mean to 
measure such a shift is to set the frequency of the outgoing light at the edge of a steep 
transmission filter and send the reflected light through the same or an equivalent filter. 
The frequency shift can be derived from the change of the transmission, provided 
appropriate normalization procedures are applied. The essential teaching of this patent 
is about improved ways of generating the outgoing radiation and locking it to a 
frequency filter. Except for mentioning the name "seed", it has no connection to our 
scheme and objective. 

b ) Documents closer related to the present invention 
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Shattil: 5 955 992 "Frequency-shifted feedback cavity used as a phased array 
antenna controller and carrier interference multiple access spread-spectrum 
transmitter" 

This patent is about controlling the emission properties (in particular the 
direction of emission) of an antenna array and the processing of received signals 
through such an array. Although it does rely on an FSF cavity but it has no modulated 
seed laser and has entirely different objectives and works on different principles. 
Ranging and phase modulation of a narrow band seed laser is not mentioned and not 
used. It does mention phase control (see e.g. column 12, line 50 and before) only in 
order to achieve mode locking, a regime which our scheme avoids. 
Phillips: 5 835 199 "Fiber-based ladar transceiver for range/Doppler imaging 
with frequency comb generator" 

This patent is about measuring jointly distance and velocity. Pulsed radiation 
is used with long pulse duration for high velocity resolution and wide spectral band 
width for high spatial resolution. A frequency shifter is used to derive radiation with a 
larger bandwidth (through discrete frequency components) from a narrow band laser. 
As in our case, the frequency shifter plays an essential role in their scheme. It is 
stated, however, that best performance is achieved when the phases of the spectral 
components of the comb are random (column 12, line 15 and 16/57) and 
incommensurate spacing of the frequency components is wanted (12/19). These 
statements alone provide clear cut documentation that not only the objective but also 
the approach is entirely different from ours, since well-set phases, controlled through 
the phase-modulation of the seed, are essential in our case. At 16/26 it is also said that 
ideally the spectrum should be of white light character. This scheme seems to rely on 
waveform detection and the ranging scheme is of the time-of-flight type. 
Takeuchi: 5 394 235 "Apparatus for measuring distortion position of optical 
fibers" 

Surprisingly the English of this text is of rather poor quality and suffers from 
ambiguous or unclear statements. The claimed invention is about an improved 
diagnostics for optical fibers. It relies on amplification through Brillouin scattering 
and is thus, in any case, limited to solid state material. In a similar (known) version of 
the scheme light with frequencies f 0 (pulsed) and f m (continuous) is coupled to one 
and the other end, respectively, of the optical fiber to be measured. When the 
difference frequency f 0 -f m equals the Brillouin shift frequency (given by phonon 

10 

BDDB01 5765466V 1 



modes in the material) the continuous laser is amplified by the other, counter- 
propagating, laser. The amplified radiation is displayed (T assume: as a function of 



time, thus correlating via suitable timing circuits any deviation from expected growth 
of the intensity to a location within the fiber, thus localising a defect region). The 
accuracy relies on a stable frequency difference. The essential teaching of this patent 
is how to best create the two frequencies. This is done by a frequency shifter which 
derives the frequency f m from the laser of frequency f 0 with the main benefit that 
frequency fluctuations of the radiation source do not affect the relevant frequency 
difference. As stated above, this technique relies on non-linear mixing within the 
medium and has no element of phase control or phase modulation as in our scheme. 

In summary, none of the cited documents relate to the physical principles underlying 
the present invention - 

The enhancement of the measurement signal obtained by the use of a modulated seed 
laser could not be expected by an average skilled person and was predicted only after very 
thorough analysis and long discussions among a number of very experienced scientists. 

I hereby declare that all statements made herein of my own knowledge are true and 
that all statements made on information and belief are believed to be true and further that 
these statements were made with the knowledge that willful false statements and the like 



so made are punishable by fine or imprisonment, or both, under § 1001 of Title 18 of the 
United States Code and that such wilful false statements may jeopardize the validity of 



the application or any patent issued thereon. 
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Abstract 

We present a new approach to the description of the output from a frequency-shifted feedback (FSF) laser seeded by 
a phase-fluctuating but stationary continuous-wave (CW) laser. We illustrate the new analysis by showing how short 
frequency-chirped pulses arise for appropriate operating conditions. We show the equivalence of two common view- 
points of the FSF laser output as either a moving comb of equidistant frequencies or as a fixed set of discrete fre- 
quencies. We also consider operation of a FSF laser when there is no external seeding laser, and instead the cavity 
radiation originates with spontaneous emission. 
© 2004 Elsevier B.V. All rights reserved. 

PACS: 42.55.-f; 42.60.Da; 42.55.Ah 

Keywords: Optics; Lasers; Frequency-shifted feedback 



1. Introduction 

Conventional lasers rely on multiple passes 
through a gain medium to reinforce a preselected 
frequency, thereby obtaining near-monochromatic 
output. For many purposes, such as excitation of 
Doppler-broadened atomic transitions, it is desir- 
able to have light spread over a broader range of 
frequencies. One means of fulfilling that objective 
is to introduce a frequency shifting element into 
the laser feedback loop such that successive passes 
of a wavepacket take place with different carrier 
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E-mail address: shore@physik.uni-kl.de (B.W. Shore). 



frequencies. The operation of such a frequency 
shifting feedback (FSF) laser has been considered 
by several research groups [1-46] and it has been 
used for a number of practical applications 
[12,21,36,41,47,48]. Although much is known 
about the behavior of FSF lasers and some of their 
surprising properties (see e.g. [20,24]), there remain 
some fundamental questions. Even the precise 
nature of the field emerging from such a laser re- 
mains a matter for discussion and confusion [33]. 

Although all authors agree on the elements that 
comprise a FSF device, there are conflicting ap- 
proaches concerning detailed description of the 
electric field E(t) that emerges from the device. In 
part this is because a FSF laser can operate in a 
variety of regimes: the output may range from 
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extreme irregularity to very regular pulse trains (cf. 
[24,46]), depending on such controllable properties 
as the amount of gain (controlled by the pump 
power) and the relative values of various time 
constants or frequency bandwidths of the device. 

In this paper we reexamine some of the prop- 
erties of this device, as deduced from equations 
that describe the time dependence of the (complex- 
valued) electric field amplitude. We also provide 
further clarification of the relationship between 
two apparently conflicting views of the output 
electric field. We are here concerned with steady- 
state operation, and will not consider the very 
interesting possibilities of irregular behavior. 

The terms "steady state" and "stationary" de- 
serve some comment when applied to pulsed la- 
sers. As an example, one can consider the 
operation of a conventional mode-locked laser 
whose output consists of a train of short pulses. 
On a short time scale (that of one pulse) the in- 
tensity, and even the frequency, depends on time 
and thus on this time scale the output appears 
nonstationary. Nevertheless, on a time scale that 
encompasses many pulses the properties of the 
output is unchanged, and one can consider this as 
a stationary process. 

We begin, in Section 2, with a general descrip- 
tion of a prototype FSF layout. In Section 3 we 
derive equations of motion describing the two 
coupled systems that comprise the laser, namely 
the electric field (including both a seeding field and 
spontaneous emission as well as the stimulated 
emission) and the atomic inversion. Our field 
equations deal with complex-valued amplitudes, 
not just with photon numbers, as in rate equation 
treatments (e.g. [24]). Thus our formulation of the 
FSF laser equations includes phases. These are 
important for the injection seeded FSF: all the 
pulse chirping effects considered in Sections 4—6 
rely on phase relationships. For the FSF laser 
seeded by spontaneous emission phase is not 
important. 

In Section 4 we present solutions to these 
equations, allowing for seeding by a noisy but 
otherwise continuous- wave signal. We show here 
some special cases and some specific examples, and 
we reconcile (see Section 6.2) two common but 
apparently conflicting views of the FSF laser (as 



discrete frequencies, Section 5 or as a moving 
comb of frequencies, Section 6). Section 7 briefly 
discusses operation of a FSF laser whose field 
originates entirely from spontaneous emission, and 
Section 8 connects the present formalism to the 
earlier rate-equation model. 



2. The FSF laser layout 

We consider a simple idealized model of a fre- 
quency-shifted feedback laser: a ring cavity of 
perimeter L in which there is a section having gain, 
some elements that provide spectral filtering, and, 
most importantly, an element that induces a fre- 
quency shift (of angular frequency A) on each 
wavepacket that passes through it. The round trip 
time is r r =L/c. For definiteness we assume that 
the frequency shift is provided by an acousto-optic 
modulator (AOM) and that this acts as a grating, 
to deflect light (into first order) that has been fre- 
quency shifted from a> to to + A. The light that 
remains in the undeflected beam (zero-order dif- 
fraction) leaves the cavity as output. 

We follow common practice and idealize the 
cavity field as having spatial variation only along 
the longitudinal coordinate z (i.e. we unwind the 
mirror action). We take the entrance face of the 
AOM to be z = 0 and regard the various optical 
elements, such as the gain medium and the AOM, 
as having negligible length. 

Fig. 1 shows a prototype schematic layout for a 
seeded FSF cavity. The input is through the AOM 
and the output is light that is undeflected by the 
AOM. With each pass around the cavity some 
light is lost. The most important loss sources are 
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seed _-ifflL_ output 





gain 






# 



Fig. 1. Schematic diagram of a ring cavity showing mirrors 
(M), frequency shifting AOM , amplifier (gain), as well as input 
(seed) and output beams. 
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the mirror transmittance and the loss to the output 
beam (through the AOM), both of which are fre- 
quency dependent. Often the cavity will include a 
filter that introduces controllable frequency-de- 
pendent loss. To simplify the presentation and to 
avoid unnecessary model-dependent details, we 
will consider only the case of a cavity whose loss 
varies negligibly within the frequency interval A. 
In this case the ordering of the elements (AOM, 
mirrors, broadband gain and filter) are not im- 
portant, and we can lump all losses into one ele- 
ment, which we call the filter. We consider this to 
be infinitesimally thin, time independent, and to 
decrease the field amplitude entering it with fre- 
quency co by a factor 2T(io) = exp[-/(co)]. We 
assume that the filter function /(co) varies slowly 
with frequency. Such broadband operation is 
found in almost all experimentally realized FSF 
lasers, and it is a prerequisite for use in high-res- 
olution interferometry. In keeping with this re- 
quirement, we assume that the effect of the losses 
can be expressed in quadratic form 

The filter function f(co) has its minimum value, 
f m , when co = cor, and it increases from that value 
quadratically with the frequency difference co — ct> f . 
The parameter T f characterizes the bandwidth 
of the filter while the central frequency co f pro- 
vides a convenient fiducial for defining frequency 
detunings. 

In order to sustain operation and produce a 
steady-state output, the cavity must include a 
section having gain, where amplification replaces 
the light lost to the output beam. As with the loss, 
we consider the gain to be localized in an infini- 
tesimally thin element, having the property that a 
field entering it with frequency co will be increased 
in amplitude by a factor exp[#(a>)]. We assume 
that the gain g(co), like the loss f(co), varies slowly 
and smoothly with frequency. 

3. The coupled field-atom equations 

As in all dynamical descriptions of a laser, we 
require an equation relating the time evolution of 



the electric field amplitude in the presence of po- 
larizable atoms (some form of the Maxwell equa- 
tions), and equations describing the time evolution 
of atomic properties in the presence of a given field 
(some form of a density-matrix equation). We here 
derive the needed coupled equations by consider- 
ing the passage of a pulse through the region of 
gain (the amplifier), and then considering the ef- 
fects of the filter and the frequency shift. 

Our model of the FSF laser differs in essential 
respects from some previous models, and readers 
familiar with those may question how one can 
consider the output of the FSF laser to be sta- 
tionary. Where we to consider (as we do not) just a 
single pulse in the FSF laser, then this pulse will 
evolve with each round trip, changing frequency 
repetitively and undergoing growth through gain. 
The output in such a scenario is clearly not sta- 
tionary, for each output pulse differs in amplitude 
and frequency from all others. 

Although it is useful to describe the frequency 
shifting properties of a FSF laser by following the 
succession of frequency shifts occurring to a spe- 
cific pulse, our mathematical treatment does not 
follow such a single pulse. Instead we consider a 
system in which the FSF laser is continuously 
seeded, either by spontaneous emission or by 
continuous injection of a seeding laser. We then 
examine the characteristics of the FSF output after 
startup transients have ceased. Because the seeding 
source is assumed to be stationary the FSF output, 
after some initial transient time development, is 
also stationary; operation of the FSF laser reaches 
a steady state in which the statistical properties of 
the output signal do not change. 

We assume that the properties of the filters, 
though frequency dependent, are independent of 
time and of the intracavity field. However, we al- 
low the gain to be influenced by both the pump 
power and by the cavity field; time variations of 
the cavity field will then induce time-dependent 
changes in the gain. 

We assume that the relaxation time of the gain- 
medium polarization is much shorter than all other 
relaxation times under consideration, so that we 
can use standard techniques of adiabatic elimina- 
tion and neglect the time derivative of the medium 
polarization. With this approximation we deal 
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only with rate equations for number densities or 
probabilities (populations), rather than with the 
Bloch-type equation for a density matrix. We as- 
sume that the population is spatially uniform 
within the (infinitesimally thin) gain medium, and 
that during one round trip within the cavity there 
is at most only a small change in the population 
(and so the gain can change only slowly with time). 

In keeping with our interest in Ti:Sapphire as 
the gain medium we idealize the amplifier as a 
collection of identical four-level atoms (cf. [49], 
whose notation we follow), with the lasing transi- 
tion between levels 1 and 2, while the pump radi- 
ation moves population between levels 0 and 3. 
Rapid decay from level 3 fills the upper lasing level 
2, and rapid decay from level 1 to the ground state 
0 ensures that the lower-level probability Pi (t) is 
negligible, and the population inversion 
w(t) = P 2 (t) -Pi(0 can be replaced by the popu- 
lation P 2 {t), 

The FSF laser is a system where the Fourier 
transform in its normal form does not appear 
particularly useful. As one learns in Fourier anal- 
ysis, sampling a signal for a time interval r gives a 
valid representation of the spectral properties of 
the signal if t is long enough so that the effect of 
the discontinuities at the beginning and end of the 
segment are insignificant. Then the Fourier trans- 
form, from which one evaluates the spectrum, will 
be independent of any change in t. 

Our goal is different. We want to obtain equa- 
tions which describe all regimes of FSF laser op- 
eration - steady state, transient, or pulse evolution. 

3.1. The field amplitude equation 

We take the center frequency of the effective 
filter, car, as a carrier frequency and express the 
time dependence of the complex-valued electric 
field E{t) at position z = 0, taken to be just prior to 
the AOM in a path around the cavity, by the 
complex envelope & (t) 

E(t) = S{t) exp(-k£>fO. (2) 

The FSF laser is one of many examples where the 
frequency content of a signal changes with time 
and conventional Fourier transform methods are 
not sufficient. This problem, of explicitly time- 



dependent spectra, has been discussed by Eberly 
and Wodkiewicz [51]. We overcome the difficulties 
by using a sliding-window Fourier transform 
(SWFT). For any fixed fiducial reference time T we 
introduce a Fourier decomposition of the envelope 
& (t), valid in the time window T - x T < t < T, 



/oo 
dtz7<?(tz7, T) exp(-iwt) 
oo 

for T - r r < t < T. 



(3) 

Note that the Fourier frequency tor is not the op- 
tical frequency; it is a detuning from the carrier 
frequency a> f . The windowed Fourier components 
are obtained from the integral 



(4) 



i r T 

g{m, T) / d/<f(0exp(+ku0. 

2k J T _ Xr 

It should be understood that our choice of the 
window duration as the round trip time is not re- 
lated to any requirements of Fourier transforms, 
nor is it to impose any phase closure condition, for 
example, a requirement that accumulated phase 
during a round trip must be an integral multiple of 
2tt. 

Because we deal with a SWFT rather than a 
conventional FT, the value of the frequency shift 
relative to the fiducial time T is not important. We 
obtain a set of finite-difference equations which 
can be used to describe the FSF laser. The small- 
ness of the frequency shift can subsequently be 
used to simplify these finite-difference equations as 
approximate differential equations. 

The intensity I{t) at z = 0 during the time in- 
terval T — r r < t < T is expressible either as the 
absolute square of the time-dependent electric field 
envelope or as a double Fourier integral 



c 

8n 



/OO roo 
J dnjj dm 2 g(m x , T)S{w 2 , T)* 
oo J -oo 



x exp[-i(n7] - w 2 )t\. (5) 

The windowed Fourier components £ (w, T), based 
on windowed Fourier analysis, serve as the basic 
quantities for theory: we derive the desired equa- 
tions by considering the connection between 
^(nj, T) and <?(m, T - x r ). Prior to arriving at the 
reference position, this Fourier component &(xd, T) 
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has made a circuit in which it was filtered, under- 
went a frequency shift, and underwent gain. It be- 
gan that circuit with frequency w — A at time 
T — x r , and after immediately shifting frequency it 
acquired a phase increment (co f + w)x T from spatial 
propagation. We express this behavior by writing 

<f (ro, T) = <${w - A,T - r r ) exp[+i(w f + w)t t 

+ g(CDf + tD, T) - f(CDf + W)). (6) 

We assume that there is little change in the gain 
during one round trip and we assume the filter has 
only a slight frequency dependence 

g(C0 f ~h tU, T T r ) — g((Of + m, T), 
f(ajf + ro + zl) ~/(oj f + m). 

To emphasize the discreteness of the periodic fre- 
quency shift it proves useful to consider a two- 
dimensional time-frequency space 

X = (ro, T), (7) 

in which successive circulations involve a change 
by the vector X 0 = (zl, t r ). Then the basic equation 
(6) for the field change in one round trip can be 
written as 

S (X) = <?(X-X 0 ) exp [G(X)] + e(X) + £{X), 

(8) 

where the effect of gain (and loss) is present 
through the function 

G(X) = i(a> f + w)x T + g(co f + m, T) - f{a> f + m) 

(9) 

and the two additional terms represent various 
sources of radiation (in Fourier space): 

• e(X) is the seeding field source (if present), 

• £{X) is the source of the spontaneous emission 
field during a time interval z T . 

Eqs. (8) and (9) are our basic equations for the 
field amplitude emerging from the FSF cavity. 
Missing from our treatment is any explicit nonlin- 
earity, such as four-wave mixing due to Kerr effect 
or gain saturation. Because we neglect such non- 
linearities our resulting theory cannot explain re- 
sults such as those of Ito and coworkers [29,32] 
where the chirped comb (though a very small com- 
ponent of the FSF laser output) was observed in an 
unseeded FSF laser. 

The source terms have the following properties. 



3.1.1. The seed laser source 

The seeding laser produces the source term 
e(X). We take this to be a source of continuous 
radiation at the frequency 

OJ s = OJf + tD s (10) 

(that is, tu s is a frequency offset from the central 
filter frequency cor of the carrier field). We allow a 
finite bandwidth to this laser by taking the win- 
dowed Fourier transform to have the form 

1 f T 

e(X) ~ e c — / dfexp[— \w s t - i(p s (t) + iwt]. 
2n J T _ Zr 

(11) 

The time varying phase cp s (t) appearing here rep- 
resents the phase of the seed laser. To permit 
treatment of the unavoidable finite bandwidth of 
any seed laser we treat this phase as a stochastic 
process. The quantity e c is the amplitude of the 
seeding laser field within the cavity. It parametrizes 
the seeded intensity I c within the cavity, 

/c=^M 2 - (12) 

3.1.2. Spontaneous emission source 
Spontaneous emission (within the gain medium) 

is described by the term £(X). This is similar to a 
Langevin force [50]. This frequency domain source 
term relates to the time domain source via Fourier 
integral within one round trip time 

Z(X) = Z( W ,T)=±J^ d^ sp (0exp(-hin7/). 

(13) 

The time- varying emission source £ sp (/) describes 
the field created from spontaneous emission 
events. These occur at random times, and succes- 
sive emissions are completely uncorrelated. This 
property will be quantified in Section 7. 

3.2. The population inversion equation 

To complete the description of the atom-field 
system we require an equation for the time de- 
pendence of the gain coefficient g(a>, T). This is 
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directly proportional to the population inversion 
w(/) = p 2 (t) - P\(t) ^ P 2 (t) of the lasing transition 
(assuming the two levels have equal degeneracies). 
We express this by the equation 

g(a>,T)=g 0 (co)w(T). (14) 

We take the frequency dependence of go{<o) to be 
of Lorentzian form 



w (T)-w(T-r r ) 



goM = 



git 



i + (cD-a> g y/rl 



(15) 



where a> g is the frequency of maximum gain and 
1/Fg is the relaxation time for the off-diagonal 
terms of the density matrix. 

To avoid unnecessarily complicated formulas 
we assume that the width of the frequency-de- 
pendent loss, as incorporated into the filter /(co), is 
much smaller than the gain width: T g » J>. With 
this assumption we can regard the gain as being 
independent of frequency, and set go(a>) = g m . 
This simplification is not essential, and the results 
are readily adapted to the case when only the gain, 
but not the filter, has significant frequency de- 
pendence. For simplification we have taken go((i>) 
and f(co) to be purely real functions of frequency. 
It can be shown that, in the quadratic approxi- 
mation, the inclusion of imaginary parts (disper- 
sion) would lead only to a redefinition of the round 
trip time. 

Within the approximations discussed here, the 
inversion w(T) changes from its value at the earlier 
time T - t r due to the following effects, each of 
which acts for time r r : 

• growth by pumping (into level 3 with rapid de- 
cay to 2) at rate R; 

• loss by spontaneous emission at rate 
y s w(T - *r); 

• loss by stimulated emission at rate yj{t)/l^ x . 
Taking these effects into account we write the 
following expression for the inversion after one 
round-trip time: 

w(r) = w(T - x r ) [l - y s / ^ dt - y s J + Rx T . 

JT-x r 'sat J 

(16) 

In the absence of lasing, I(t) = 0, this equation 
gives the stationary inversion w(T) = w(T — T r ) = 
w m = R/y s . Then the inversion equation reads 



= -y, Ur - r r ) - w m + W {T - t r ) f T ^ ^ 

I JT-x t 'sat T r 

(17) 

In traditional treatments of population inversion 
one considers the limiting case of an infinitesimal 
time increment (here the fixed increment x r ) 



w(T) - w(T - T T ) 



dw 
df 



and thereby obtains an equation for the time de- 
rivative of the inversion. However, in the present 
work the discreteness of the round trip time is an 
essential part of the physics of the field, and so we 
here retain also a similar discrete form for the 
evolution equation for the inversion. 

The time integral over intensity can be taken as 
the definition of the average intensity during a 
round trip 



I(T)=- f dtl(t) = f° dm/for), 

T r JT-x T J -oo 

where, in the frequency domain, 

T) = I" du/^m, T)S(w + w\ T)* 

87t 7-oc 

1 - expf-in/x,-] 



(18) 



itZ7'T r 



expfin/r]. 



(19) 



With this definition we can write the equation for 
the inversion as 



w(T) = w(T - r r ) 



1 -y^t-j T,T r 



(20) 

This equation completes the necessary set of 
equations for the FSF laser. 

3.3. Steady-state operation 

We are interested in steady-state operation of 
the FSF laser, as contrasted with irregular or 
transient behavior. We therefore carry out a sto- 
chastic average (• - •) of the inversion and the av- 
erage intensity 
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vv=<w(r)) = (>v(r-T r )>, I=(I(T)). 



(21) 



4. The FSF laser with (noisy) seeding 



This step rules out possible relaxation oscilla- 
tions of the inversion, but it does not require 
that the phase be constant. With this approxi- 
mation we deal with the average gain coefficient 

gin*. 

In this work we consider a laser whose spectrum 
has a width much larger than the free spectral 
range dco axia \ = 2n/x T . Under this restriction we 
can replace w(T) and w(r + t r ) with vv and I(T) 
with 7 in Eq. (20). Then we find that the stationary 
average inversion vv and the average intensity / are 
related by the equation 



vv = 



1+7//. 



(22) 



sat 



3.4. Solving the coupled equations 

The dynamics of a laser with frequency-shifted 
feedback can be described by two stochastic finite- 
difference equations: Eq. (8) for the time-depen- 
dent spectral component &(X) and Eq. (20) for the 
inversion w(T). The procedure for solving the 
equations is as follows. 

1. We find the solution S(X) of Eq. (8) with vv 
treated as an unknown parameter. 

2. Using £{X) in Eq. (19) we obtain I{w y T). 

3. We do the stochastic averaging of I{T) over re- 
alizations of the spontaneous emission source 
£(X) and (or) realizations of the stochastic 
phase q> s (t) of the seeding laser, and from these 
we obtain 7. 

4. We solve Eq. (22) for vv. 

5. Using the expression obtained for vv we obtain 
&{X) and the intensity I(T) of items 1 and 2. 
We also obtain the spectral density of the inten- 
sity, I(m, T). 

The following sections discuss implementa- 
tions of this program with various approximations 
and simplifications for the following important 
cases: 

• A laser with injection seeding, but without 
spontaneous emission, e(X) ^ 0, £(X) = 0. 

• A laser with spontaneous emission, but without 
seeding, e(X) = 0, Z(X) ^ 0. 



An important simple case of the FSF laser oc- 
curs when a seeding laser provides the only source 
term. We take this to be a CW laser, of frequency 
(o s , but we allow finite bandwidth by including 
phase noise (p s (t). 

4.1. The electric field 

We here neglect the spontaneous emission 
source ^(X), so that Fourier components of the 
laser field are determined by the equation 

g(X) = S{X -^ 0 )exp [G(^)] + e(X). (23) 

Here we use the function G(X), defined by Eq. (9), 
with stationary and frequency-independent gain 
coefficient g m w, so that it becomes 

G{X) = i(co r + w)x T + g m w - f(a)f + m). (24) 

As in Eq. (11) we take e(X) to have only phase 
variation (i.e. constant amplitude). 

The solution to the difference equation (23) is 
the infinite sum 



£{X) = e{X) + ^2 ex P 



/=0 



(25) 

i<Ps(<') 
(26) 



x e(X — nXo), 

where 

exp [ - ito s / / 

+ i( w -nA)t']. 

From the windowed Fourier transform of expres- 
sion (25), over the time interval T - x r < t < T, 
one finds (after straightforward but lengthy al- 
gebra) the time-dependent envelope of the field to 
be 

g{t) = e c exp[-hz7 s f - i(p s (t)] 

oo roo 

+ £c ^2 I dm exp ( —iwt) 

x r n {m) x q„(m), (27) 

where r n {w) incorporates all consequences of gain 
and loss 
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r n (w) = exp 



\n(co T + tzj s + [n + l]A/2)r r 



/=0 



(28) 



and ^n(w) expresses all the effects of phase fluctu- 
ations (seed-laser bandwidth) 

1 f T 

g n (rn) = — j dr'expf-K^f' - nx T ) 

+ i(ro-a7 s -nA)f\. (29) 

The function q„(xn) is maximum near w = m s -j- nA. 
The width of this maximum is roughly l/t r H- <5oj s , 
where 2n/x r is the free spectral range and Sco s is the 
bandwidth of the seeding laser. We assume that 
the width of the filter J> is much larger than these 

n » l/T r + (5a> s . (30) 

This equation says that the function r n {w) is a 
slowly varying function of the detuning w com- 
pared with the narrow function q n (™)- Thus one 
can replace r n (w) by its value at the point 
m = m s + nA. After integration over w and using 
the expression (1) one can write the electric field 
envelope &{t) as a sum over terms with constant 
real-valued amplitudes a n 



n=0 



- i(ro s + «/j)f]. 



(31) 



Section 4.2 provides a simple interpretation of this 
expression, and of the indexing integer n. The 
constant real-valued amplitudes are expressible as 
a n = exp^,,) with 

S n = n g m w-fm 

6xnl + 6w s nA + 6w s A + 2* 2 J 2 -h 3 A 2 n + J 2 

i2r 2 

(32) 

As defined earlier, / m is the filter loss at maximum 
transmission (a> = a> f ) and w s is the detuning of 
the seed laser frequency from a> f . The phase of the 
nth component has a stochastic portion <p s (t) and a 
time-independent part that is quadratic in n 



<P n = -T r n[a> s + (n + \)A/2]. 



(33) 



The steady population inversion vv affects only the 
amplitudes a n of Eq. (31), not the phases of the 
Fourier expansion (31). The inversion depends 
upon the steady mean intensity 7 through the 
condition of steady-state operation, Eq. (22). In 
turn, the mean intensity 7, 



/ = 8^(')l 2 = /c£l a «l 2 - 



(34) 



n=0 



involves vv through the dependence of a n shown in 
Eq. (32). Eqs. (22) and (34) combine to give a 
single transcendental equation for either vv or 7. 

4.2. Simplified interpretation 

For the specific case of a seeded FSF laser the 
expression (31) for the envelope of the intra- 
cavity field E(t) = <f(f)exp(-icOff) can be inter- 
preted in the following way. At any time t the field 
E{t) inside the cavity at z = 0 is the superposition 
of fields entering the cavity at earlier times that 
underwent at least one round trip. For simplicity 
we assume the incident field to be perfectly 
monochromatic at frequency a> S9 meaning that 
there is no random phase, <p s (t) = 0. Then the 
incident seed field (prior to passing through the 
input mirror) can be written as 



E in (t) = ^ s exp(-iov). 



(35) 



In this case the amplitude e c of the seeding source 
is £ c = ^ s ^"i n , where ^~ in is the transmittance of 
the input mirror for the incoming light. 

The cavity field entering the AOM at time / is 
expressible as the sum 



(36) 



where component E n (t) at time t has, by definition, 
made n round trips and has undergone n frequency 
shifts. The field emerging at time t as the output 
beam is evaluated by considering the transmit- 
tance of each beam, 3~ n = J(w s + «J), through 
the AOM without frequency shifting 



(37) 



n=0 
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The amplitude of an individual component is de- 
termined by the effective gain Q(w) 9 



w 2 ( w 2 \ 



(38) 



defined as the difference between saturated gain 
and filter losses. Here Q m is the maximum of the 
effective gain and wo is the larger of the two fre- 
quency detunings for which the effective gain is 
equal zero 

Q m = g m w -f mj w 0 = y/2Q^r r . (39) 

It is well known that for a conventional laser, in 
CW operation, the saturated gain is equal to the 
loss. However for a FSF laser the gain is not sat- 
urated completely, even when the laser operates in 
the stationary regime. 

The field components can be evaluated from 
the following consideration. The component E n {t) 
entered the cavity at time t — nx T> with frequency 
co s . It decreased in amplitude by the factor ST XTy * 
shifted frequency by A, and underwent attenua- 
tion or gain [depending on the sign of the effec- 
tive gain Q(co s + ^)] by a factor R\ = 
exp(£?(co s + 1^)). On the next round trip it un- 
derwent frequency shift from a> s + A to co s + 2A 
and again underwent attenuation or gain, this 
time by the factor ^ 2 = exp(£?(co s + 2A)). It 
makes n such round trips, on each of which it 
changes amplitude by a factor <%i(a>) and acquires 
a phase shift cox T because it propagates the dis- 
tance T r /c. It follows that the various components 
are 

E 0 (t) = <^i„exp[-ia> s f], (40) 
£,(f) = ^s^in^i exp[-i(ca s + A)t] 

x exp[+i(a> s + A)x T ], (41) 
E 2 (t) = <f s ^ in ^i^2exp[-i(a> s + 2A)t) 

x exp[+i(co s + A)x T + /(a; s + 2A)x T ], (42) 



E n {t) = S s y xn &^ 2 • • • @ n exp[-\(co s + nA)t] 



x exp 



\y^{qj s + mA)x v 



(43) 



(44) 



= exp 



(a>s + nA) 2 

2r 2 r 



(45) 



The first exponential factor describes the time 
variation appropriate to frequency co„ = co s + nA. 
The last exponential factor describes the spatial 
phase Yin (o> n /c)L acquired by propagating n full 
round trips through the cavity, each of length 
L = r r /c but at a different frequency a>„. 

Note that in the absence of any frequency shift, 
in the absence of the gain, and with 
@(co s ) = 0t < 1, the individual components have 
the form E„(t) = S s ^ m ^" exp[-iw s r + inco s x r ] and 
the field has the construction 



(46) 



= ^s^m exp[-iay] 
1 ) 1 - ^exp[+icL> s T r ] ' 

The maximum of this field occurs when 
x T co s = 2nq, for any integer q. This is the standard 
Fabry-Perot cavity condition for axial modes. 



5. The discrete frequency model 

Comparison of the heuristically derived Eqs. 
(35) and (44) with the field presented in Eq. (31) is 
done by introducing an envelope (?df(0 as in Eq. 
(2), 

E(t) = exp[-iw f f]<f DF (0. (47) 

Using Eq. (10) and the definitions (32) and (33) of 
S n and $„ we find 



n 

x exp[— i(zd s + nA)t) 
= exp[— ia> s f] y~]A n exp[-inAt]. 



(48) 



n=0 



where the reflectivity & n is 



This is just the amplitude of Eq. (31) with q> s = 0 
(i.e. a monochromatic seed) and e c = ^^"in- 

The amplitude <2W(0 defined by the summation 
(48) is periodic, with the period T p = 2k/ A, and 
hence the expression is in the form of a Floquet 
series: a time-dependent prefactor exp[— ico s t] 
multiplying a periodic function, expressed as a 
Fourier series with constant complex-valued 
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Fourier coefficients A n = ^ exp[-i# n ]. Therefore 
the field in the seeded FSF cavity is a comb of 
discrete frequency components, separated by the 
increment A y and starting from the seed frequency 
a> s . The phase <P n varies quadratically with n, 
whereas the exponent of S n , defined in Eq. (32), 
varies cubically. However, in practice the variation 
of S n is significant only near its maximum 

S ma * = £w where h a 1 ** 0 is quadratic. 

The magnitude of A n is governed by the real- 
valued exponent S n the sign of which is determined 
critically by the saturated inversion iv. In turn, w 
depends on the total intensity and hence on the 
seed detuning xd s and on other parameters of the 
laser. We will show in Section 5.1 that, for wide 
range of parameters, w is determined mainly by the 
filter width /> and the frequency shift A. Using the 
definitions of Eq. (39) for the peak effective gain 
Q m and n7 0 we rewrite the real-valued exponent S n 
of Eq, (32) as 



S„ = 



0 «(l- 



6ot? H- 6m s nA + 6w s A + 2n 2 A 2 + 3 A 2 n j- A 2 
6^0 



(49) 

If the seeding laser is tuned near the lower fre- 
quency where the effective gain is zero, so 
tds^ - G7 0 , then starting from n = 0, the value of 
S„ grows steadily with n as long as the gain exceeds 
the loss, i.e. Q(w s + nA) > 0. Eventually, for 
n > (wo + \ w s\)/d 9 the loss dominates the gain, 
and the exponent S n falls with further increase of n. 
Thus there is a maximum of S n , determined from 
the solution of the equation (d/dn)S„ = 0. This 
value, 5 max , occurs for the index n = [n max ], where 
[. . .] means the integer part. For large wo^> A the 
peak occurs at 

« (TZ70 + KD/A (50) 
This is a value of n for which Q(w $ + nA) =0. The 
variation of the magnitude a n = exp[S„] of the 
frequency component with index n is therefore 
very nearly of Gaussian shape centered around 
n = H max ; the exponential argument S n departs 
from quadratic form in n only when n is large or 
small compared with « max . 

Fig. 2(b) illustrates the variation of the com- 
ponents a n with n, normalized to unit value at the 
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Fig. 2. (a) The relative effective gain function Q/Q m vs. n> 
showing the definition of w 0 . (b) The magnitudes a„ of the field 
components E n (t) vs. n, normalized to unity at the maximum. 
The parameters \/{n w f = 240. For clarity only every second 
component is plotted. One can see, from comparing the two 
frames, that the peak of the component distribution occurs 
where the effective gain changes from positive to negative. 



maximum. The peak of a n is very close to the lo- 
cation of the value Q{w s + nA) = 0. For reference, 
Fig. 2(a) shows the effective gain Q{w) associated 
with these components. Fig. 5 of [8] shows an 
experimentally observed spectrum which looks 
exactly like our Fig. 2(a). 

5.7. Gaussian approximation 

The S n of Eqs. (32) or (49) follows a quadratic 
dependence on n for appropriate conditions. This 
gives the amplitude a n a Gaussian form. To iden- 
tify those conditions, we consider the field given by 
Eq. (31) when the seeding laser is detuned far away 
from the central frequency of the laser spectrum. 
In this case the first frequency component is very 
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weak and the shape of the laser output spectrum is 
close to Gaussian. To quantify this relationship we 
express S n as a Taylor series around n maxy where 
dS„/dt = 0, retaining only the quadratic term. The 
magnitudes a n can then be written as 



I ( n — flmax \ 2 



(51) 



where the width parameter n w obeys the relation- 
ship 



K) 2 = 



wo A 



(52) 



and /7 max is given by Eq. (50). The common am- 
plitude s/ is expressible as 



s4 = exp 



2xUq + tu^ — 3tu s n7Q 

6^ 



]■ 



(53) 



where the frequency parameter y is defined as 

y = (^) ,/3 . (54) 

The Gaussian approximation of Eq. (51) is valid if 
n w is large enough, meaning 

n w » rj/xnl. (55) 

This condition depends, through w 0 and Eq. (39), 
on the still-unknown saturated inversion u>, 



- A- , 

vv = — + 



(m) 2 



g m 2 gm (r r ) 



2 ' 



(56) 



To find vv we assume that the Gaussian approxi- 
mation is valid, so that the mean number n w of 
discrete components is large, n w » 1. Then we 
replace the summation over n in Eq. (34) by an 
integration 



(57) 



Taking into account the condition \A\ «: T f , we 
have from Eq. (22) the result 



exp 



2w\ ±wj-3 



3V 3 




(58) 



The small parameter 1 introduced here is 

2/3 



P ~ V2/S3, (t,-\)\Aj 



(59) 



(60) 



with the threshold parameter r\ being defined as 

ga 

Eq. (58) gives an implicit relationship between the 
value of xdq (and hence of vv ) and the parameters y, 
/? and ro s . Fig. 3 shows the numerical solution of 
this equation for different values of detunings iu s of 
the seeding laser as a function of the parameter /?. 
The dependence on ft is very weak. 

A simple approximation, based on inspection of 
this figure, is to take the parameters of the 
Gaussian profile to be 

With this approximation the amplitudes of the 
Fourier components depend only on the detuning 
A and the filter bandwidth T f , through the com- 
bination y. 

More detailed calculations, taking into account 
the parameters / c , I s and t], make only a small 
change. When P is extremely small (so 1 / ln(l / /S) is 
also small) the transcendental equation (58) can be 
solved approximately for the specific choice of 
detuning trr s = 0. In this case one obtains the ap- 
proximation 

y- (62) 



TJJq 



Wq 



3 1 1 

2^J 




Fig. 3. The shift w 0 of the position of the spectrum maximum 
(solid lines) and the width t of the spectrum (dashed lines) in 
units of y as a function of the parameter /?. The detuning of the 
seeding laser is: (1) w s = -8y; (2) tu, = 0; (3) w s = —2y. 
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The maximum of the spectrum is located at 

Wmax = <Of + Wq. (63) 

The mean number of components (the "width 
parameter" n w ) is given by the formula 

y/2 y 



[(3/2)ln(l//?)] ,/6 ^ 



(64) 



and the mean intensity / is 

7=*/«(if-l). (65) 

The spectral width, T — An*, and the mean in- 
tensity are both nearly independent of the seeding 
laser intensity. 

We have derived here, from first principles, 
formulas for the width and the peak of the dis- 
crete-frequency model. The important parameters 
are A, the frequency shift per round trip, and the 
filter width T f ; this latter has only a weak depen- 
dence on /?, which incorporates the saturation in- 
tensity, on the threshold parameter y\, on the 
intensity of the seeding laser in the cavity and on 
the ratio of the filter width to the AOM frequency. 
In principle the values of these quantities are 
known for any particular laser, and so the opera- 
tion of the FSF laser can be fully described. 

5.2. Explicit Fourier formula 

When the Gaussian approximation is justified 
for the amplitude of Eq. (48), the construction for 
the electric field E(t) reads 



E{t) = Gxp[-ico s t]e c s/ x ^ exp 
+ \x r n[oj s + (n + l)A/2] 



(n — rt max ) 



exp[— in At]. (66) 



Because the frequency width of the spectrum t is 
much smaller than the filter width T f one can ne- 
glect the frequency dependence of the AOM effi- 
ciency near a> max . Then the output field E oui (t) is 
directly proportional to the intracavity field E(t). 

Because we assume n w to be very large we can 
approximate /i max by the integer /i max ~ [n max \. For 
purposes of display, the largest-magnitude Fourier 
component, which occurs for n = /i max , serves as a 



useful reference point for expressing the others. 
We rewrite the argument of the exponential, 

S n - /#. = - ( " ~ " max)2 + h T n[co s + {n + 1)4/2], 

in terms of the difference between n and n max , using 
the variable m 



= ~{m/n w ) , 



>W+m = ~ m[oi s T r + (* max + l)Ar r /2] 



-m 2 (Ax T /2). 



(67) 



(68) 



We eliminate the occurrence n max from this for- 
mula by shifting the time t — ► t — to with an ap- 
propriately chosen initial time t 0 appropriately. 
The choice 



^o = T r [(n max + l)/2 + ^] 



(69) 



available when A ^ 0, eliminates the linear varia- 
tion with m from the phase, and eliminates also the 
dependence on the seed frequency co s . This choice 
leads to the expression 

E(t) = exp[-i(o; s + n max A)(t - f 0 )] 



x A' m exp[-imA(t - to)} 

m—-oo 

with 

A' m = A 0 exp[-(m/n w ) 2 - \m 2 {Ax r /2% 



(70) 



(71) 



where the new normalization constant A 0 incor- 
porates all constant phases and the factor e c j2/. 
With this choice the seed frequency co s appears 
only in the common time-dependent prefactor 
exp[-ia> s (f - t 0 )], but not in the summation. The 
prefactor does not affect the intensity of the output 
radiation, |£(/)| 2 , and it gives only a constant 
contribution to the instantaneous frequency Q(t), 
defined by 



Q{t) = — arctan(Re[£(0]/Im[£(0]). 



(72) 



The right-hand side of Eq. (70) is a Floquet series 
in which the coefficients are Gaussians with a 
phase that varies quadratically with the index m. 
Apart from normalization, the amplitudes depend 
on only two parameters: the phase change in one 
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round trip, Ax Ti and the width of the effective gain 
distribution, n w . 

5.3. Examples of pulsing output 

It is instructive to examine illustrative examples 
of the output, Eq. (70). For that purpose we in- 
troduce the dimensionless time variable x = t/T p , 
thereby expressing time in units of the basic repe- 
tition period T p = 2n/A. We then deal with the 
function 



(73) 



+M 

F(x) = F 0 ^ exp[-i/n2n(x - x ref ) - (m/n w ) 2 

m=-M 

-im 2 (m T /T p )}. 

The function F(x) depends on three parameter 
combinations: n Wy x re f = to/T p and r r /T p . The first 
of these determines the width of the Gaussian 
distribution of Fourier amplitudes, and thereby 
controls the number of components that can add 
coherently. The relative time * rcf merely adjusts the 
position of the pulse peaks along the x axis. The 
most important parameter combination is 
r r /r p = 2uAt t . Whenever the parameter t t /T p is an 
integer, the field is a sequence of short pulses, one 
per period. Such pulses are described by the 
function 

+M 

F(x) =F 0 ^2 exp[-iw27i(x - x rcf ) - (m//i w ) 2 ]. 

m——M 

(74) 

The duration of an individual pulse is inversely 
proportional to the number of Fourier compo- 
nents that contribute significantly to the sum. The 
present work parameterizes this number by « w , 
the width of the Gaussian distribution of Fourier 
components. Thus as the n w increases, the pulse 
durations can become shorter. The time-depen- 
dent frequency of such pulses, defined by Eq. 
(72), remains constant during the presence of a 
pulse. 

Fig. 4 shows examples of these short pulses, for 
several values of n w . Earlier such pulses were re- 
ported in [20]. For uniformity of plots we nor- 
malize the field by taking the quantity l/F 0 to be 
real and equal to the sum (we take the summation 
limits to be M = 4« w ) 




0.00 0.01 0.02 0.03 0.04 0.05 
t/Tp 

Fig. 4. The relative intensity /(/) of the output field vs. time t y in 
units of the period T p = 2n/A. Plots are for n w = 10, 20, 30 and 
40. 



l/^o= Yl exp[-(m/* w ) 2 ]. 



(75) 



= -M 



More generally, trains of short pulses occur 
whenever x T /T p is the ratio of two integers. These 
become less intense as they become more numer- 
ous in one period. Fig. 5 presents examples of such 
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Fig. 5. Intensity / vs. time /, in units of T p , for values of 
x T /T p = I, 3/2,4/3, and 5/4. 
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Fig. 6. (a) The relative intensity / of the output field vs. time, in 
units of the period T p , for n w = 20 and x r /T p = 1.01. (b) The 
instantaneous frequency vs time, in units of T p . 

pulses, for several choices of i T /T p . As the differ- 
ence between r r /T p and an integer increases, the 
pulses become longer, and the instantaneous fre- 
quency develops a chirp. The chirp rate increases 
until adjacent pulses begin to overlap. Fig. 6 il- 
lustrates the chirp. With further increase of \x r /T p \ 
the overlapping pulses interfere, and a more 
complex structure develops. 

The predicted chirp rate for the frequency chirp 
shown in Fig. 6 can be estimated from the fol- 
lowing argument. The spectral width of the short 
pulses does not depend on the difference between 
round trip time t r and the time interval T p — 2n/A. 
Thus, for x r = T p we have short pulses with dura- 
tion T p /n m and with the spectral width n^A. If x r is 
not equal to T p (as it is in Fig. 6) the duration of a 
pulse is about T r , but the spectral width is the 
same, n w A. This width is due entirely to chirp, and 
so one can estimate the chirp rate as n^Ajx^ 

6. The moving comb model 

The preceding section presented a description of 
the FSF laser as a static comb of frequency com- 



ponents. An alternative viewpoint of the FSF laser 
output is that in stationary operation the output 
electric field E MC (t) consists of a moving comb of 
chirped frequencies [14,18,29,36,39]. There are an 
infinite number of frequency components ("teeth") 
in this comb, each with a time varying amplitude 

«.(0 



£mc(0 = S 0 B n (t)exp[-i Pn {t)]. 



(76) 



Associated with the real- valued amplitude B„(t) is 
a phase p n {t), whose time derivative is the instan- 
taneous frequency of the nth component 



^(0=^(0- 



(77) 



The moving comb model assumes that each fre- 
quency co n (t) increases linearly with time; for the 
nth component one can write the time-dependent 
frequency as 



0) n (t) = CO max + — (' ~ «T S ). 



(78) 



That is, starting from a very large negative value in 
the remote past (as t — ► — oo), this chirped fre- 
quency increases steadily, exceeding the value co max 
when t = nx s . The frequency continues to increase 
indefinitely thereafter, while the amplitude B n (t) is 
vanishingly small in the remote past, when its in- 
stantaneous frequency lies outside the range for 
which the gain exceeds the loss. As time progresses 
the chirp brings the frequency into the frequency 
domain of net amplification, and the amplitude 
grows. The amplitude reaches a maximum, by 
definition, when t = ht s ; the frequency (o„(t) then 
is a> max . As time increases the instantaneous fre- 
quency grows beyond the region of net gain, and 
losses cause the amplitude to diminish. In the re- 
mote future, when the instantaneous frequency has 
become very large, the amplitude again becomes 
negligible. 

It is generally assumed, e.g. [32], that the am- 
plitude B„(t) has a Gaussian form 



B n {t) = exp 



(n 0 T T ) 2 



(79) 
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peaking at / = nx % and with a Gaussian width pa- 
rameter of n 0 T T . We take the phase /?„(/) to be 



Pn{t) =27^- "T S ) 2 + ">max(' - 



(80) 



From Eq. (80) we determine the instantaneous 
frequency to be equal to the value given by Eq. 
(78). Combining the amplitude with the phase, and 
summing over all possible comb components, we 
obtain the field 



<?mc(0 =<*o £ exp 



(t - *T S ) 2 
("0T r ) 2 



- i — (f - /rr s ) 2 - iaw(f - it s ) 



(81) 



This is the fundamental expression for the field in 
the moving comb picture, when individual ampli- 
tudes are idealized as Gaussians. 



6.1. Interpretation 

It is customary to depict this moving-comb field 
by presenting a density plot of the (real) amplitude 
B rt (t) as a function of time t (along the horizontal 
axis) and of instantaneous frequency a> (along the 
vertical axis). Such a plot appears as a set of 
slanted ridges, each corresponding to one of the 
frequency components. Fig. 7 illustrates a portion 
of such a comb. The intensity distribution shown 
in the picture has the following properties. 

• For fixed frequency (a horizontal line), one finds 
a set of peaks, spaced evenly with time interval 
r s = 2rc/J, the repetition period of the field. 

• At a fixed time (a vertical line), one finds a set of 
ridges separated in frequency by the mode spac- 
ing, 2u/t t . This set of ridges forms a comb in 
frequency space. As time increases, the peaks 
shift toward higher frequency. 

• The nth ridge peaks at time t — /rr s , at which 
moment the instantaneous frequency is co max . 
The ridges all become negligibly small as the in- 
stantaneous frequency becomes very different 
from a> max - 

• Each ridge has the same slope, (d/df)co„(f), 
equal to the frequency-chirp rate y c = A/r r . 




Fig. 7. Schematic diagram of a portion of intensities of FSF 
laser output as a function of frequency co and time t. At any 
fixed time there are a finite number of equidistant frequency 
components separated by the free spectral range, 2ji/r r . At fixed 
frequency there is an infinite succession of pulses, with period 
2n/A. 



This model, of a moving comb of frequencies, 
explains a number of observed properties of FSF 
lasers. However, it has been questioned by some 
authors. In particular, it seems puzzling that, al- 
though the frequency shift occurs only within a 
single short segment of the ring cavity (namely 
the length of the AOM), the output field shows 
only a smoothly increasing instantaneous fre- 
quency [33]. 

The model has been used to describe the 
output of an unseeded FSF laser [37], that is, one 
in which the internal field originates with ampli- 
fied spontaneous emission. Presumably a set of 
discrete-frequency components such as shown by 
a fixed-time slice of Fig. 7 derive from a single 
spontaneous emission event. These are thereaf- 
ter phase coherent with this initial phase. But 
many spontaneous emission events occur, each 
at random times, and so it is not obvious, for 
example, what sets the time to at which any 
component has a specified reference frequency. 
It is likely that frequency feedback may offer 
a similar means of stabilizing a comb pattern, 
by selecting one discrete comb out of a contin- 
uum of combs that differ only in the time t 0 . Any 
irregularities or noise will introduce a departure 
from the regular pattern depicted in Fig. 7. A 
seeded FSF laser deals with this ambiguity by 
providing a definite phase and frequency of the 
input signal. 
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6.2. Unification of models 

In the model presented in Section 5 there is no 
continuously varying frequency, no comb appar- 
ently steadily advancing with time. Instead, the 
frequency takes a succession of constant values, 
being shifted by A with each round trip. The 
output field £ D f(0» where DF refers to the dis- 
crete-frequency model, is the superposition of the 
original seed field and fields that have made all 
possible round trips within the cavity. 

A glance at Eq. (81) shows that the field asso- 
ciated with the moving-comb model is unchanged 
when / is replaced by t + mr s , where m is any in- 
teger. This field therefore shares with the field of 
the discrete-frequency model a periodicity, of pe- 
riod t s = 2n/A. We therefore know that this field 
can be represented as a Fourier series. We write 
this series representation of the moving-comb 
model in the form 



(82) 



where f M c is a reference time. The Fourier coeffi- 
cients C m are determined by integrating over one 
period 

(83) 



A /" ^,+t, 
#oC m = — £ dtE MC (t). 



For the moving-comb model defined by Eq. (81) 
the Fourier components can be evaluated in closed 
form. They are 



C m = 



2n 



V /(/2 0 Tr)" 2 +i(^/2T r ) 



x exp 



(tOmax ~ mA) 



4((« 0 T r )- 2 + i(A/2x t )) 



- imAt MC 



(84) 



We are interested in the case where there are many 
Fourier components, which occurs when no » 1. 
We also assume that the frequency shift A is 
comparable to the cavity mode separation, mean- 
ing that Ax r /2n is comparable with unity. In this 
case, with the choice of reference time as 
*mc = oj max T r /A, the Fourier components can be 
written as 



C m = C 0 exp 



[-(=^2)1 



x exp[+im 2 (4T r /2)]. 



(85) 



Here the factor C 0 incorporates the common 
phase, the first exponential describes a Gaussian 
distribution of amplitudes, and the final expo- 
nential describes the (quadratic) dependence of 
phase upon m. 

Although the discrete frequency model of the 
previous section deals with a seeded laser, there is 
no such reference frequency for the moving-comb 
model. Nevertheless, the two pictures have the 
same Fourier transforms, and hence they describe 
the same output field, if only one assumes that the 
frequency of the seed laser of the discrete-fre- 
quency model is some large integer multiple of A. 
The connection between the two models is com- 
pleted with the identification of the Gaussian 
widths, « 0 = n w - 

6.3. Comments 

Our analytical expressions of the earlier sec- 
tions provide a first-principles derivation of the 
width and the peak of the discrete-frequency 
model and, from the proven equivalence of the 
two models, the moving-comb model as well. 
The quadratic /n-dependence of the phases of the 
Fourier components of the moving-comb model, 
C m , given in Eq. (85) are identical with the phases 
of the amplitudes A m of the discrete-frequency 
model. With simple assumptions about the gain 
and loss properties of the device, one obtains in 
both cases Gaussian distributions of amplitudes. 
Thus the two models are identical when they both 
describe a FSF laser seeded by a CW field. Then 
Eq. (85) is exactly the field of a seeded FSF laser 
whose seed frequency is equal to some (very 
large) integer multiple of A. 

The equivalence of the two descriptions means, 
in particular, that even though the moving comb 
model appears to have a continuously varying 
linearly chirped frequency, in fact the field con- 
tains only discrete frequencies when it is seeded. 

The moving comb model has been applied not 
only to a seeded system but to a system that grows 
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from spontaneous emission. In this case there is no 
identifiable seed carrier frequency co s . Neverthe- 
less, the system retains the same periodicity es- 
tablished by the presence of a fixed frequency shift 
A at regular intervals. 

In all cases only these discrete frequencies nA 
are present in the field. This means that unless 
the field passes through some nonlinear element, 
then it cannot acquire any additional frequencies. 
This implies that for use in interferometry that 
depends on measuring a beat frequency whose 
value is proportional to delay time between two 
paths [28,29,36] such a field can not be used 
directly. 



7. The FSF laser with spontaneous emission 

The procedure used to solve Eq. (23) is also 
applicable to the treatment of a FSF laser based 
on the growth of spontaneous emission - the laser 
acts as a regenerative amplifier of spontaneous 
emission. The relevant equation is 

g{X) = $(X - X 0 ) exp [g(X)] + £{X). (86) 

The solution to Eq. (86) is 



S{X) = Z{X) + £ exp 



f^G(X-lX Q ) 



1=0 



x £(X - nX 0 ). 



(87) 



We now use the properties of our assumed noise 
model to evaluate stochastic averages. 

7.7. Stochastic averages 

We are interested in the intensity /(w, T) of Eqs. 
(18) and (19) averaged over stochastic realizations 
of the Langevin force £(X) 

J(w, T) = (/(a?, T)). (88) 

For the steady-state cases considered in most of 
this paper, this average is independent of T; in 
Section 8 we consider more general situations. To 
carry out the stochastic average (•• •) and obtain 
the spectrum J(w, T), we assume that the sponta- 
neous emission process £ sp (0 is a delta- correlated 
(Wiener-Levy) stochastic process 



<MM£ P ('2)> = - '2). (89) 

The parameter £q that characterizes the sponta- 
neous emission intensity depends on the inversion 
w, the spontaneous emission rate y s , and the ge- 
ometry of the laser cavity. Because the width of the 
filter is much larger than the free spectral range 
/> » l/t r , the result of this averaging is the 
spectrum 



1 f w-lA \ : 

2 1 r t ) 



n-\ 
/=0 



gmW - f m 



(90) 



After integrating Eq. (90) over w we obtain the 
averaged intensity 



/-f-oo 
dm J(m, T) 
CO 

r t 2 oo 

£ 

n—\ 



16tt 3 /2 

-A 



n 12 



(91) 



7.2. Simplification 

We assume, as before, that zl < f f . In this case 
we can replace the summation over n in (91) by 
integration and write 

7 = VW, 

A 



[2[ gm w-f m } + A 2 /(\2r 2 t )}-^ 2 , m 



Vnr r 

where the function F(x) is defined as 

F(x)=x i/3 / C "diexp[r 2 (l -x 2 / 4 )] (93) 
Jo 

and the magnitude of the spontaneous emission 
source is characterized by the parameter 



SP " Vrt \ 2" A ) 



1/3 



(94) 



By taking into account that g m w~ f m one can 
transform Eq. (22) for vv into 
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p sp F(X) = l, 



(95) 



where we have introduced the (small) parameter 



ftp = 



•sp 



U(n -l) 



«: l. 



(96) 



Under typical operating conditions the parameter 
/? sp is very small: /? sp ~ 10" 18 . Then we can use the 
asymptotic form of the function F(x) for large 
x » 1 



K 



1/2 



jc 1/3 exp 



(97) 



In this case the initial inversion w 0 can be found 
from the condition 



^-^exp(y ) = y«/3 ) 



31/222/3 

where 



(98) 



y =wx=jj - /«) + A 2 mm V2 - 

The solution to Eq. (98) for fi sp <C 1, and thus for 
.ln(l/&p) » 1, is 



y ^ ln 7T- 

Psp 

Then the requirement for inversion is 



2(g m w -f m ) = - 



nr 2 r 



3J 1 



1 2/3 



(99) 



(100) 



Taking into account the condition J T f we find 
that g m w ~ / m , as is expected. The spectrum 
J(m, T) can be written now as 



7(n7 ' r) = i^°£ exp 
j 2 



i2r 2 




(101) 



In Eq. (101) one can replace the summation over n 
by integration. Taking into account the condition 
ln(l//? sp ) » 1 we can carry out the integration, 
with the result that the spectrum takes the form of 
a Gaussian 



_ r . 1 (nj — <5a>) 
y(m,r) = 1^-1)—— exp 



whose width T is 

r = y/V^ 

and whose peak is offset from m = 0 by 

<5cu = cry 

with 



=(K) 



1/3 



(102) 
(103) 
(104) 

(105) 



Recall that m is the frequency offset from the 
filter center cof, meaning that the spectral maxi- 
mum occurs at aj max = cor -f ay. This result for the 
continuous spectrum is to be compared with that 
of Section 5.1 for the distribution of discrete 
frequency components within a Gaussian enve- 
lope. 

As in the case of the seeded laser we have here a 
first-principles derivation of the width and the 
peak of the spectrum. Again there is only a weak 
dependence on fi sp and the, only important pa- 
rameter is y = [Ar*] ]/3 . It is noteworthy that pa- 
rameter a can be easily determined from the ratio 
dco/t = cr 3/2 . In principle these are known for a 
particular laser, and so the operation of the FSF 
laser can be determined from these known quan- 
tities [24]. 



8. Rate equations 

A previous paper [24] has used a rate-equation 
model to describe the operation of the FSF laser. 
This model is based on two sets of coupled equa- 
tion, one for the population inversion and another 
for the spectral density of photons. It is interesting 
to compare our approach based on the Fourier 
expansion of the electric field with the rate-equa- 
tion model. 

Considering the FSF laser with spontaneous 
emission and using Eq. (86) for S{X) in the model 
of spontaneous emission, Eq. (89), we obtain an 
equation for the averaged spectral density J (ro, T): 



LP. Yatsenko et ai I Optics Communications 236 (2004) 183-202 



201 



J(w, T) =J(w - A,T -x r ) 

x exp(2[^o(n7)w(r) JjL. 

(106) 

This equation was derived earlier using the con- 
dition of broad loss T f » l/r r . Here we do not 
assume steady-state operation and we allow the 
averaged spectral density J(w, T) to be time de- 
pendent (on the time scale of inversion relaxation). 
In traditional treatments of laser dynamics one 
assumes changes of J(w, T) are small during one 
round trip, J{xd, T) ~ J(w — A,T — x r ). In this 
limiting case Eq. (106) can be written as 

= J(w,TK2[g 0 (w)w(T) -f(w)}) +^^- 

(107) 

With the same assumptions the differential equa- 
tion for the inversion reads 

(108) 

Here y s , as defined in Section 3.2, is the sponta- 
neous emission rate and fjj, defined in Eq. (89), 
parameterizes the spontaneous emission intensity. 
Apart from notational changes, these equations 
differ from the earlier rate equations of [24] only in 
the presence of the partial frequency derivative in 
Eq. (107) for the spectral density. In [24] the au- 
thors discretized the frequency space. For a 
smooth spectrum, which we assumed here, this 
discretization replaces Eqs. (107) and (108) with 
partial difference equations, leading to the ap- 
proximate solutions discussed in [24]. 

9. Summary and conclusions 

We have presented a formalism for treating the 
coupled equations linking electric field output and 
population inversion for a frequency-shifted feed- 
back laser. The equations take full account of the 
phase of the field, and include a possible seed laser, 



which may have arbitrary phase modulation or 
fluctuation. 

Using these equations we have derived simple 
expressions for the steady-state output electric 
field. These have a simple interpretation in terms 
of multiple passes of a wavepacket through the 
system, gaining or losing energy while undergoing 
a periodic frequency shift A. 

We show that the distribution of spectral com- 
ponents follows a Gaussian form, as is often as- 
sumed. The width of this distribution determines 
the minimum duration of pulses. The characteris- 
tic parameters of the Gaussian amplitudes are al- 
most entirely determined by the frequency shift A 
and the filter bandwidth T f , through the combi- 
nation y = (zir f ) ,/3 . 

The form of the output varies with the relative 
phases of these spectral components, and depends 
critically on the ratio of the cavity round-trip time 
T r = L/c to the repetition period T p = 2n/A. We 
have presented illustrative examples of the steady- 
state output in which there occur various forms of 
short pulses. Under appropriate conditions these 
pulses exhibit a linear frequency chirped. 

We have also shown the equivalence of two 
popular models for the output field of the FSF 
laser, one in which there is a moving comb of 
frequencies and the other in which there is a fixed 
discrete set. Our analytical expressions provide a 
derivation from first principles of the width and 
the peak of the spectrum of the moving comb 
model. We also show that the rate equations used 
previously in modeling the FSF laser can be re- 
covered from the equations presented in this 
paper. 
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Abstract 

The potential advantages of chirped pulses for very precise measurement of distance, through frequency-domain 
ranging, has prompted consideration of frequency shifted feedback (FSF) lasers as sources of interferometer light. 
We here derive theoretical limitations to the spatial accuracy one can expect in such applications, by considering ana- 
lytical expressions for the electric field emerging from a frequency shifted feedback (FSF) laser seeded by a CW laser 
whose finite bandwidth originates in phase fluctuations. We also consider consequences of fluctuations in cavity size. 
We show that, for surfaces flat within the laser footprint, such a system can provide the subwavelength accuracy of 
conventional interferometry but without dependence on material-dependent phase shifts. Although noise has been 
important for previous uses of FSF lasers in optical ranging and interferometry, we here show that a frequency mod- 
ulated seeding laser can be used to better advantage than noise. 
© 2004 Elsevier B.V. All rights reserved. 

PACS: 42.55.-f; 42.60. Da; 42.55. Ah 

Keywords: Optics; Lasers; Frequency shifted feedback; Interferometry; Profilometry; Frequency-domain ranging 



1. Introduction 

Rapid and nonintrusive measurements of dis- 
tances from millimeter to kilometer with accura- 



* Corresponding author. Present address: Department of 
Physics, 618 Escondido Cir, Livermore, CA 94550, USA. Tel.: 
+ 1 925 455 0627. 

E-mail address: bwshore@alum.mit.edu (B.W. Shore). 



cies of microns or less are now feasible using 
techniques of optical frequency domain ranging 
(OFDR) wherein a measurement of distance is ob- 
tained from a measurement of frequency differ- 
ences [1-5]. The technique can be regarded either 
as a form of chirped optical radar [6] or as a form 
of interferometry using chirped laser input. As has 
been demonstrated, [7,8] the needed interferometer 
seed laser can use the technique of frequency- 
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shifted feedback (FSF) [9-17]. In such a device 
feedback introduces a fixed frequency shift A with 
each passage of a wave packet around the cavity 
(during the round-trip time x r = Lie in a cavity of 
perimeter L). By interpreting the output as a mov- 
ing comb of frequencies [18,19,15,20-22] Nakam- 
ura et al. [7,8] analyzed the use of FSF for 
distance measurements. A key parameter of their 
approach is the chirp rate y c = Alx T : the difference 
in lengths of the two interferometer arms h can 
be found from the measured beat frequency a) B 
as h = caW2y c . 

An alternative view is also possible. Under 
proper conditions the output of a seeded FSF laser 
can be a regular train of pulses [9,10,18,23,13,24] 
and the individual pulses can exhibit a frequency 
chirp [26]. The apparent similarity with the longer- 
wavelength pulses used in chirped radar for accu- 
rate ranging suggests that operation at optical 
wavelengths could offer accuracy adequate for 
nanoscale structures. However, this simplistic view 
of FSF laser output is misleading. The chirp of 
individual pulses from a FSF laser is not con- 
nected with the average chirp rate y c but, as was 
shown in [26], is instead determined by the devia- 
tion of the frequency increment A from cIL. 

Although the literature on FSF lasers is exten- 
sive, there remain some fundamental questions 
[25]. Some of these concern the ultimate limits on 
accuracy attainable when such a laser is used in 
interferometry. Earlier we have developed a first- 
principles treatment of the FSF laser [26]. In the 
present paper, we reexamine some of the proper- 
ties of this device, with particular emphasis on 
interferometric applications. Our results provide 
useful estimates of the limiting accuracy one can 
achieve using a FSF laser for distance measure- 
ment. We clarify the importance of noise in pro- 
viding the beat frequency that is used in the 
frequency-domain ranging, and we suggest an 
alternative approach, using controlled frequency 
modulation rather than uncontrolled noise, to pro- 
duce the desired signal. 

The present paper is organized as follows. 

In Sections 1.1 and 1.3, we present a simple dis- 
cussion of relevant aspects of interferometry and 
of frequency-domain ranging based on chirped 
pulses. 



In Section 2, we discuss the properties of the 
FSF laser, and the form of the output field. The 
analytic expressions given here provide the basis 
for our discussion of FSF interferometry that 
follows. 

In Section 3, we consider interferometry using 
the field presented in Section 2. We point out the 
importance of phase fluctuations in providing the 
frequencies needed for range measurement. 

In Section 4, we consider the FSF laser seeded 
by a phase-noisy cw input field. We derive expres- 
sions for the beat-frequency spectrum. The width 
of this spectrum determines the accuracy with 
which the peak position (the beat frequency) can 
be measured, and hence the accuracy with which 
one can measure the distance. 

In Section 5, we consider the effects of unavoid- 
able fluctuations in the cavity length, and the 
effects these have on the accuracy of beat- 
frequency measurement. 

In Section 6, we consider using a frequency- 
modulated seed to the FSF laser. We show that 
this has advantages for improved accuracy. We 
consider both phase sensitive and phase insensitive 
detection of the interference intensity. 

Section 7 discusses the accuracy obtainable, 
using either a noisy seed or deliberate frequency 
modulation to provide the needed frequencies. 

Section 8 summarizes our results. 

An appendix provide details of our noise model. 

LI. Interferometry 

Interferometry offers a means for obtaining ex- 
tremely accurate distance measurements, such as 
are needed for determining the heights of nanom- 
eter-scale features on lithographed surfaces [27], 
In a typical Michelson implementation, a laser 
beam passes through a beam splitter from which 
one beam travels a fixed distance in a fixed refer- 
ence arm. The other beam travels to a sample, 
where it is reflected. The field of the returning 
beam is combined with that of the reference beam 
and the resulting intensity pattern is recorded; 
changes of intensity provide a measure of phase 
difference between the two beams: Apart from a 
fixed (but important) phase increment upon reflec- 
tion, the phase is proportional to the difference in 
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lengths of the two paths. Thus a measurement of 
intensity, converted to a measurement of phase, 
gives a measurement of distance. 

One standing difficulty with interferometric 
techniques applied to surfaces comprising several 
materials, often of unknown composition, is the 
presence of a material-dependent phase shift of 
the reflected field. For reliable measurements of 
surface heights to tolerances of a few nanometers, 
it is essential that this material phase shift be 
accurately known [28]. This shift is particularly 
difficult to determine when the surface comprises 
granules embedded in a substrate, unless the 
granules are either much smaller than a wave- 
length, so that the surface can be characterized 
by a complex index of refraction for an effective 
medium, or are much larger than a wavelength, 
so that they can be resolved readily. Even then 
it is difficult to evaluate the phase shift. A recent 
paper offers one means of overcoming this limita- 
tion [27]. As we will show, chirped pulse interfer- 
ometry also offers the potential to overcome these 
handicaps. 

1.2. Ranging 

Conventional interferometry with a single 
wavelength X provides subwavelength accuracy of 
path differences, but because phase has an intrinsic 
periodicity, it can only provide this accuracy mod- 
ulo A, i.e., it cannot distinguish between distance h 
and h + A. Various techniques are used to over- 
come this limitation, such as the use of two differ- 
ent wavelengths. 

For measurements of distances of meters or 
miles, radar techniques are often used. Typically 
these transmit a train of pulses and measure the 
time delay before a return pulse is observed. The 
requirement for very precise measurement of time 
intervals makes it difficult to obtain positions of 
even stationary objects with micron accuracy. 
The technique of chirped radar [29] replaces time 
measurements with beat-frequency measurements, 
but has much in common with interferometry. The 
technique described here can be used for measure- 
ments of large distances, as is done with optical 
ranging, as well as for small distances; all distances 
are measured with equal accuracy. 



1.3. Chirped-pulse interferometry 

It is easy to show that if one has a chirped fre- 
quency as input to a Michelson interferometer, see 
Fig. 1, then the Fourier transform (i.e., the spec- 
trum) of the beat intensity includes a frequency 
equal to the product of the chirp rate and the delay 
time between the two arms. Thus one can antici- 
pate construction of an interferometer in which a 
measurement of frequency gives a measurement 
of distance. To show this connection consider a 
simple model of a field whose carrier frequency 
varies linearly with time (a chirped frequency) at 
rate y c during the time interval 0 < / < T 0 . Assume 
a constant amplitude &. Then the complex-valued 
field at a fixed position is 

£•(0 = (fexp[-i(a> 0 + yj/2)t). (1) 

We use this field as input to a Michelson interfer- 
ometer in which one arm has length z and the 
other has length z + A, where h is to be determined. 
Fig. 1 illustrates the general layout of the interfer- 
ometer. Our interest here is in sub-micron metrol- 
ogy of surfaces, such as occur with lithography, 
wherein the reference beam can be reflected from 
the substrate that holds the sample; in other work 
the arms are often of very different lengths [7,8]. 

Upon reflection the electric field of beam j ac- 
quires a phase shift <f>j 9 which depends on the prop- 
erties of the surface. The difference in arm lengths 
gives a relative time delay 

T = 2h/c. (2) 




Laser source 



Reference arm 



Beam splitter 



i 



Sample 



Fig. 1. Schematic diagram of interferometer. Laser source 
passes through beam splitter, and then along two arms, one of 
which contains the unknown height to me measured. The height 
h is relative to the reference arm length, as shown. 
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For simplicity, we assume perfect reflection at the 
reference and target surfaces and T < T 0 . Then the 
interferometer output intensity is 

/(/, T) = \E(t)+cxp(\A<t>)E(t + T)\ 2 

= \B(t)\ 2 + \E{t+T)\ 2 

+ 2Re[exp(iA0)£*£(r + T)]. (3) 

where A(f> is the difference between the material-de- 
pendent phase shifts 4>j in the two arms. We are 
interested in the interference-term signal 

5(0 = 2Re[exp(iAtf>)E*£(r + T)) 

= 2\S\ 2 cos[(2Ay c /c)f + A0 + 9{T)\. (4) 

This expression shows that a measurement of the 
beat frequency [the oscillation frequency of S(t)] 

co B = 2nv B = y c T = {yjc)2h (5) 

converts, using the chirp rate y Cy directly into a 
measure of T and thence to a measurement of dis- 
tance h. The material-dependent phase difference 
A<p adds a constant to the geometrical phase 
increment 

<P(T)=to 0 T-y c T 2 /2 (6) 

to produce the interferometer phase (P(/) + A</>. 
Thus any determination of T from the phase of 
the interference term must rely on values of the 
material phase shift A<f>. By contrast, measure- 
ments based on beat frequency determination are 
insensitive to this phase. 

As will be explained, interferometry with the 
FSF laser also offers a means of avoiding depend- 
ence of distance measurements upon the material 
phase shift A<£, and hence it offers the potential 
for accurate surface topography measurements 
without knowledge of surface characteristics. Fur- 
thermore, by giving a measurement of heights of 
both reflecting and partially reflecting surfaces, it 
can reveal structure underlying a partially trans- 
parent layer with known refractive index but of 
arbitrary thickness. 

The FSF laser is particularly suitable for optical 
ranging: it provides a strictly linear chirp that can 
have a very large chirp rate, say 100 MHz in 5 ns, 
or 5 x 10 17 s" 2 . The FSF laser also functions with- 
out an external seed, driven by spontaneous emis- 
sion in the gain medium. The work of Nakamura 



et al. [7,8] is done with such a laser. Our work dem- 
onstrates that, when the properties of the laser are 
dominated by an external seed, then much im- 
proved accuracy of distance measurements can 
be achieved. The bandwidth of the seeding laser 
determines the accuracy, but the origin of the 
bandwidth is also significant. We show that band- 
width from controlled phase modulation is partic- 
ularly useful. 

L4. Accuracy and resolution 

When the laser beam covers a surface pattern 
having several heights, then a histogram of distance 
values displays several peaks. To resolve these 
using frequency chirping the difference in heights 
must be greater than half a wavelength. This limita- 
tion can be understood in the following way. 

From the relationship (5) between the measured 
feature height h and beat frequency to B it follows 
that the height error Sh is proportional to the 
beat-frequency error, Sh = (c/2y c )Sco B . This error 
is approximated as dco B « 2n/T 9 where T is the 
duration of a pulse. During a pulse the frequency 
changes by y c T. This change cannot be larger than 
the optical frequency to = 2nclX 9 and hence the fre- 
quency error cannot be less than dco « (yjc)k. 
From this limit it follows that the resolution limit 
is set by the optical wavelength A, 

Sh > X/2. (7) 

When the laser beam covers a flat surface then, as 
we will note, the error in a measurement of surface 
height relative to a reference surface can be signif- 
icantly less than a wavelength. This error, the 
accuracy of interferometer measurement, is deter- 
mined by the accuracy with which one can meas- 
ure the location of the peak value of a spectral 
profile, i.e., by the error 5co B in the beat frequency 
co B . In turn, this is set by the resonance width r rcs 
of the spectral profile. As we will note (see section 
7), accuracy of 100 nm appears possible. 



2. The FSF laser 

In brief, a FSF laser contains a feedback cavity 
(which is typically either a closed loop, of perime- 
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ter L, or a Fabry-Perot cavity of length L/2) in 
which there is a section having gain, one or more 
spatial or frequency filters and, most importantly, 
an element that induces a frequency shift (of angu- 
lar frequency A) on each wavepacket that passes 
through it. Typically the frequency shift is pro- 
vided by an acousto-optic modulator (AOM), act- 
ing as a grating that feeds the first order diffraction 
back into the cavity, while the zero-order light 
emerges as the output of the ring cavity. Fig. 2 
illustrates an example of such a ring cavity [25]. 

The properties of such a laser have been exten- 
sively reported, e.g., [30-37,13,24,14,38]. The 
important parameters of the FSF cavity itself, 
apart from gain and loss, are as follows: 

A = the frequency shift (rad/s) per round trip; 
r r = Lie = the round trip time in the cavity; the free 
spectral range is 2n/x r . 

t s = 2n/A = the repetition period of the field; 
y c = A/x r = the average chirp rate of the instanta- 
neous frequency of the cavity field. 

Another set of parameters describe the net gain 
within the cavity and the loss due to reflection and 
filters. These are distributed over some finite range 
of frequencies. We treat them as localized in an 
infinitesimally thin segment of the path. Fig. 3 
shows the distribution of gain and loss as a func- 
tion of frequency (discretized by ri). 

A key parameter controlling the beat spectrum, 
from which follows the achievable accuracy of 
beat frequency measurement, is 

n w = the effective number of discrete frequency 
components within the output spectrum from the 
seeded FSF laser. (The output bandwidth divided 
by the the frequency shift). 

AOM 

seed _J^»— _ output 





gain 









Fig. 2. Schematic diagram of ring cavity showing mirrors (M), 
frequency shifting AOM (FS), beam splitter (BS), amplifier (G), 
and input and output beams. 
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Fig. 3. (a) The relative effective gain function Q(n)/Q,„ vs. n. (b) 
The magnitudes a„ of the field components £,,(0 vs. n, 
normalized to unity at the maximum. The parameter 1/ 
(n„) 2 = 240. The frequency range of positive net gain is 2<p 0 , 
and Q m is the maximum of Q{n). For, clarity only every second 
component is plotted. One can see, from comparing the two 
frames, that the peak of the component distribution occurs 
where the effective gain changes from positive to negative. From 
[26] Fig. 2. 

T f = the bandwidth of the frequency limiting ele- 
ments of the cavity (a filter). 
y = (zJ 2 r f ) 1/3 = a frequency parameter largely 
determining the properties of the FSF laser 
output. 

A FSF laser can operate in a variety of regimes: 
the output may range from extreme irregularity to 
very regular pulse trains or even cw, depending on 
such controllable properties as the amount of gain 
and the relative values of various time constants or 
frequency bandwidths of the device. We are here 
concerned only with steady-state operation (i.e., 
infinite pulse trains or cw). 

2.7. The moving-comb model of FSF output 

One viewpoint of the FSF laser output, used 
extensively by Nakamura et al. [15,16], is that in 
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stationary operation the output electric field 
£* MC (/) consists of a moving comb of chirped fre- 
quencies. There are an infinite number of fre- 
quency "teeth" in this comb, each with a time 
varying amplitude B n (t): 



E MC (t) = S 0 B a (t)cxp[-ip„(t)]. 



(8) 



Associated with the real-valued amplitude B n (i) of 
tooth n is a phase p„(t) whose time derivative, the 
instantaneous frequency (o nith varies linearly with 
time (a frequency chirp) 



co„(t) = ^P.(0 = *W +^*(*-«t s ). 



(9) 



Typically the amplitude B„(t) is taken to be of 
Gaussian form, 



B„(t) = exp 



(t - nx s f 



(n«.T r ) 



(10) 



peaking at t = nx s and with a Gaussian width 
parameter n w x T . 

Because the output of the FSF laser can be 
viewed as a comb of chirped frequencies, it is nat- 
ural to consider using it as a source of radiation 
for chirped-pulse interferometry. However, this 
use is not as direct as one might think. In reference 
[26] we showed that this output of the FSF laser 
can also be written in the form of a Floquet series 

£mc(0 = &oY, C " QX Pl-im A (t - *Mc)h (11) 

m 

where t M c is a reference time. We take 
'mc - o) max r J A and n w » 1 and thereby write the 
Fourier coefficients as 



C m = C 0 exp 



(m - (co ma jA)) 2 



exp 



+\m 



Atj. 



(12) 

Here the factor C 0 incorporates the common 
phase, the first exponential describes a Gaussian 
distribution of amplitudes with width n w> and the 
final exponential describes the (quadratic) depend- 
ence of phase upon m. This decomposition shows 
that the field comprises a discrete set of time-inde- 
pendent frequencies, separated by the frequency 



shift A. By contrast, the field needed for interfer- 
ometry should comprise a continuum of frequen- 
cies if it is to include the beat frequency 2y c T for 
arbitrary delay time T. 

In [26] we have shown that, because of the dis- 
creteness of the frequencies actually present in the 
moving comb model, it can be viewed as a discrete- 
frequency model (our term) seeded by a cw laser 
whose frequency is some integer multiple of the 
frequency shift A. 



2.2. FSF cavity seeded with variable phase 

The frequencies required for frequency-domain 
interferometry can be supplied in several ways. To 
observe a beat signal at some frequency it is essen- 
tial that the seed laser spectrum contain this fre- 
quency. Noise fluctuations provide one potential 
source of the desired bandwidth. Previous authors 
[15,16], have relied on random amplitude fluctua- 
tions to provide seed bandwidth. Here we discuss 
the use of a modulated seed pulse, of carrier fre- 
quency co s . When we allow a variable phase <p s (0 
on this input to the FSF cavity, the output field 
has the form 



E{t) = J E n(t) exp [-i(w s + nA)t], 



(13) 



n=0 



where the complex-valued time-dependent ampli- 
tude is 

E n {t) = fic^ exp [-\<P n - iq> s (t - «T r )]. (14) 

The factor e c provides a normalization: it is related 
to the seeding-laser intensity inside the cavity I c 
through the formula 



(15) 



The dimensionless amplitude a n is, in the limit 
of interest here (large n„, see [26]), 

2" 



a n ~ stf exp 




(16) 



where amplitude sf, width n w and central position 
n m are determined by laser parameters (see [26]). 
Here we assume that n w ^> 1 . The phase has a time- 
independent contribution that is quadratic in n, 
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*, = - Tr ^ + (fl+lW2l. (17) 

The field of Eq. (13) also receives a contribution 
from time varying phases q> s (t), either random or 
deterministic. 



3. Interferometry with a FSF laser 

We are interested in the homodyne beat signal 
created by the interference between two electric 
field amplitudes with relative time delay T- 2hlc 
and with a phase difference A<f> originating in the 
difference between the material dependent phase 
shifts of reflected light in the two arms. This signal 
is proportional to the real part of the complex-val- 
ued product Qxp(iA<t>)E(t)£r(t + T): 

S(t) = exp(iA<f>)E{t)E*(t + T) + c.c. 

= ^2 E„(t)E',(t + T) exp [iA<£ - i(a> s + nA)t 

n,/=0 

+ i(co s + lA){t + T)] + c.c. 
= exp[iA<£ + \<o s T\ £ £ E n {t)E'„_ m {t + T) 

n=0 m=-oo 

x exp [\(n - m)AT - imAt) + c.c. (18) 

(Note that we need not consider the spectrum of 
E(t)E*(t) or E(t + T)E*(t + T) because these con- 
tain no information about the distance.) This sig- 
nal contains many harmonics m, originating with 
combinations of the various frequencies in the 
sum. We will assume that the frequency 
7> c = TAIx T produced during the time T is much 
less than the frequency shift A. This implies that 
T <§: r r and, in turn, that the possible beating fre- 
quency is small. For application to surface profil- 
ometry, where one is interested in measuring 
small variations in height, it is always possible, in 
principle, to adjust the reference arm relative to 
the sample so that T «: r r . For use with large dis- 
tances, in optical ranging, this is not always possi- 
ble and it is necessary to resolve the phase 
ambiguity of successive chirped pulses. That is, 
the time T is determined by interferometry only 
within some multiple of the repetition time t s . To 
resolve this ambiguity one can use a second 
AOM frequency to determine the order m to be 



used with Eq. (18), as was done by Nakamura 
et al. [8]. 

Without loss of generality we will here consider 
only the terms in (18) for which / = n or m = 0. The 
signal is therefore 

x exp [iAtf> + i(co s + nA)T] + c.c. (19) 

The only time dependence of this signal origi- 
nates in phase variation of the seeding laser, as de- 
scribed by <p s (0 in E n (t). 

Note that in the absence of phase variation, 
<P$(0 — 0, only discrete frequencies contribute to 
the electric field: these differ from the seed fre- 
quency by n increments of the AOM frequency 
shift A. In particular, the interferometer generates 
no new frequencies. The time shift T leads only to 
a phase shift of beating at the frequency A and at 
the higher harmonics. 

We emphasize that when there is no variation of 
the laser phase (e.g., no noise and no deliberate 
modulation), then there is no beat signal. Experi- 
mental observations [14] have confirmed that, 
when feedback in the FSF laser causes mode lock- 
ing (and a consequent removal of noisy phase fluc- 
tuations), then any interferometric beat signal 
disappears. Thus to understand the characteristics 
of the beat signal, and the accuracy attainable in 
interferometry, we must consider the effect of noise 
of the seed laser. As we will show, it is possible to 
produce a beat signal whose power spectrum has a 
very sharply defined maximum at the frequency 
cob = Ty c . We also consider the possibility of intro- 
ducing frequency modulation of the seeding laser, 
rather than relying on noise. We will show that this 
technique offers promise of high accuracy. 

Our task now is to compute the power spectrum 
Js(oj) of the beat signal S(t). This is determined as 
the Fourier transform 

J s {oj) = 2 I dr G(t)coscot 
Jo 

= /°°dTG(T)exp(iorr) (20) 

J -oo 

of the autocorrelation function C7(t), 



588 



G(T)==<S(/)S(f + T)>. (21) 

Here and elsewhere <• • •> denotes a stochastic aver- 
age. Because we deal with steady state, the auto- 
correlation function is independent of time t. 

4. FSF interferometry with a noisy seed 

In the presence of phase variation of the seed la- 
ser the autocorrelation function G(z) can be writ- 
ten as 

G(t) = £ £ MV| 2 |a/|V + (T - /tr + nx r , T) 

11=0 /=0 

x exp [i2A</> + i(2cu s + n A + IA)T] 

+ £ £ \<k\ 4 \an\ 2 \a,\ 2 F-(T - h t + HX r , T) 

n=0 1=0 

x exp [+i(/i - I) AT] + c.c. (22) 
where we define 

F±{t\,t 2 ) = (exp[-i^ s (0+i<P s (' + 'i) 

T i<p s (f + f 2 ) ± i<Ps(' + '2 + 'i )])• (23) 

Note that the only dependence upon the material 
phase shift A<£ occurs with the first summation, 
in the exponential associated with F+. 

We now consider the Fourier transform of this 
correlation function. From the above expressions 
we obtain the formula 

j s (a>) = F + (a>, T)&(a> + y c T)&(-a> + y G T) 
x exp [+iAtf> + i2co s T) + F + (co, T) 
x <?(a> -y c T) &(-<*> -y c T) 
x exp [-iA</> - i2w s r] + 2F_(co, T) 
x £f(co - y c r)jS? (-a> + y c r) (24) 

where all of the dependence upon phase fluctua- 
tion occurs in the factors 
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4.1. The profile function 



F± (co, T) = / dt F ± (t, T) exp [ion] , 
while the remaining factors 
j£f (o) = |fl„| 2 exp [+icoT r *] 



(25) 



(26) 



n=Q 



serve as profile functions independent of any phase 
fluctuation. 



Using formula (16) we write the profile function 
££(to) as 



£?(oj) ~A 2 ^2 ex P [+ iwT r"l exp -2 



n=0 



(27) 

When the number of frequency components be- 
comes large this is a very narrow function of oj. 
In this situation, with n w » 1 and n m » 1 the sum- 
mation in this expression can be replaced by inte- 
gration from -oo to +oo, with the result 

<£(oS) ~ n„A 2 J- exp [+\coT T n m ] exp g— 

(28) 

In this case, we can extract a phase and an ampli- 
tude from the profile function j£? (o>) and introduce 
a normalized Gaussian profile function 

\2~ 



^gH = exp 



(C0T r /1 w ) 



8 



(29) 



Then the spectrum takes the form 

J s (o) = 7r(/i vW 4 2 ) 2 |€ c | 4 [F + (a>, r) cos [2A 0 r + 2zl# 

x^ g (co + y c r)if g (a;-y c r)-hF.(co,r) ■ 

x [j? g (a>-y c r)] 2 ]. (30) 

In this formula, like the general expression of Eq. 
(24), the contributions from F+ involve the prod- 
uct of profile functions J£? g (co) centered about the 
values +y c T and -y c T. The width of 
^fg(co) is \Zl/(w w T r ). When this width is much less 
than y c r, as it is for large n w , then the overlap of 
these functions will be negligible, and the only con- 
tribution to the spectrum J^co) will come from the 
last term. Under these conditions, we obtain the 
formula 

J s (a>) = n(*U 2 ) 2 |ec| 4 /M«, r) - y c T)] 2 - 

(31) 

Because the spectral width of F_(a>,7) is usually 
quite large compared with the width of ££ % {oS), 
Eq. (31) expresses a very narrow resonance, in 
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which position and width depend on laser param- 
eters (through if) and the magnitude is deter- 
mined (through F-) by the noise characteristics. 

4.2. Modeling seed- laser noise 

The function F_(a),T) depends on the fluctua- 
tion spectrum. To evaluate this dependence we as- 
sume that the time derivative of the seed-laser 
phase (p s (t) is a zero-mean stochastic process £(/) 
[see Appendix, Eq. (A.4)], 

m = «0- ( 32 ) 
In Appendix A, we discuss using an exponentially 
correlated (Ornstein-Uhlenbeck) process, with 
zero mean and exponential correlation 

<«> = 0, «(0«(O> = GDexp(-G\t - 

(33) 

to carry out the stochastic averaging in expression 
(23). We use here the results for a limiting case, a 
delta-correlated process G — ► oo, for which simple 
analytic expressions can be obtained. 

As shown in Appendix A, the width of the func- 
tion F-(a>,T) is larger than \IT. The spectrum gi- 
ven by Eq. (31) is appreciable only within a 
narrow interval, of width 6co « y/E/(n w x r ) near 
the resonance value co B = y c T. Because the delay 
is small, T T r , the width obeys (for any n w ) the 
inequality 5a> 1/r, and we can take F_(a>,T) to 
be F-(y c T,T), Moreover, from the conditions 
Ar r « 1 and T<^z T we have y c T<^ 1/r, and so 
we can take F_(co,T)^F_(y c ,T) =F_(0,7). Using 
expression (A. 13) we obtain the result 

F_(0, T) = {l/D)[l - (1 +2D7)exp(-2Dr)]. 

(34) 

Thus the expression (31) can be written as 
Js(o>) = J m exp [-(co - yjf(x t n w /2f\ , (35) 

where 

J ni = n{ nw A 2 )%\\\/D)[\ -(1 +2Z?r)exp(-2£»r)]. 

(36) 

These expressions describe a narrow resonance in 
the spectrum of the beat signal. The maximum of 
this resonance coincides with the frequency 



o>b = y c T and the width of this resonance is equal 
to 



r res = 2/(x T n w ). 



(37) 



Neither of these parameters depends on the fluctu- 
ation properties. The only dependence on the fluc- 
tuations appears in the amplitude of the signal, 
which depends on the intensity of the noise as 
parameterized by D. For small DT<g: 1 the magni- 
tude / rcs is proportional to D J 2 . There is an optimal 
value of DT= 0.89664 for which the amplitude is 
maximum, with value / rcs = 0.596977t[/i„A 2 ] 2 . 



5. Fluctuations of the optical length 

Whereas uncontrolled fluctuations of the seed- 
laser frequency can be useful, any fluctuations of 
the optical path will degrade the accuracy. In this 
section we estimate the effect of such fluctuations. 

Consider the influence on the resonance profile 
(31) caused by fluctuations in the optical length of 
the cavity. We write the optical length as 
L(t) - L 0 + <5£(0, where SL(t) is a zero-mean fluc- 
tuation of the optical length. Taking into account 
these fluctuations it is easy to show that Eqs. (13), 
(14) for the field can be rewritten as 

oo 



n=Q 



- i(co s + nA)t)exp[i5&„(t)], 



(38) 



where 



14 I C 



8L(t - /T r ) 



--^X>('-/Tr). 



n-l 



(39) 



/=0 



For the conditions of the observations of the 
narrow resonance (31) the correlation function 
(21) can be written as 



oo oo 



G(r) = M 4 Yl £ l a »| 2 H 2 ex P - ; ) Ar ] 
n=0 /=o 

x F_(r - h t + m r , T)F 3 (T, x, l,n), (40) 
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where 

F 3 (7, t, /, n) = (exp \-\5<P n {t) + \S<P n (t + r) 

- i<5tf>,(r + t + r) + i50t(t + t)]). 

(41) 

To evaluate the functions F 3 (r, t, /, «) we assume 
that the fluctuation 5L(t) of the optical length is 
an Ornstein-Uhlenbeck process [39,40] whose 
parameters we denote by Gl and D L > 

(3L) = 0, 

< tt (W)> = GiD £ exp(-G L |* - *U (42) 

We use formula (A. 5) to evaluate stochastic aver- 
ages of exponentiated noise and obtain for 
F 3 (T, x, /, ri) the result 

F 3 (r, t, /, n) = exp {-[*„„ (0) + g„(0) 

+ Wlh (43) 

where 

n I 

SnM =J2Y1 T ' G ^ eX P [-G L |f + (iff - t)t r |]. 

(44) 

We assume that fluctuations of the optical length 
are relatively slow, G L x T «C 1. This is always true 
for fluctuations caused by variation of pump inten- 
sity or vibrations. Then we replace the summation 
over m and k in Eq. (44) by an integration, and for 
the most interesting case nG L x vy lG L x x » 1 write 

g nl (t) ~ G L D L / df, / df 2 exp(-G £ |f + t t - h\) 
Jo Jo 

^2D^-^[exp(-G,M) 

+ exp(-G L |f + nx r - lx T \)l (45) 

where t< is the lesser of t\ = nx r and t 2 - lx x —t. 
Using Eq. (45) for small T <gi \IG L one can obtain 
the expression 

F 3 (7\t, /,*) = exp {-v[2 - cxp(-<? L |T|) 

- exp(-G L |r + nx T - /r r |)]}, (46) 

where 



v = T 2 D t G L . (47) 

Note that we made no assumptions about the 
intensity D L of the fluctuations. Thus the quantity 
v can be large, and we cannot generally expand the 
exponential. For the spectrum J&jS) we obtain the 
formula 

M*>) = £ WW f°dTF„(T,n 

x exp (— v[2 — exp(— G £ |t 

+ (/-«)T r |)-exp(-G £ |T|)]) 

x exp [iarr + ico(/T r — nx T ) + i(« — /)AT]. 

(48) 

As follows from the above consideration of the 
delta-correlated fluctuations of the seeding laser 
the function F_(x,T) can be considered as a very 
narrow function of t. Therefore we use the approx- 
imation (16) for the amplitudes a„ and write the 
spectrum Js(<*>) as 

1 f°° 

J s (co) = Jres-7= / d*exp (-x 2 ) 

V n J-oo 

x exp{-v[l exp(-G L T r /7 w |x|)]} 

x exp [-\xn w (co - y c T)x T ] (49) 

with y res given by (36). 

This is a general result valid for any choice of 
the noise parameters. For fluctuations originating 
in various mechanical or technical variations, G L 
is about 10 4 to 10 5 s*" 1 . Then for r r = 1 ns we have 
n w G L T T = nMO~ 4 -lO~ 5 )- Thus there are two dif- 
ferent possibilities. The first one is a laser with 
moderate numbers of the frequency components, 
n w «C 10 4 , and the second one is a laser with broad- 
band gain (Ti-Sapphire laser) with n w ^ 10 4 to 10 5 . 
For given laser properties the two possibilities cor- 
respond to correlation times \IG L much longer 
than or much shorter than n w x r . We consider these 
two cases separately below. 

5. 1. Long correlation time 

For the first case, when G L nw x T 1, one can 
expand exp(— G L x r nw \x\) as 

exp(-G/:T r rt w |*|) = 1 - G L x T n w \x\ (50) 
and obtain the resonance shape as a Voigt profile: 
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1 f°° 

J s {oj) = J Tts —7= / dx(-x 2 - vG L x T n w \x\) 

V K J-oo 

x exp [—ixn w (a> — yT)x T ] 



sftn w x t vG L + i(a> -yT-x) 

(51) 

We now consider two limiting cases of the expo- 
nential argument vG L Xrn w . 

• For vG L Xrii w 1 the fluctuations of the optical 
length do not influence the resonance shape. As 
noted earlier, in this case a narrow resonance 
exists in the spectrum of the beating signal. 
The maximum of this resonance coincides with 
the frequency a> B = y c T y and the width is 
r TCS = 2l(x,n w ). 

• For vGzXrn w » 1 the spectrum J^cS) is 
Lorentzian, 



J s (co) = J n 



1 



\ + (a>-y c T) 2 /(vG L ) 2 ' 



(52) 



The maximum of this resonance again coincides 
with the frequency co B = Vc^but the width r rcs 
of this resonances is larger than 2/(x r nw ), and is 
equal to r res = vG L = T 2 D L G 2 L . The amplitude 
Jres of the resonance is equal to J ns — J^lj 
(y/nvG L x T n w ) and is much smaller than / rcs . 

These formulas show that even with slow fluctu- 
ations there will be a limit to the accuracy, i.e., a 
finite bandwidth of the beat spectrum. 

5.2. Short correlation time 

Consider now the second case, when 
G L n w x T ^>\. For large detunings \oj— yl\^\l(n 
w x T ) the main contribution in the integral over x 
in Eq. (49) and so this integral can be written as 

exp(-v) 



J s (oj) =J l 



-2Re 



y/nG L x T n„ 

J poo 
I dy{exp[+vexp(-.>')] - 1} 
o 

(<o-y c T)] 



x exp — ly 



(53) 



This integral describes a resonance which is close 
to Lorentzian with width r rcs = G L . The ampli- 
tude of this resonance is equal to 



Js(0)=J n 



2exp(-v) 



y/iiG L x r n w 

Jroo 
I dy{exp(+vexp(-j))-l}. (54) 
o 

This expression has the limiting values 
2 f v for v 1, 



Js(0)=J» 



y/nG L x r n 



v for v <3C 1 . 



(55) 



The maximum value of this function of v is 
1 3\2J TC J(GiXrnJ), occurring for v = 1.503. 

For small detunings \co — y 7] ^ \ln w x T the main 
contribution to the integral over x in Eq. (49) 
comes from large |x| ^ 1 . The result can be written 
as 

J s (co) ~ exp(-v)J res 

x exp [-(a> - y c T) 2 (x r n w /2) 2 ] + J S (0). 

(56) 

Thus for relatively fast fluctuations, meaning 
HwGL*r ^> 1 , and when there is only partial correla- 
tion between changes of the phases of the fre- 
quency components of the laser spectrum, the 
beat-frequency spectrum Js(co) is the superposition 
of a narrow resonance (56) and a wide resonance 
(53). The ratio of the amplitudes of the narrow 
and wide resonances is y/nx T n w /(2T 2 D L ) » 1 for 
small v = T^DlGl 1 and is exponentially small, 
Qxp{-v)y/nG L x x n w v/2 for large v»l. (Fig. 4) 
shows the beat signal for various v and G L . 



6. FSF interferometry with an FM seed 

Earlier discussion emphasized that some modu- 
lation of the seed laser is essential if an interferom- 
eter beat signal is to be observed. The needed 
modulation will automatically occur if the seed la- 
ser is noisy. But one may also consider introducing 
structured modulation, i.e., a frequency-modu- 
lated seed laser. Here we analyze this possibility. 
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Let us consider the case when the seeding laser 
is frequency modulated. We take the phase of <p s {t) 
to be, instead of a stochastic quantity, a sinusoid 

ft(0 = «* + 0Mn(O-O. (57) 

where Q m is the modulation frequency, P is the 
modulation index and q> 0 is the initial phase. Then 
the FSF laser field (13), (14) reads 



E(t) = ^E m (t)aty[-\{to M + nA)t] % 



n=0 



where the amplitude is now 

£«w = £ c Yl a » ex p h*« - 

n 

-tfsin [Q m (t - /it,)]]. 



1.0 

0.5 

0.0 
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Fig. 4. The evolution of the beat-frequency spectrum Jsico) in 
the case of short correlation time (G L Xrn w = 20) with increasing 
of the optical length fluctuations intensity: (a) no fluctuations 
v = 0; (b) v = 1; (c) v = 5 and (d) v = 10. 



The homodyne beat signal S(t) given by Eq. (18) 
can be rewritten as 

5(0 = exp(iA*)|et| 2 £ a n a] 

nJ=Q 

x exp {-i(#„ - #/) - ijSsin [Q m (t - nx T )\ 
+ ifi sin [Q m {t + T - /x r )]} exp \-\{a) s + /i J)/ 
+ i(o>, + /zl)(f + r)]+c.c. (60) 

We are interested in the part of the signal S(i) that 
oscillates with the seed modulation frequency Q m . 
More particularly, we consider slow frequency 
modulation Q m — y c T <§c A. To evaluate such mod- 
ulation we need only the terms in Eq. (60) for 
which l-n. Then we have the expression 

oo 

5(0 = exp(izltf>)|£ c | 2 Yl l a »! 2 ex P H-^* + nA ^ 

x exp { +2i0 sin [Q m T/2] 

x cos [Q m (t - nx T + T/2)]} + c.c. 

Using the formula 

exp(iMcos0)= iV„(M)exp(i/i4>), (61) 

n—-oo 

where J n (x) is the nth-order Bessel function, we ob- 
tain the result 

5(O=exp(iA0)|e c | 2 ^K| 2 

xexp [\lQ m {t - nx T + T /2)] 
x exp [+i(co s + nA)T] 

x f2 i'M2pu* [GmT/2]) + c.c. (62) 

The processing of this signal can proceed in two 
ways: either by means of a phase sensitive detector, 
or else with a detector that is insensitive to phase. 
As we will show next, there are notable advantages 
for the latter. 



6.1. Phase sensitive detectors 

Consider phase-sensitive detection of the modu- 
lation signal. Typically this produces, as a function 
of modulation frequency Q m and synchronous 
detector phase & d9 a signal 
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(63) 

With the same assumption which were made to ob- 
tain Eq. (35) the signal of Eq. (62) can be written 
as 

= -«H-4 2 |£ c | 2 ^|/i(20sin [Q m T/2]) 

x cos [(-«„ + yj)x t n m + Q m T/2 
-<P d + (o s T + A<f>] 



x exp 



8 



• (64) 



An examination of the argument of the exponen- 
tial reveals that the signal Sphasc is a narrow reso- 
nance function of the detuning of the modulation 
frequency Q m from the frequency y Q T. This signal 
depends, through the argument of the cosine, on 
the phase 

(-ftn + 7 c T)r T n m + Q m T/2 - <P d + co s T + A<f>. 

(65) 

This dependence can lead to a resonance shape 
that differs significantly from a Gaussian form. 

6.2. Phase insensitive detection 

Phase sensitive detection of the interference 
beat signal introduces a dependence upon the 
material phase shift A<f>. Thus it will be much more 
convenient to use amplitude detection. Then the 
signal is proportional to the square of the ampli- 
tude of the variable components, 



Sampl ~ \/ (S<p d=0 ) 2 + (S<p d=n /2) 2 

= J, (20 sin [Q m T/2))\e c \ 2 n^ 2 yj^ 
(ft, - yT) 2 x 2 n 

x exp 



8 



(66) 



Here there is no dependence upon the material 
phase shift A$. Note that this signal can easily 
be obtained using a standard RF-spectrum ana- 
lyzer. Regarded as a function of modulation fre- 



quency the signal appears as a Gaussian profile. 
The width of this intensity profile is the same as 
that of the amplitude profile shown in the beat-fre- 
quency spectrum without cavity fluctuations of 
Fig. 4(a). 

To make use of the formula given in Eq. (66) 
for a determination of A, one would vary the FM 
frequency Q m to find the value £2 max for which 
Sampi is a maximum. The delay time is then 



T = Q m /y c 



(67) 



Once an approximate delay time is established, the 
strength of the peak value can be optimized by 
adjusting the modulation index . For application 
to surface profilometry, the reference arm of the 
interferometer would be adjusted so that the beat 
frequencies have values within a convenient range, 
much larger than the width r rcs , and the modula- 
tion index /? would be adjusted to optimize heights 
within a range of values. 



7. Achievable accuracy 

We will show here that the accuracy of interfer- 
ometry based on the FSF laser can be comparable 
to the accuracy of conventional interferorhetry, 
i.e., a fraction of a wavelength. The analysis starts 
by expressing the error in step-height measurement 
Sh in terms of beat frequency error 5a> as 



dh = ^-3o) = ^-3oj. 



(68) 



We can make this distance error small, for fixed er- 
ror in beat frequency, by increasing the chirp rate 
y c , by making the cavity circumference L smaller, 
or by increasing the frequency shift per round trip 
A. Ultimately, however, the achievable error is set 
by the accuracy with which one can measure the 
beat frequency. The accuracy of this measurement 
is determined primarily by the width of the beat- 
frequency spectrum r rcs . We take the error to be 
some fraction e of this width: 



So) = ef n 



(69) 



and write the height error as proportional to a 
coherence length h c , 
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dh = eh c . (70) 

We can then use the result r rcs = 2ln w x T and y c = Al 
t r to write this parameter variously as 

cr^s = _c__ _ c 

n w A Sco L 1 



h c = 



2y c 



(71) 



where Sa> L = n w A is the width of the FSF spec- 
trum. In reference [26] we give several estimates 
of this width. We show, for example, that 3co L 
and n w are related to T f , the bandwidth of the fil- 
tering mechanisms within the FSF cavity, through 
the approximation 

8a> L = n w A « (Ar 2 f ) l/ \ (72) 
Thus we can write 



(AT?) 



•2x1/3* 



(73) 



For example, with filter bandwidth of 100 nm for 
light with central wavelength 800 nm (a bandwidth 
of 47 THz), and a frequency shift per round trip of 
200 MHz, this length is 63 \im. 

To estimate e we here consider the error present 
in various realizations of FSF interferometry. 

7.1, Phase noisy seed 

First we consider the FSF laser to be seeded by 
a cw laser subject to phase noise. When the reso- 
nance width r rcs is sufficiently small, one can con- 
sider a direct measurement of the beat frequency 
using a frequency meter to count cycles within 
some observation time interval T 0 . This time 
should be significantly longer than the resolving 
time for the resonance width, T 0 ^> l/r rcs . If we 
consider a measurement time of T 0 = 10~ 4 s, then 
it is necessary that the resonance width be larger 
than 10 kHz, meaning a width of a few tens of 
kHz, in order that T 0 » l/r res . In this case c is 
roughly 

6^ . 1 . (74) 

V^resT'o 

This same accuracy is obtainable with other meth- 
ods. For example in the frequency domain, using a 
spectrum analyzer that scans the frequency, one 
can measure the frequency of the peak in the spec- 
tral distribution. The result is again Eq. (74). 



7.2. Frequency modulated seed 

Next we consider using a frequency modulated 
seed and measuring the modulation frequency 
Q m for which the output signal of Eq. (66) maxi- 
mizes. There are several possible contributions to 
the error in measuring the frequency of this peak 
value. 

To use the resonance (66) for the measurement 
of the delay time T one can additionally slowly 
modulate the frequency Q m . We assume that 
G m {t) = + 8Q m cos{2nft) where / is a slow 
modulation frequency and neglect optical length 
fluctuations. Then for small modulation amplitude 
SCl m and detunings (dQ my \Q ( ^ -y c T\ «: r res ) the 
magnitude of the detected useful signal varies line- 
arly with modulation detuning Q 1 ^ — y c T, 

Useful - J\ (20 sin [Q m T/2])\e c \ 2 n w A 2 ^ 
^ (O m - yJ){x T n w f 



(75) 



The noise signal S noisc can be estimated using (34) 
for delta-correlated phase fluctuations and small 
DT<: 1 as 



Snoise = AH\/^ 

~ V^\e c \ 2 n^ 2 y/FDy/S^, 



(76) 



where 8co d = 1/ T 0 is the filter passband of the 
detection system. We take the frequency error Sa> 
to be the value of the detuning Q m - y c T where sig- 
nal and noise are equal, S noisc = S uscfu i. This leads 
to the estimate 



Sco ^ 



2(r res ) 2 v / r^y^ 



y, (20 sin [Q m T/2))5Q m 



(77) 



For rough estimation one can take 8Q m ^ r res 
and, with appropriate choice of the modulation in- 
dex 0, y,(20sin[Q M r/2]) ~ 0.5. Then 



(78) 



8a>~4r^T 2 D5(D d , 
which leads to the result 

€ ~ 4Ty/D5^ d = J— x ATy/WZ, (79) 

V 1 0* res 
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Table 1 

Estimates of height errors showing values for three methods of measurement: (1) Phase noisy seed , (2) FM seed, with spectral analysis, 
(3) quantum noise limit 



measurement time [ms] 


0.1 


1 


10 


100 


T f = 5 GHz, h c = 3 cm , n w = 7 










h error [nm] noisy seed 


5x 10 6 


1.6 x 10 6 


5x 10 s 


1.6 x 10 5 


h error [nm] FM seed 


8x 10 4 


2x 10 4 


8x 10 3 


2x 10 3 


h eiTor [nm] quant, lim. 


80 


25 


8 


2 


r t - 47 GHz, h c = 0.8 cm, n xv = 31 










h error [nm] noisy seed 


2.5 x 10 6 


8x 10 s 


2.5 x 10 s 


8x 10 4 


h error [nm] FM seed 


1.6 x 10 4 


5x 10 3 


1.6x 10 3 


5x 10 3 


h error [nm] quant, lim. 


17 


5 


1.7 


0.5 


T f = 47 THz, h c = 77 nm, n w = 3 x 10 3 










h error [nm] noisy seed 


2.5 x 10 5 


8x 10 4 


2.5 x 10 4 


8x 10 3 


h error [nm] FM seed 


1.6x 10 2 


51 


16 


5 


h error [nm] quant, lim. 


0.2 


0.05 


0.02 


0.005 



where D is the bandwidth of the seed laser. For 
example, with D = 2n x 1 MHz and h = 1 cm, the 
e of Eq. (79) is 2 x 10~ 4 times smaller than Eq. (74). 

7.3. Quantum noise 

Spontaneous emission in the seed laser itself sets 
a lower bound on the achievable bandwidth D. As 
discussed in textbooks on laser theory, this Schaw- 
low-Townes limit can be as small as a few Hz. This 
bandwidth imposes a fundamental limit on the 
accuracy of the FM method: when all "technical" 
noise has been eliminated, there remains an irre- 
ducible bandwidth attributable to white noise. As 
an example, for a linewidth of D = 2n x 1 Hz, 
T= 0.5 ns and r = 2n x 10 kHz the factor 
4Ty/DT\Z is 10" 6 . 

Another ultimate limit on observation comes 
from shot noise in the detector. We estimate this 
to be 

e "\Z^(?) = 7f=^ x V^' (80) 

where P is the laser power and r\ is the detector effi- 
ciency. i* P hot is the photon counting rate of the 
detector (laser power times efficiency divided by 
photon energy fica). With reasonable estimates of 
these parameters (say P = 1 W, r\ = 0.5) this value 
of e is 2 x 10~ 7 times smaller than Eq. (74), and an 
order of magnitude smaller than the Schawlow- 
Townes limit mentioned above. 



7.4. Examples 

Table 1 gives examples of the accuracy possible, 

using the preceding formulas and the following 
parameters 

X - 800 nm A = 2n x 200 MHz 

h = 1 m D = 2n x 1 MHz 

P=1W r] = 0.5 



The three sections of the table illustrate results 
for several choices of bandwidth T f , height h and 
number of discrete components n w . As can be seen, 
modulation of the seed offers here a hundredfold 
improvement in accuracy compared with the reli- 
ance on phase noise. The accuracy is still far from 
the quantum limit. 



8. Summary and conclusions 

The FSF laser offers useful application for high- 
accuracy interferometry. It allows determination 
of distance by measuring frequencies, and is insen- 
sitive to material phase shifts that occur upon 
reflection. As we have shown, it is essential that 
the laser output have available a range of frequen- 
cies that include the beat frequency y c T, for any T 
of interest. Although the reliance on noise can pro- 
vide the needed frequencies, we suggest an alterna- 
tive based on frequency modulation. 
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The basic principle is that of chirped radar, here 
implemented with a non-trivial modification in the 
form of a FSF laser. The greater chirp rate and ex- 
treme linearity of chirp make this far superior to 
conventional methods. 

The position error Sh can be made small, for 
fixed error in beat frequency, by increasing the chirp 
rate y c , by making the cavity circumference L smal- 
ler, or by making the frequency shift A larger. 

Our theory provides analytic expressions for the 
intensity of the beat signal as a function of fre- 
quency. For any seed frequency, the beat signal 
has a spectral distribution sharply peaked about 
the value o>b = {2yjc)h and thus a measurement 
of the position of the spectral peak provides a 
measurement of h. The theory also provides a con- 
nection between the spectral width r rcs and laser 
parameters, specifically the effective number of dis- 
crete frequency components n w within the output 
spectrum of the FSF laser (i.e., the output band- 
width divided by the frequency shift A). It is this 
width that ultimately restricts the possible resolu- 
tion by limiting the accuracy with which the beat 
frequency can be determined. Under appropriate 
conditions the error in beat frequency determina- 
tion can provide subwavelength accuracy compa- 
rable to, or exceeding, what is obtained with 
conventional interferometry. 

To achieve the accuracy 5h = 100 nm with a fre- 
quency shift of A = 2n x 100 MHz one needs 
n w « 10 5 . This value is very large, but may be 
achievable with a Ti-Sapphire laser or fiber laser. 
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Appendix A. Stochastic averaging 

A. 1. The noise model 

We consider in this appendix the evaluation and 
behavior of the functions F±(co i T) 9 defined as Fou- 
rier transforms 



F ± (a>,T)= / dTF ± (r,r)exp[ia;T] (A.l) 

J -OO 

of the functions 

F±(tuh) = (exp[-i?,(0 + i9s(' + fi) 

Ti<p s (t + t 2 )±i<P s {t + t 2 + t l )]), (A.2) 

where (• •) denotes an averaging over realizations 
of the stochastic phase <p s (t). These functions 
incorporate all of the effects of noise upon the 
interferometer beat spectrum. 

We assume that the time derivative of the phase 
<p s (/) is a stochastic process £(/)> 

m = «0. (A.3) 

with zero mean and exponential correlation [£(t) is 
an Ornstein-Uhlenbeck process] 

<£) = 0, mm) = GDexp(-G\t - t'\). 

(A.4) 

We then evaluate stochastic averages of exponenti- 
ated noise, using the formula 

<exp(i#)> = exp[-i<* 2 )]. ' (A.5) 

For F_(7*,t) one obtains the result 

F_ (r, , t 2 ) = exp (- 1 [g(t, t )+g(t + t u t + ti) 

+ g(t + tt +t 2 ,t + h+t 2 ) 

+ g(t + t 2 ,t + t 2 )-2g(t,t + ti) 

+ 2g(t, t + t 2 + h)- 2g(t, t + t 2 ) 

-2g(t + t u t + t 2 + h)+2g(t + ti,t + t 2 ) 

-2g(t + t 2 + tut + t 2 )]j, (A.6) 

where, with f< as the lesser of t\ and t 2 , 
g(t u t 2 ) = (<p s (ti)<P s (h)) 

= r d/' r d^«(o«(o> 

J— OO J — OO 

= GD f" dr" P dr"exp(-C|f" - <\) 
Jo Jo 

= 2Dt < - -[1 - exp(-Gr,) - exp(-G/ 2 ) 
+ exp(-G|f 2 -f,|). (A.7) 
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The functions F±(tiJ 2 ) are even function of t\ and 
/ 2 . For positive t x and t 2 one obtains the results 

F_(f,,r 2 ) = exp -2£>r< - ^{-2 + 2exp(-Gr,) 
+ 2exp(-G/ 2 ) - exp(-G|/ 2 - 
-exp(-G[* 2 + *,])}], (A.8) 



-2D(2/, - /<) - -{-2 + 2exp(-G/,) 



F + (/i s / 2 ) = exp 

- 2exp(-Gf 2 ) + exp(-G|f 2 - 
+ exp(-G[/ 2 + /,])}]. 



(A.9) 



These expressions, when used with formulas for 
Js(o)), allow one to calculate the spectrum of beat- 
ing (24) for arbitrary values of the phase fluctua- 
tion parameters G and D, the delay time T and 
the laser parameters. 

A.2. Slow fluctuation 

The preceding subsection dealt with an Orn- 
stein-Uhlenbeck process governed by parameters 
G and D. To get a simpler estimation we consider 
further two limiting cases: (1) the phase changes 
slowly, G D and (2) the noise is delta correlated, 
G — > oo. Bearing in mind the optimal conditions 
for narrow resonance (24) observation we will deal 
only with the function F_(t,7). In the first case, of 
slow phase fluctuations, the expression for F_ 
becomes 



F_(T,r) 



exp [-l/3ZX7V(3r - t)], t ^ T, 
exp [-l/3DG 2 T 2 (3t — T)], x > T. 

(A. 10) 

For weak fluctuations, DG 2 J 3 <1, the character- 
istic width of F_(t,T) as function of t is <5t = 1/ 
(DC 2 !*). Hence in this case the Fourier compo- 
nent F_(a>,7) is a narrow function of co with char- 
acteristic width Sco = DG 2 T 3 ± < f . As the 
fluctuation intensity increases the function F_(t,7) 
becomes narrower and hence the width of F-(a>,T) 
increases. For DG 2 T* » 1 we obtain a Gaussian 
profile for F_ 



F.(co,T) = T _^l_ exp . (A.ll) 

The width of this is Sco = 2VDG 2 T 3 ± » I . 

A. 3. Rapid fluctuations 

Here we assume the noise fluctuates rapidly, 
and consider the limit G — » oo when the correla- 
tion time is infinitesimal, i.e., the noise is delta cor- 
related. Then the expression for F_ becomes 

F-(t u t 2 ) =«p[-2D/<]. (A.12) 

When co ^ 0 we obtain for the spectral function 
F_(<x> ,7) the formula 

F. { ^T)=2R^ l -^^ 2DT + la}T) 



-exp(-2£>r) 



2Z> - ico 
1 - exp (/a>r) 



-ico 



(A.13) 



For small DT <zi 1 this function has the width 
<5a>~J;. With increasing DT^> 1 the width of 
.F_(aj,7) increases. The function is Lorentzian 



F_(cu,r) = — , 



with width 5co = £>/2. 



(A. 14) 
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Abstract 

This paper demonstrates that a frequency-shifted feedback laser, when seeded by a phase-modulated narrow-band radiation field, is a 
powerful tool for distance measurements to accuracy better than 10 um and resolution better than 100 jam, for distances of a few meters. 
In such measurements the unknown distance forms one arm of a Michelson interferometer, in which the intensity of the output signal is 
modulated at the phase-modulation frequency of the seed. The amplitude of the output-signal modulation exhibits a resonance for every 
distinct signal delay, i.e. for each distinct distance within the laser spot on the target. The use of a phase-modulated input seed allows one 
to use a very narrow-bandwidth filter when measuring the return signal. The results reported in this paper are in excellent agreement with 
previous theoretical predictions [L. Yatsenko et al., Opt. Commun. 242 (2004) 581] for the resolution limit and high signal-to-noise ratio 
for this new technique. 
© 2006 Elsevier B.V. All rights reserved. 

PACS: 42.55.-f; 42.60.Da; 42.55.Ah 
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1. Introduction 

Laser based techniques - variants of radio-frequency 
RADAR - find increasing use for a variety of ranging 
applications, including measurements of distances with 
accuracy approaching or less than optical wavelengths. 
Such distances may be meters or kilometers, as is the case 
with LIDAR, or millimeters or centimeters, as is the case 
with typical workbench interferometers. 

We here describe, and demonstrate, a new technique for 
accurate distance measurements using a frequency-shifted 
feedback (FSF) laser seeded by a phase-modulated (PM) 
laser. As we will note, FSF lasers differ significantly from 
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ordinary lasers, and PM seeding has particular advantages 
when the laser is used for ranging. 

In all FSF lasers, a feedback mechanism introduces a 
fixed frequency shift A A om with each passage of a wave 
packet around the laser cavity (during the round-trip time 
T r ). After its first demonstration by Kowalski et al. [I], 
FSF lasers were further studied and developed by Berg- 
mann and co-workers [2]. Later Nakamura et al. further 
studied such systems and demonstrated the use of FSF 
lasers for ranging [3,4]. They pointed out that the output 
of a frequency-shifted feedback (FSF) laser can be 
regarded as a succession of chirped pulses that can be used 
to measure distances, in a manner akin to chirped 
RADAR. In such operation one measures a frequency off- 
set proportional to the time delay t = 2h/c between a 
probe and a reference pulse, and from this one evaluates 
the incremental distance h. 
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In a previous article we provided a theoretical descrip- 
tion of ranging measurements using a PM-seeded FSF laser 
[5]. The present article provides an experimental demon- 
stration of that work. We show, in particular, that the 
operation of the FSF laser with PM-seeded radiation offers 
significant advantages over previous approaches to ranging 
with a FSF laser. The PM-seeding technique encodes the 
signal into a narrow bandwidth. As we demonstrate, this 
enhances the signal-to-noise ratio while maintaining the 
high resolution associated with the large bandwidth of 
the FSF radiation. 

2. Laser ranging techniques 

Laser-based techniques are widely used for distance 
measurement; for a detailed critical review of the usual 
techniques for distance measurement see [6], In this section, 
we briefly discuss the more popular methods. These fall 
into three main categories: triangulation, time-of-flight 
and interferometry. 

An optical triangulation system [7] uses a laser as one 
corner of a triangle and a camera (or other recording 
device) as the second one. The distance to the object can 
be determined using simple trigonometry, if one knows 
the distance between the laser and the camera and the tri- 
angle angles. The technique is applicable for distances of 
a few millimeters to tens of meters. The accuracy decreases 
with range, and for commercial devices can be about 0.01% 
of the range. 

A time-of-flight (TOF) system [8] measures the time t 
between a pulse emission and the return of the pulse 
reflected from the object. The distance h is half of the prod- 
uct of this time and the velocity of light, h = tc/2. The pre- 
cision of such range measurements depends mainly on 
noise-generated timing jitter, walk (timing error created 
by variations in the time-pickoff circuit), nonlinearity and 
drift. By carefully designing the system to minimize system- 
atic errors, and by averaging successive measurements, one 
can achieve a precision of a few millimeters, although the 
accuracy of commercial devices is a few centimeters. A 
TOF system has the advantage of a weak dependence of 
the accuracy on distance. Furthermore, the technique does 
not need highly-reflective "cooperative" targets, such as 
mirrors. 

Conventional interferometers (e.g. [9]) use fixed-fre- 
quency continuous-wave light. Distances are evaluated by 
measuring (via superposition of the test beam with a refer- 
ence beam) the phase <f> = 2coh/c induced by propagation of 
the test beam through an incremental distance h to a 
reflecting (mirror) surface. Interferometric measurements 
provide only values of relative distances: they give the dis- 
tance that the mirror has moved from its initial position, 
rather than an absolute positional measurement. Further- 
more, distances obtained using a single-wavelength inter- 
ferometer have an ambiguity of half a wavelength. Most 
interferometric systems use mirrors as reflectors and fre- 
quency-stabilized lasers. The accuracy of the interferomet- 



ric measurement can be very high (much better than 1 nm), 
but to reach this high precision it is essential to minimize 
changes in the index of refraction of the air that occur as 
a result of changes in air pressure and temperature. 

3. Chirped-pulse interferometry 

For measurements of absolute distances without the 
interference phase ambiguity, chirped-pulse interferometry 
(or the technique of frequency-modulated continuous-wave 
laser radar) is often used [6]. The operation of chirped- 
pulse interferometry can be understood by considering a 
laser whose frequency co varies linearly with time (a linearly 
chirped frequency) according to the formula: 

where y c is the chirp rate. The chirped radiation passes into 
a Michelson interferometer having one fixed (reference) 
arm and one variable-length (measurement) arm, as indi- 
cated schematically in Fig. 1. A beamsplitter (BS) provides 
the beams for each arm. After reflection from mirrors, 
these beams are brought into coincidence on a photodetec- 
tor. When the measurement arm differs in length from the 
reference arm by h the beam in the latter arm incurs a delay 
T=2h/c. As a result, its frequency is shifted by 
Aco = y c T = Ihyjc relative to the frequency of the reference 
beam. The two overlapping beams produce an amplitude- 
modulated field on the detector, with the beat frequency 
co B — Aco. A measurement of this beat frequency provides 
a measurement of the distance, evaluated from the 
formula: 

h = cco B /2y c . (2) 

Single-frequency laser diodes are usually used as the light 
sources in chirped-pulse interferometry because one can 
easily modulate their frequency by modulating the diode 
current. The maximum measurable distance of a fre- 
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Fig. 1. Schematic diagram of a Michelson interferometer. Chirped-pulse 
laser light passes through a beam splitter (BS), and then along two arms, 
one of which contains the unknown height to be measured. Mirrors (or a 
reflecting surface) and the beam splitter bring these beams into coincidence 
on a photodiode detector. The height h is found relative to the reference 
arm length, as shown by dotted partial circumference. 
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quency-modulated CW laser diode system is ultimately lim- 
ited by the coherence length of the laser diode, which typ- 
ically is of the order of several tens of meters. 

In chirped-pulse interferometry the measurement accu- 
racy is inversely proportional to the tuning range of the 
laser frequency. As an example, using continuously tunable 
laser diodes with a tuning range of 500 GHz, a spectral 
width of 25 MHz and a repetition time 10 us (a chirp rate 
y c = 5 x 10 16 ), one can expect accuracy of 2.7 urn for 1 cm 
distance and 43 urns for 1 m [10]. The large modulation 
bandwidth of laser diodes makes them suitable for accurate 
and fast measurement of short distances with this 
technique. 

Practical realization of the potential advantages of 
chirped-pulse interferometry is hindered by the following 
problems. To achieve high accuracy the frequency sweep 
must be linear. A sweep nonlinearity will cause the beat fre- 
quency to change during the sweep and will thereby 
decrease the measurement accuracy. This effect is particu- 
larly troublesome for measurements of large distances. 
Furthermore, if there is no prior knowledge of the range, 
the detection system for a single-mode sweep requires a 
large bandwidth. For large distances, one needs to measure 
very high beat frequencies, unless one is able to adjust the 
length of the reference arm. All these problems can be 
solved. For example, some elegant methods have been sug- 
gested to minimize the errors due to sweep nonlinearity 
[11]. Nevertheless, the most attractive feature of the fre- 
quency-modulated CW laser diode system, namely its sim- 
plicity, disappears. 

These general principles of chirped-pulse ranging apply 
also when the laser is a FSF laser. However, there are sig- 
nificant differences, as we shall note. 

4. The FSF laser as an interferometer source 

A FSF laser differs from an ordinary laser by the pres- 
ence of an element that shifts the optical frequency (by a 
discrete amount) on every round trip of a wavepacket 
within the cavity. In our case, the frequency-shifting ele- 
ment is an acousto-optic modulator (AOM), wherein 
acoustic waves create a travelling grating. The arrangement 
is such that the first-order diffracted light (altered in fre- 
quency) continues within the cavity, while the zero-order 
emerges as the output. This construction eliminates the 
usual cavity interference that, by repeated passages of a 
wavepacket, builds up discrete fixed-frequency modes. As 
a result, a FSF has no longitudinal mode structure. The 
output radiation from a FSF laser has been discussed 
extensively [2,12]. It can be interpreted as a succession of 
frequency chirped components, a "moving comb" of fre- 
quencies, each of which is linearly chirped in accord with 
Eq. (1) [13]. Experiments have demonstrated that, as one 
might anticipate from such radiation, it is possible to use 
this output to conduct distance measurements [3,4]. 

Fig. 2 illustrates this comb-like behavior. Within the 
gain profile this comb has a large number of "teeth", each 
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Fig. 2. Schematic diagram of the "moving comb" in the frequency-time 
space. The frequency distance between "teeth" - the free spectral range 
(FSR) of the FSF laser - is 2nh f . The time separation is 2ti/A AO m- Every 
"tooth" has a Gaussian envelope. 



with a time varying amplitude B n (i). The result is the elec- 
tric field 

+OQ 

E MC (0 = <*o Yl fl »(') ex P I-Wi (01. (3) 

n =—oo 

where p n (t) is the phase of the nth component ("tooth"). 
The frequency of this component, a> ni undergoes a linear 
chirp 

O)„(0 = lp.(/) = aw + ^ (/ - «t s ), (4) 

where A AO m is the acousto-optic modulator frequency, r r is 
the round-trip time, t s = 2k/Aaom is the time separation 
between teeth and a> max is the frequency of the, FSF laser 
spectrum maximum. The chirp rate is y c = A AO m/ti- 

As in any laser, light in a FSF laser passes through a 
gain medium, where amplification occurs. Any spontane- 
ous emission from the excitation responsible for the gain 
will, through multiple passes (at successively shifted fre- 
quencies), contribute to the circulating field and hence to 
the output. For better control of the output one can inject 
a predetermined field into the cavity using a seed laser. This 
has been our approach. 

According to Eq. (4), in the moving-comb model the 
output from an FSF laser is a set of frequencies (comb 
teeth), separated by the cavity axial mode interval 
A c = 27i/r r , which are moving with the rate y c 

<o„(t) = ctw + y c t - "&c- (5) 

These frequencies are uniformly chirped, at rate y c , so in an 
interferometer a time delay T produces a number of beat- 
ing signals between different teeth of the reference and re- 
flected fields. In the absence of fluctuations the phases of 
all these beat signals are such that they sum to zero. How- 
ever, fluctuations break this balance, and one can therefore 
observe a resonant increase in the spectral density of the 
RF spectrum near these beat frequencies. Although such 
signals are weak and noisy, nevertheless, they have been 
successfully used for distance measurements of up to 
18 km [4]. 
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With the phase-modulation technique described in this 
paper, the maximum amplitude modulation of the FSF 
laser output takes place when the modulation frequency 
coincides with one of these beat frequencies. The simplest 
beats, those of lowest frequency, are those that occur 
between teeth having the same n. Such beats occur, at fre- 
quency ct) B = y c r, when the delay T is less than the time 
separation between adjacent teeth, t s = 2n/A AOMt so that 
the beat signal is less than the cavity frequency A c . How- 
ever, other beats occur, between tooth n and a tooth 
n + 5/i, with integer 5/i, at the beat frequency 



o) B (dn) = 6/i Ac + y c T. 



(6) 



The existence of this family of beat frequencies allows us to 
choose one that falls within a frequency range for which 
electronic signal processing is convenient. It is only neces- 
sary that we determine 5/i in order to deduce the delay T. 
This we can do by observing the effect of a slight alteration 
of the cavity frequency from A c to A c + S c , as produced by 
a change of the cavity optical length. This changes the ob- 
served beat frequency by 

5ct> B (5n) = 5/i x <5 C (7) 

(with the neglect of the small change in chirp rate). From 
knowledge of <5 C and a measurement of 5o> B (5/i) we can 
deduce Sn. 

Using this technique it is possible to measure large dis- 
tances without changing the reference arm or measuring 
very high-frequency beats. By making frequent measure- 
ments of 8/? it is possible to track a moving target. The very 
regular progression of chirped frequencies in the moving 
comb model provides an ideal linearity of sweep. 

The following sections describe the realization of this 
new technique for ranging with a PM-seeded FSF laser. 
The technique allows a dramatic enhancement of the sig- 
nal-to-noise ratio (up to three orders) while preserving all 
the attractive features of the FSF laser expected from the 
moving comb model. 

5. Seeded FSF laser 

The output of any FSF laser can be used for optical 
ranging, whether the laser field grows from spontaneous 
emission within the cavity or from a fixed-frequency exter- 
nal-seed laser. The present FSF laser differs from those of 
previous demonstrations by using a seed laser whose phase 
&(t) — co s t + cp s (t) is modulated sinusoidally, 

<Ps(*) = <p 0 + ps\n(Qt) (8) 



as was suggested in [5]. Here, co s 
Q is the modulation frequency, 
and (p 0 is the initial phase. 

Such a phase-modulated seed 
laser output, consisting of a set 
nents with the frequencies co s 
The phases of the components 
contribution that is quadratic 



is the seed laser frequency, 
fi is the modulation index 

generates a specific FSF- 
of discrete optical compo- 
+ "A A om (" = 0,1,2,...). 

have a time-independent 
in n and a time-delayed 



contribution from the phase-modulated seed laser phase 
<p s (/ — /ir r ). The amplitudes of the components are propor- 
tional to a Gaussian function of n y cxp[-(n - n m ) z /n^], 
where /i w is the effective number of optical components 
within the Gaussian profile of the laser spectrum and n m 
is the central position of the profile. These two character- 
izing parameters are large, /i w » 1 and n m » 1, and can be 
evaluated from measurable laser characteristics (for details 
see Eqs. (13)— (17) in [5]). When this output is used as 
input to a Michelson interferometer whose arms differ in 
length by A, there occurs in the output from the interfer- 
ometer an amplitude modulation with the frequency G. 
The amplitude S AM of this modulation is 

Sam =SoJi (2p sin [QT/2\) exp [-(Q - y c T) 2 T 2 n 2 w /a] , (9) 

where J\(x) is the Bessel function of order 1 and So is a nor- 
malization factor that incorporates the laser output power. 
This signal has a Gaussian profile that reaches a maximum 
when the modulation frequency G is equal to 

G M = y c T = y c 2h/c. (10) 

Thus, by measuring the modulation frequency for which 
the interferometer output maximizes we measure h. The 
constant of proportionality between the distance h and res- 
onance-peak frequency G M is c/2y c . For our setup this had 
the value = c/2y c = 0.02 mm/kHz. 

The accuracy with which one can determine the reso- 
nance maximum, and hence the accuracy 6/i with which 
one can measure A, is set by the width of this resonance, 



AG = 



2\/2 

T r rt w 



and can be written as 



6h = cAG 



2y c 



(11) 



(12) 



Here, the small parameter e 1 depends on the signal-to- 
noise ratio and, as shown in [5], its value can be smaller 
than 10~ 3 . It is worth noting that for such accuracy it is 
important to account for a shift of the resonance (9) due 
to the influence of the frequency dependent Bessel function. 
For a small modulation index p «: 1 this shift can be esti- 
mated as 



3Q = cot 



AQT 



AG. 



(13) 



Using values typical for our setup, T= 2.5 ns, G M = 2n • 
12.5 MHz and AG = 5 kHz we obtain the estimate &G ^ 
0.8 x 10" 3 A&. 

In the present work we scanned the frequency of phase 
modulation and measured the signal 5 A m as filtered by a 
narrow bandpass RF filter set to match the phase modula- 
tion frequency G. From the envelope of filtered signals we 
evaluated the frequency of maximum signal, and from this, 
using Eq. (10), we obtained the distance increment h. 
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6. Experimental setup 

To demonstrate this technique of PM-seeded FSF for 
optical ranging we used the setup sketched in Fig. 3. The 
FSF laser portion consists of an optical fiber incorporating 
a gain section, closed into a ring configuration by means of 
mirrors and the first-order diffraction output of an AOM 
modulator. This is seeded by a PM laser, and its output 
feeds a Michelson interferometer. 

The FSF laser has been described in [15]. Briefly, its 
parameters are as follows: 

• A AO m/27t = 80.3 MHz, the frequency of the AOM, 

• r r = 10.7 ns, the round trip time, 

• yjln = A A oM/2rcT r = 7.5 x 10 15 Hz/s, the chirp rate. 

The spectrum of our FSF typically is approximately a 
Gaussian, with width (FWHM) around 300 GHz, but 
under appropriate conditions this width can be extended 
to 3000 GHz (its shape is then no longer Gaussian). The 
spectral width also can be decreased to 60 GHz by inserting 
a Fabry-Perot etalon (FPE) with FSR 1000 GHz into the 
cavity. The FSF laser can operate in the wavelength region 
1030-1075 nm. For the distance measurements shown in 
this work we used the wavelengths 1040 nm (with 300 GHz 
and 60 GHz bandwidths) and 1055 nm (for 1000-3000 GHz 
bandwidths). 

The seed laser is a single-frequency external cavity diode 
laser constructed on the Littman-Metcalf scheme, as 
described in [14]. The phase modulation of the seed laser 
is produced by an electro-optic modulator (EOM) driven 
by a FM-generator that sweeps the frequency linearly. This 
generator is synchronized with the radio-frequency spec- 
trum analyzer (RSA) that contains a RF-filter. This system 
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Fig. 3. Schematic diagram of the layout used for distance measurements 
using a PM-seeded FSF laser. A rectangle encloses the interferometer 
portion of the apparatus. A signal generator drives both an electro-optic 
modulator (EOM), which modulates the phase of a seed laser, and adjusts 
the frequency of an RF filter which precedes the beat signal detector. P v is 
the optical spectrum of the FSF laser, Sam is the amplitude modulation 
signal. 



allows us to detect envelopes of the amplitude modulation 
and thereby measure distances. 

7. Example of optical spectra 

When the seed laser power is low, spontaneous emission 
dominates the FSF-laser field. When the seed power is very 
high its influence dominates and distorts the shape of the 
spectrum and, in turn, the envelope of 5am This change 
in optical spectrum is evident in Fig. 4, which presents typ- 
ical spectra of the FSF laser used in the present work, with 
and without an external seed laser. The frequency and the 
power (near 100 uW) of the seed laser have been chosen to 
maximize the amplitude modulation signal Sam- The spec- 
trum without a seed is close to a Gaussian in shape, but the 
seeded spectrum has marked asymmetry induced by the 
seed laser. 
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Fig. 4. Output spectral power /% (arbitrary units) versus frequency offset 
Av 8 from the seed frequency, with external seed (thick line) and without 
seed (thin line). 
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Fig. 5. Amplitude modulation signal S AM versus modulation frequency 
Q/2k for signals obtained with external seed (thick line, single scan) and 
without seed (thin line, averaged over 50 scans). 
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The use of a external PM seed allows one to limit the 
detector signal to the PM frequency, thereby eliminating 
noise from other frequencies. This is possible because the 
spectral width of the amplitude modulation signal is very 
small (in our experiments it is less than the apparatus res- 
olution of 50 Hz). This property allows us to use a narrow 
RF-filter for detecting the desired signal. The benefit of this 
filtering is seen in Fig. 5, which compares the signal 5 AM of 
an unseeded FSF laser and one seeded by a PM signal. The 
signals have nearly the same width Aft, but whereas the 
first signal is weak and noisy (only after averaging 50 scans 
does it appear close to Gaussian), the second one has a lar- 
ger amplitude and does not require averaging, 

8. Bandwidth-resolution relationship 

In keeping with all Fourier-transform frequency-time 
pairings, one expects an inverse relationship between the 
bandwidth of the FSF laser and the resolution with which 
one can measure delay time and distances [13]. Fig. 6 illus- 
trates this relationship by presenting three examples of the 
spectrum of the seeded FSF laser, each with a different 
width, together with the corresponding signals of ampli- 
tude modulation. The upper plots are for the FSF laser 
with FPE 1000 GHz inside the cavity. The middle ones 
are for the usual FSF laser operating near the wavelength 
corresponding to the maximum of the Yb 3+ gain profile. 
The lower ones are for the FSF laser operating near the 
wavelength corresponding to the local minimum of the 



Yb gain profile, where joint action of AOM dispersion 
and gain frequency dependence results in nearly flat effec- 
tive gain and consequent larger spectral width. Theory pre- 
dicts [13] that the inverse of the width Aft of the signal 
profile Sam should be directly proportional to the band- 
width r of the FSF laser. We found this dependence, 
shown in Fig. 7, fits a straight line even if the spectrum is 
not Gaussian (as in the lower left hand frame of Fig. 6). 
This linearity demonstrates that the FSF laser maintains 




0 500 1000 1500 2000 2500 3000 



Av s [GHZ] 

Fig. 7. Dependence of the inverse width 2n/AQ of the amplitude 
modulation signal upon the FSF-laser spectral width. Circles are 
experimental points, straight line is a linear fit. 
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coherence of the acousto-optic modes across the entire 
spectrum of the FSF laser. 

9. Demonstration of spatial resolution 

The optical ranging technique relies on the observation 
of a return signal reflected from a target. In our basic 
experiments we use a smooth flat reflecting surface, ori- 
ented to reflect the light directly back into the interferom- 
eter. More generally, the target may have a distribution 
of surface heights within the laser-beam spot, producing 
multiple reflections; each height provides a different delay 
time and a different resonance peak in the signal S/^m- 

If 

these peaks are sufficiently well separated to be resolved 
they can be identified as originating at distinct heights. 
The relative magnitude of a signal provides a measure of 
the relative area, within the laser spot, at the corresponding 
height. 

To demonstrate the resolution available with our sys- 
tem we fabricated samples that provided two distinct 
heights. These samples were flat rectangular reflectors 
whose surfaces were interrupted by a single step. As 
the laser beam moves across such a surface it encounters 
first one height, then two heights, and then one height 
again, as illustrated in Fig. 8. In accordance with theory, 
the relative positions of the peaks remains fixed (these 
depend only on the fixed height of the step on the sam- 
ple), while the heights of the peaks vary in accordance 
with the relative portion of the beam spot that is occu- 
pied by the corresponding height. 

We used samples having two height steps: one of 1.0 mm 
and 0. 1 mm. The height difference of 1 mm is readily 
resolved even with a laser whose spectral width is 1 nm 
(near 300 GHz). Fig. 9 shows how the amplitude modula- 
tion signal 5 A m changes as the laser beam moves across 
the 1 mm step sample past the border between the two fac- 
ets. Relatively low peaks at the positions b, c and d in com- 
parison with a and e can be explained by the fact that 
sample facets are not absolutely parallel and so there are 
slightly different optimal adjustments for every facet. 

The step height of 0.1 mm on the second sample was 
resolved only by the FSF laser with spectral width of about 
7 nm (near 2000 GHz). We measured the step height to be 
1 lOum, with an accuracy of better than 10 um (see Fig. 10). 
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Fig. 9. Signals S AM versus modulation frequency Q/2n for the sample 
with step 5h = 1 mm for several positions of the laser beam, as shown in 
Fig. 8. 
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Fig. 10. Signal S AM versus modulation frequency £2/2n for the sample 
with step 8/i= 1 10 nm. The upper border shows the distances correspond- 
ing to frequencies, according to Eq. (10). 



10. Conclusions, outlook and summary 

We have demonstrated that a FSF laser, when seeded by 
a PM signal, is a powerful tool for distance measurements, 
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Fig. 8. Sample for the demonstration of spatial resolution: two reflecting facets differing in height by 5/i. Successive frames show positions used in making 
observations. 



V. V. Ogurtsov et al. I Optics Communications 266 (2006) 266-273 



273 



at least in range of a few meters, to an accuracy better than 
10 urn and resolution better than 100 um. 

In such measurements the unknown distance forms one 
arm of a Michelson interferometer, in which the output 
signal is modulated at the PM frequency of the seed; this 
output signal exhibits a resonance for every distinct signal 
delay, i.e. for each distinct distance within the laser spot on 
the target. The use of a phase-modulated input seed allows 
one to use a very narrow bandwidth filter when measuring 
the return signal The results reported here are in excellent 
agreement with previous theoretical predictions of a high 
signal-to-noise ratio and resolution limit for this new tech- 
nique [5]. 

Our current resolution is about 30 wavelengths, based 
on the criterion of resolving two resonances with minimal 
full width at half maximum equal 2.5 kHz (see Fig. 7). 
This is smaller than what is achievable with frequency- 
modulated continuous-wave optical radar [6] but it is 
larger than the expected limit for frequency-ranging tech- 
niques (about one wavelength). Our resolution was limited 
by the width of the FSF laser spectrum. This limit can be 
improved by using lasers with broad, almost flat effective 
gain (defined as the difference between frequency depen- 
dent gain and losses). Our current accuracy is at least 
5 um (assuming that the maximum can be determined 
within 10% of the width). Further improvements in band- 
width and signal-to-noise ratio should enable us to reach 
accuracy of 1 um. 

The high accuracy of the present technique makes possi- 
ble submicron resolution of surface structures whose shape 
is known, so that results can be fitted to theoretical profiles. 
To achieve such high resolution the spectrum of the FSF 
laser must be a (broad) Gaussian, to avoid signal oscilla- 
tions that originate from the Fourier transform of other, 
truncated spectral distributions. 

The present limitation on the distance range of measure- 
ments with a FSF laser, seeded by a PM signal, is the 
coherence length of the seed laser. This value has to be 
determined using the seed laser spectrum width without 
slow (on milliseconds scale) drifts. For our laser this coher- 
ence length is larger than 1 km. 
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Abstract 

We present experimental characteristics of an Yb 3+ -doped fiber ring laser operating with frequency-shifted feedback (FSF) through 
an acousto-optic modulator (AOM) and seeded by both a stationary continuous-wave (CW) laser and spontaneous emission. We show 
the spectrum and output characteristics for operations with several effective gain bandwidths, as established by Fabry-Perot etalons 
inside the cavity. Observation using a high finesse Fabry-Perot interferometer shows that, as expected from earlier work, although 
the spectrum of the FSF laser without seeding is continuous, when seeded by a CW-laser the spectrum consists of a comb of discrete 
modes, each offset from the seed by an integer number of AOM frequency shifts. The experimental results are in excellent quantitative 
agreement with the theory developed earlier [L. Yatsenko, B.W. Shore, K. Bergmann, Opt. Commun. 236 (2004) 183]. 
© 2006 Published by Elsevier B.V. 

PACS: 42.55.-f; 42.60.Da; 42.55.Ah 

Keywords: Optics; Lasers; Frequency- shifted feedback 



1. Introduction 

In conventional lasers, the light makes multiple passes 
through a gain medium inside a cavity. The cavity length 
defines discrete allowed longitudinal modes for which con- 
structive interference supports gain [1], By contrast, fre- 
quency-shifted feedback (FSF) lasers have within their 
cavity some element that introduces a frequency shift dur- 
ing each successive pass. In such a laser no single frequency 
has any special advantage for gain, and growth (lasing) can 
occur from any initial seed frequency. Consequently, the 
output spectrum from an FSF laser can cover a much 
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broader range of frequencies than do conventional lasers, 
see e.g. [2], 

The first FSF laser system was suggested in 1965 by Fros- 
ter et al. [3], who demonstrated modelocking operation of a 
He-Ne laser coupled to an external linear cavity by means 
of an acousto-optic Bragg diffraction cell driven at half of 
the laser-mode separation frequency. After interaction with 
the cavity, the frequency-shifted radiation locked all the 
longitudinal oscillating modes. A modern concept of a 
FSF laser, with an acousto-optic cell inside the cavity, 
was suggested in 1970, when Streifer et al. [4] proposed a 
new scheme using a dye laser electronically tuned by an 
acousto-optic filter. After these pioneering works, the idea 
of a FSF laser was picked up and extended by many 
authors, who suggested and investigated numerous FSF- 
laser systems, such as: a dye laser that was electronically 
tunable over 78 nm [5]; diode lasers with an external cavity; 
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a diode laser with frequency-shifted feedback that gener- 
ated output over a large frequency range [6,7]; and bulk 
solid-state FSF lasers such as Ti:sapphire [8-10], Nd:YV0 4 
[11] and many others. 

The reason for this interest is easily understood, for FSF 
lasers have many peculiar spectral features: the spectrum is 
modeless, nearly Gaussian in shape (shifted away from the 
minimum of losses), with a very weak dependence upon 
pump and seed power. The output exhibits a chirped struc- 
ture, which can cover a large, continuous range of frequen- 
cies. These features under the many applications of FSF 
lasers, such as optical frequency domain ranging [12,13], 
optical pumping [14,15], velocity control and cooling of 
atomic beams [16], and measuring the separation of widely 
spaced optical frequencies to an accuracy of ±5 kHz [29]. 

One of the most intriguing application of FSF lasers is 
to distance measurement. As has been discussed by Kasa- 
hara et al. [11], the FSF-laser output consists of periodi- 
cally generated chirped-frequency components whose 
chirp rate is faster than 10 17 Hz/s. Phase relationships are 
maintained among various discrete frequency components. 
This structure allows one to overcome the conventional 
limits on operating range and spatial resolution and has 
thereby enabled distance measurements of 3.7 km with a 
resolution of 9.4 mm [12] and 18.5 km with a resolution 
of 20 mm [13]. 

A new stage in the study of FSF lasers began with the 
appearance of fiber lasers. These have a number of advan- 
tages, such as compactness, great range of cavity length, 
simplicity of exploitation, and great amplification per round 
trip. As with the usual FSF laser, a fiber FSF laser can oper- 
ate in different regimes. The Kerr-type nonlinearity in an 
amplifying fibre may initiate the spontaneous formation 
of short pulses [17,18]. These pulses can be controlled by 
properly designing the gain-bandwidth of the fibers 
[19,20]. Fiber FSF lasers can also demonstrate continu- 
ous-wave, modeless radiation, that can be tuned over a 
broad range [21], and to broadband multiwavelength oper- 
ation; for example, the appearance of 34 lasing wavelengths 
within a bandwidth of 28 nm [22]. Injection of a narrow- 
bandwidth seeding laser into the cavity will generate a 
discrete comb of optical frequencies, separated by the fre- 
quency A associated with the total round trip time [23]. 

Fiber FSF lasers have already found many applications, 
such as measurements of group velocity [24] and polariza- 
tion-mode dispersion [25] in fibers. They are also useful for 
optical reflectometry. The first works in that direction have 
already been reported: Yoshida et al. [26] described a FSF 
laser using erbium-doped fiber as a gain medium that 
allowed distance measurements with a spatial resolution 
of 1.5 mm. Ogurtsov et al. [27] have improved this method 
and have demonstrated a ytterbium-doped fiber FSF laser 
that is able to measure meter-long distances with an accu- 
racy of 10 um. Such results support the conclusion that 
fiber FSF lasers offer opportunities for new devices that 
will enable distance measurements to be made rapidly 
and with high accuracy over a large range of distances. 



Early publications dealing with FSF-laser operation 
concentrated on characterizing the output under various 
conditions of operation [17,7,30]. As that work made clear, 
a FSF laser can exhibit a variety of output scenarios, from 
steady to pulsed, depending on controllable parameters 
such as the pump power and the frequency shift. 

Although FSF lasers have been discussed for many 
years, and many interesting properties have been recog- 
nized [32,8], only recently has a theory provided a complete 
description of the phase and intensity of the output of a 
FSF laser, as seeded both by spontaneous emission within 
the gain medium and by a monochromatic continuous- 
wave (CW) laser injected from outside the cavity [33]. 
There the spectrum of an externally seeded FSF laser is 
predicted to have some remarkable properties. The present 
paper provides experimental verification of this theory. 

We begin, in Section 2, with a description of our appara- 
tus, which includes an Yb 3+ -doped fiber ring FSF laser with 
associated pump laser, a seeding laser and diagnostic 
devices with their parameters. Section 3 presents several 
characteristics of the FSF laser without any elements that 
deliberately limit the bandwidth. Because the loss profile 
is not known, these results could not be compared directly 
with the theory. However, data taken after inserting a 
well-characterized Fabry-Perot etalon (FPE) into the cavity 
provided a clear demonstration of the theory. Section 4 
describes measurements on such a system, with and without 
an external seeding laser. To investigate the detailed struc- 
ture of the FSF-laser spectrum, we inserted into the cavity 
a narrow bandwidth filter, consisting of four FPEs. The 
data show that without external seeding the spectrum of 
the FSF is continuous, but that with such seeding it consists 
of a comb of discrete frequencies, differing by the frequency 
of the acousto-optic modulator (AOM), as predicted by the- 
ory [33]. 

2. Apparatus 

The experimental arrangement comprises four intercon- 
nected parts: (a) the FSF ring laser itself; (b) the pump laser 
that produces the gain in the active fiber; (c) the CW seed 
laser; and (d) output diagnostics. Fig. 1 shows these four 
parts and their individual components, each numbered 
for reference. 

2.7. The FSF laser 

The FSF laser itself consists of several elements, labeled 
in Fig. 1 . Gain occurs in an ytterbium-doped fiber (2) [Jena 
325sA4 IPHT] with length 0.45 m; the Yb 3+ fiber is fusion 
spliced onto a passive fiber (1) with length 1.15 m. Light 
from the fiber passes through a dichroic mirror (4), a turn- 
ing mirror (5), and an output mirror (8), whose partial 
transparency (about 5%) diverts some of the light to the 
diagnostics section. Following this mirror, most light con- 
tinues within the FSF, passing into an AOM (9). The first- 
order diffracted light from the AOM, shifted in frequency 
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diagnostics 

Fig. I. Layout of our apparatus, showing (a) the FSF ring laser, (b) the 
pump laser, (c) the seed laser, and (d) the apparatus for diagnostics. Heavy 
lines, with arrowheads, show the closed path of FSF laser light. Dotted 
lines, with arrowheads, show the pump light path. Thin lines with 
arrowheads show seed and diagnostics light paths. 

by v AO m = 80 MHz, then returns to the fiber (1) for subse- 
quent round trips. Aspherical lenses (3) and (10), with focal 
lengths 4.5 mm, collimate the light and couple the laser 
radiation into the fiber. 

In such an arrangement, after each round trip the fre- 
quency of the intracavity radiation decreases by v AO m- 
The maximum efficiency of the AOM is about 87%, so 
some 1 3% of the intracavity power is lost with each pass. 
A set of FPEs (7) are usually inserted to decrease the spec- 
tral width of the FSF-laser and to provide reproducible and 
measurable loss characteristics. 

As indicated in Fig. 1, the light in the ring laser circu- 
lates in a counterclockwise direction. An optical diode (6) 
inside the cavity diminishes the clockwise wave. The optical 
length of the ring FSF laser is 3.1 m, and its free spectral 
range (FSR) is about 96 MHz. 

2.2. The pump laser 

Excitation and consequent gain in the Yb 3 ~*"-doped fiber 
is produced by pump radiation of wavelength 980 nm from 
a fiber-pigtailed diode laser (11) [Multiplex, Inc. MTX9P], 
which enters the FSF ring path through a dichroic mirror 
(4) and thereafter couples into the Yb 3+ -doped fiber by a 
lens (3). The maximum output power of this pump laser 
is 200 mW. An optical diode (12) reduces the influence of 
back reflection of the laser radiation on the pump laser. 

2.3. The seed laser 

The external seed radiation for the present work comes 
from a single-frequency external-cavity diode laser (13), 



constructed according to the Littman-Metcalf scheme 
[34]. It consists of a laser diode, an aspheric lens, a diffrac- 
tion grating and a mirror. The spectral width of the laser 
radiation is less than 1 MHz. The seed-laser radiation 
passes through a cylindrical-lens telescope (14) and an opti- 
cal diode (15) and enters the FSF laser cavity through the 
zeroth order of the AOM (9). The maximum power of the 
seed-laser radiation on the face of the AOM was about 
1 mW. Only 2-3% of that was injected into the FSF laser 
fiber. As result, the injected seed power is much smaller 
than the intracavity laser power in the operation regime, 
and the seed laser can be considered as weak. 

2.4. The diagnostics 

Output from the FSF laser emerges through the dichroic 
mirror (8). It passes, through beamsplitters (16), (19) and 
(21), into several monitoring devices. A spectrometer (20), 
with a resolution of 50 GHz, monitors the optical spectrum 
of the FSF laser. Spectral characteristics requiring higher 
resolution are measured with the aid of two scanning 
Fabry-Perot interferometers (22) and (24), two photodi- 
odes (23) and (25) and a digital oscilloscope (26). The first 
interferometer has a FSR of 800 GHz and a finesse of about 
80, while the second one has a FSR of 7.9 GHz and a finesse 
of about 400. The time dependence of the FSF-laser output 
is monitored with a fast InGaAs photodiode (17) whose 
bandwidth is 5 GHz. The radio-frequency spectrum of the 
FSF-laser intensity is recorded by a radio-frequency spec- 
trum analyzer (18), which is sensitive over a range from 
9 kHz to 2.6 GHz. 

3. Characteristics of the unfUtered FSF laser 

Two measurable basic properties hold particular interest 
when characterizing the operation of the FSF laser. One is 
the output power, Pout* integrated over all frequencies, and 
the other is the spectral distribution of this output power. 
Both of these properties are affected by the power of the 
pump laser, by the gain (and loss) bandwidth within the 
FSF-laser cavity, and by the presence of an external seed 
laser (supplementing or competing with the ever-present 
spontaneous emission seed). Each of these factors is con- 
trollable; displays of the resulting FSF attributes - and 
comparison of these with theory - provide useful indicators 
of the FSF-laser operation. 

3.1. Output power of the unfiltered FSF laser 

A traditional diagnostic tool for characterizing laser 
operation is measurement of the total output power of a 
laser as the pump power is varied. Fig. 2 illustrates this 
dependence of output on pump input for the simplest 
configuration of the fiber FSF ring laser. For these results 
there was no external seed laser (only spontaneous emission 
within the range of frequencies where gain occurs), and no 
FPEs (element (7) in Fig. 1) to limit the gain bandwidth 
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Fig. 2. Output power of unseeded FSF laser versus pump power for 
operation without intracavity FPEs. The inset shows details of the 
threshold. Representative error bars are shown on a single element of each 
plot. 



within the ring cavity. To eliminate additional frequency- 
dependent losses caused by polarization effects in the fiber, 
we adjusted a polarization controller created by fiber coils 
to obtain a smooth laser spectrum. The results shown here 
do not depend on whether there is an injected seed. 

This curve is a typical laser input-output curve [1]. It has 
a threshold for laser operation at about 32 mW; for higher 
pump powers the FSF laser output is a linear function of 
the pump power. The inset in Fig. 2 enlarges the threshold 
region of the figure. A clear threshold exists for this laser, 
in contrast to the reported behavior of a TirSapph FSF 
laser with an intracavity filter [28], which showed a smooth 
onset of lasing rather than a well-defined threshold. 

When no intracavity filter is present, the FSF fiber laser 
generates only CW radiation for pump powers up to the 
maximum available, 200 mW.This is the preferred operat- 
ing regime for practical applications of FSF lasers to 
frequency ranging [27]. However, if some frequency selec- 
tivity of the cavity is introduced, an oscillatory regime 
exists that is typical for FSF lasers [8]: for low pump pow- 
ers there occur low frequency (about 50 kHz) pulsations of 
the output power (with 100% modulation). These oscilla- 
tory regimes have been described in detail for a Ti: sapphire 
laser [28]. 

3.2. Spectrum of the unfiltered FSF laser 

To obtain the spectral distribution of the FSF laser out- 
put, we used a scanning Fabry-Perot interferometer (FPI) 
(item (22) and (24) of Fig. 1). Fig. 3 shows the results 
obtained when the FPI had a FSR of 800 GHz and a 
finesse of 80. The figure shows typical results for the oper- 
ation of the FSF laser without external seeding or with a 
weak seed laser that cannot significantly distort the spec- 
trum of the FSF laser. When the FSF laser is not externally 
seeded, the observed spectral bandwidth of the FSF laser is 
260 GHz (full-width at half-maximum, FWHM). When the 
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Fig. 3. Spectrum of the FSF laser with (thick line) and without (thin line) 
seed laser. The arrow marks the frequency of the external seed laser. The 
difference frequency Jv, expresses the offset from the seed laser frequency. 
The pump power is 90 mW, the seed laser power is 60 uW. There are no 
selective elements in the cavity (no FPEs). 



external seed is weak, the spectrum is slightly broader and 
somewhat shifted towards the seed laser frequency. The 
width of the spectrum, T, is much smaller than the effective 
gain bandwidth determined jointly by the Yb 3+ -doped fiber 
gain line and by dispersion of the AOM. We estimate the 
effective gain bandwidth to be about 1 0 4 GHz. 

As will be noted in the following sections, the spectral 
narrowing and the asymmetric spectral profile are charac- 
teristic of FSF lasers, and are in accord with recent theory 
[33]. 

3.3. Dependence of the spectral width of the unfiltered FSF 
laser on pump power 

When the FSF ring laser operates without any band- 
width-narrowing filters (the FPE (7) of Fig. 1), the spectral 
width is relatively insensitive to the pump power. Fig. 4 
shows the dependence of the width f and the asymmetry 
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Fig. 4. Width T (□) and asymmetry parameter T 3 (•) of unfiltered FSF 
laser spectrum [FWHM] versus pump power. 
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parameter T 3 on the pump power. These parameters were 
obtained by fitting the experimental spectrum to the asym- 
metric Gaussian 



Py(Av) = P 0 exp 



(Av-Av 0 ) 2 (Av-Avpf 



1 



(i) 



The spectral width increases only slightly as the pump 
power increases from threshold (32 mW) to 160 mW. 



4. Characteristics of the filtered FSF laser 

The operation of a seeded FSF laser can be understood 
as follows. The intracavity field consists of many compo- 
nents equidistantly separated in frequency. The first one 
is the CW field of the seed-laser radiation coupled into 
the cavity. The frequency of the seed laser is normally set 
very close to the minimum of the frequency-dependent 
intracavity losses. The frequency of each subsequent com- 
ponent is larger by the AOM frequency. The amplitude 
of the component is a result of equilibrium between four 
competing processes: 

(1) income from the adjacent component that has been 
frequency shifted by the AOM, 

(2) gain, saturated by the total intensity of all spectral 
components, 

(3) frequency dependent losses, (those present intrinsi- 
cally and those added deliberately as spectral filters), 
and 

(4) loss due to AOM frequency shifting to the next 
component. , 

The self-consistent theory of the dynamics of the FSF 
laser [33] predicts the formation of a stationary spectrum 
in which the total intensity does not fully saturate the gain. 
Thus, starting from the initial seed frequency, the amplitude 
of each component increases as long as the gain is larger 
than the losses. At some frequency the gain will no longer 
suffice to overcome the losses. Beyond this frequency the 
component amplitudes will decrease. Thus, the output spec- 
trum will exhibit a maximum. 

A similar picture applies when spontaneous emission 
provides the only seed. We can then imagine that the field 
in the cavity results from incoherent seed sources with dif- 
ferent frequencies continuously distributed near the maxi- 
mum of the effective gain profile (the difference between 
gain and frequency-dependent losses). In this case the spec- 
trum is predicted [33] to be continuous with a shape similar 
to that observed with a seed laser. 

The position of the spectral maximum, and the width of 
the output spectrum, are both predictable if one knows the 
spectral characteristics of the gain and loss mechanisms 
affecting the FSF laser. These are difficult to measure accu- 
rately - indeed, they vary irregularly and uncontrollably - 
but they can be controlled and accurately measured by 



introducing well-characterized spectral filters. Towards 
that end, we investigated a FSF-laser cavity containing a 
relatively broadband filter, created from FPEs (element 
(7) of Fig. 1). The effective FSR was 1 THz. The following 
paragraphs describe the observed properties of that system. 

4. 1. Loss profile and output spectrum of the filtered FSF 
laser 

4.1.1. Broadly filtered FSF laser 

Fig. 5a shows the measured relative transmission profile 
(dots) of the FSF cavity when a broadband FPE was pres- 
ent to act as a filter. The full line is a Gaussian fit to these 
values; this is an excellent approximation. 

Fig. 5b shows the measured output spectrum of the FSF 
laser in more detail, as recorded using a scanning FPI with 
FSR of 800 GHz and finesse of 80. The vertical dashed line 
marks the maximum of the output spectrum. Physical intui- 
tion, validated by theory [33], predict that this maximum 
should occur where the saturated gain equals the losses 
from a cavity round trip [33]. This behavior allows us to find 
the full width of the frequency interval in which the saturated 
gain exceeds losses. We.measured this width to be 180 GHz. 
It is significant that Fig. 5, which is entirely experimental, is 
basically the same as Fig. 2 from the theory paper [33]. As 
predicted by theory [33], the experimental spectrum is close 
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Fig. 5. (a) Transmission profile of broad filter, and (b) spectrum of seeded 
FSF laser with this filter inside the cavity. The detuning Av is defined as 
the difference between the observation frequency and the frequency of 
maximum filter transmission. 
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to a Gaussian near the maximum, but shows asymmetry in 
the wings of the profile. As also predicted, its maximum coin- 
cides with the zero of the effective gain. 

4.1.2. Spectra of a narrowly filtered FSF laser with and 
without a seed laser 

The only difference between the experimental results 
shown in Fig. 5 and the theory presented in Fig. 2 of 
Ref. [33] is that the experimental spectrum does not reveal 
the discrete structure of the output spectrum of a seeded 
FSF laser. This is predicted to consist of a multitude of dis- 
crete components, each separated by the frequency incre- 
ment v AO m = ^/2tc = 80 MHz of the AOM. To resolve 
these components it is necessary to limit the bandwidth 
of the cavity by inserting a very narrow bandwidth filter. 
Such filtering leads to a laser spectral width smaller than 
the FSR of 7.9 GHz for the available FPI. We accom- 
plished this filtering by inserting into the cavity four FPEs, 
with FSR 1000, 200, 50 and 15 GHz. The resulting spectral 
output has a bandwidth of 2.8 GHz, in contrast to the 
60 GHz bandwidth shown in Fig. 5. 

The lower frame of Fig. 6 shows the output spectrum of 
the narrow-band FSF laser, as recorded using a scanning 
FPI with a FSR of 7.9 GHz and a finesse of 400. Arrows 
mark the spectral bandwidth (FWHM) of 2.8 GHz. The 
individual discrete spectral components expected from the- 
ory are clearly visible. For comparison, a solid line shows 
the spectrum observed when there is no external seed laser; 
as expected, the latter is a continuum, with a bandwidth of 
2.7 GHZ. 

The upper frame of Fig. 6 provides an expanded view of 
the individual spectral components. Vertical lines, rising 
from the lower border, mark increments of the AOM fre- 
quency. Fiducial marks along the top border occur at the 
spacing expected for longitudinal modes in a cavity whose 
length equals that of the FSF laser. The observed spectral 
components do not have the frequency spacing of longitu- 
dinal modes; they occur at multiples of the AOM fre- 
quency, as predicted. 

It is worth noting that the experimental evidence of a 
discrete stationary spectrum, comprising a set of frequen- 
cies separated by the AOM frequency, does not conflict 
with the well-documented moving-comb model of a FSF 
laser [11]. In the moving-comb model the output of a 
FSF laser is a set of discrete frequencies ("teeth"), whose 
separation is the cavity axial mode interval A SLX = 2tcv lm , 
and which shift steadily in frequency at the rate y c = A/x T 
where x T is the round-trip time. In [33] we proved that, in 
a FSF laser without fluctuations, both models are equiva- 
lent. This becomes immediately obvious if one notices that 
the moving comb model [11] describes a process that is 
periodic in time, with the period 2n/A = 1/vaom- This 
means that the process can be considered as a discrete Fou- 
rier series with the frequency difference A between compo- 
nents, as fixed by the AOM. Although the FSF laser, when 
seeded by CW radiation, is periodic, it will not be periodic 
when the field grows from a short pulse. 



(a) 



Av 8 /v LM 




0.0 -0.5 -1.0 -1.5 -2.0 -2.5 -3.0 -3.5 -4.0 
Av B [GHzJ 

Fig. 6. Spectrum of the FSF laser observed through a narrow bandwidth 
filter: the power spectral density P 9 (arbitrary units) versus Jv„ the 
frequency offset from the frequency of seed laser, (a) expanded view of (b), 
with Av t in units of the AOM frequency. Thin lines are with external seed 
laser, thick lines are without seed laser. 



To display a moving comb of frequencies one needs a 
spectral recording device with bandwidth 8a> smaller than 
the cavity axial mode interval A ax (to resolve adjacent com- 
ponents). The signal integration time l/5co has to be 
shorter than the time interval \/A between the occurrence 
of teeth near the observation frequency. Thus, the band- 
width 5co must satisfy the condition A < bco < A ax . Because 
in our experiment we have the condition A > A ax , we do not 
see any comb-like structure. In [31] the moving comb has 
been observed for a very small FSF frequency (1 MHz or 
less) realized by shifting 80 MHz up by one AOM and 
79 MHz down by a second AOM. 

These restrictions on the bandwidth requirements are 
relevant for displaying a moving comb from which one 
could directly measure the changing positions of frequency 
components, i.e. the chirp. Alternatively, one can use a 
Michelson interferometer and interpret the results to infer 
a chirped frequency of a moving comb with fluctuating 
parameters, as was done by Kasahara et al. [1 1]. We discuss 
this further in a forthcoming paper submitted to Optics 
Communications. 
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Because the output bandwidth is important for applica- 
tions, it is useful to estimate the frequency range over 
which a seeded comb can extend. For a broadly filtered 
(and especially for an unfiltered) FSF laser, direct observa- 
tion of a stationary comb is impossible with our scanning 
FPI for any total laser output bandwidth that is larger than 
the FSR of the FPI. Nevertheless, we have observed [27] 
that, when using a technique requiring coherence between 
comb components, the distance resolution is directly pro- 
portional to the bandwidth T of the FSF laser (see Fig. 7 
in [27]). This observed linearity demonstrates that the 
FSF laser maintains coherence of the acousto-optic modes 
across the entire spectrum of the FSF laser, right up to the 
spectrum width of about 3500 GHz. This corresponds to 
about 80,000 frequency components (comb teeth). 

4.2. Output power of the spectrally filtered FSF laser 

Fig. 7a shows that the output power from the FSF laser 
depends upon the pump power when a broadband FPE is 
inserted into the ring cavity but no external seed laser is 
present. In contrast to the behavior of the unfiltered FSF 
laser (see Fig. 2), the output of the FSF laser .with the 
FPE spectral filter in place exhibits a regime of CW opera- 
tion (labeled C, this exists for P Q < P out < Ptj as the pump 
power increases), a regime of low-frequency oscillation 
(labeled O, this exists for P oxlt < P 0 as the pump power 
increases and for P om < P& as the pump power decreases), 
and a third regime, one characterized by modelocking 
(labeled M, this exists for P out > P\j as the pump power 
increases and for P oul > P D as the pump power decreases). 
Unlike the unfiltered laser, for which a CW regime exists 
for all available pump power (up to 200 mW), the FSF 
laser exhibits mode locking when a broadband intracavity 
FPE filter is in place and the pump power is high. (When 
the filter bandwidth is narrow, no FSF modelocking regime 



exists.) The transition to mode locking is apparent in the 
output when viewed by the radio-frequency (RF) analyzer. 
Fig. 8 shows the RF analyzer signal near the axial mode 
frequency of the ring cavity (about 95.88 MHz). In the 
CW regime of FSF laser operation, there is no axial mode 
beating. However, near the frequency corresponding to an 
axial-mode difference there occurs a resonance increase in 
the spectral density of laser output power fluctuations. This 
resonance arises because any fluctuation in the laser intra- 
cavity power persists much longer than a round-trip time, 
thereby contributing periodically to the output power fluc- 
tuations. The width of the resonance is determined by the 
fluctuation life-time. This picture changes completely in 
the mode-locking regime: the amplitude of the signal grows 
by more than two orders of magnitude and the beating sig- 
nal is almost monochromatic (the width of the spectrum in 
Fig. 8(b) is determined by the apparatus function). The 
narrowing of the RF spectrum is clear evidence that, in 
the modelocking regime, there exist frequency components 
separated by the axial mode frequency. We observe a signal 
arising from the beating between these phase-synchronized 
components. The time duration of the pulses generated in 
the modelocking regime is determined, as usual for fluctu- 
ation-free pulses, by the total laser spectral bandwidth. It is 
about 10 ps. 

Unlike behavior reported with a Ti:Saph FSF laser [8], 
the fiber FSF exhibits no discontinuity in the spectrally 
integrated total power at the onset of modelocking. For 
an unseeded laser, the interval of pump powers where the 
CW regime exists is quite narrow (from P 0 — 42 mW to 
P v ~ 49 mW). The modelocking behavior shows a pro- 
nounced hysteresis, evident in Fig. 7. As the pump power 
increases, mode locking begins when the pump power is 
Pu, but when the power is decreasing the modelocking 
ceases when the pump power is P D . For pump power less 
than P& the CW regime does not exist at all. 




0 20 40 60 80 100 120 140 160 180 0 20 40 60 80 100 120 140 160 180 

P p [mW] 

Fig. 7. Output power P oul of the FSF laser versus pump power P ¥ when (a) the cavity includes a broadband filter and no seed laser, and (b) with external 
seed-laser. 
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Fig. 8. Radio frequency analyzer signal near the axial mode difference 
frequency of the ring cavity (about 95.88 MHz). Frame (a) is for the CW 
regime, and frame (b) is for the self mode-locking regime. 



4.3. Seeded self mode- locking 



nounced when its frequency is 30 GHz higher than the fre- 
quency of minimum loss. 

4.3.2. Dependence on the seed-laser power 

As the seed power increases, the hysteresis changes and 
the pump powers P v and P D , for which self modelocking 
begins and ends, each increase. Fig. 10 shows the depen- 
dence of P\j and upon the seed laser power. Horizontal 
arrows mark the values observed when there is no seed. 
Thus, the seed laser suppresses self modelocking. This 
shows that the self modelocking regime is qualitatively dif- 
ferent than pulsed operation of a seeded FSF laser consid- 
ered in [33]. 

4.3.3. Spectra of the FSF laser in the self-modelocked regime 
Fig. 11 illustrates the behavior of self mode-locking in 

the fiber FSF laser. A large shift and distortion of the spec- 
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Fig. 9. Pump powers Pu and P 0 corresponding to beginning (□) and 
ending (•) of self modelocking versus detuning Av s0 of the seed laser 
frequency from the cavity loss minimum. Arrows mark the laser operation 
threshold and the values obtained when there is no seeding. 



4.3.1. Dependence on the seed- laser frequency 

The operation of the filtered FSF laser depicted in Fig. 7a 
relies on signal growth from spontaneous emission. The use 
of an external seed laser provides better control over the 
output characteristics, evidenced in the values of Pu and 
P a where self modelocking begins and ends. The seed laser 
increases thresholds of the beginning and the end of the 
mode-locking interval, thereby enlarging the regime of 
"useful" CW operation (see Fig. 7). 

Fig. 9 shows the pump powers Pu and Pu correspond- 
ing to the beginning and ending of self modelocking as a 
function of the detuning of the seed laser frequency from 
the frequency where the cavity loss is minimum (the central 
frequency of the inserted FPE). The thresholds P v and P D 
are reproducible, as indicated by the size of the error bars 
in the figure. 

For reference, the figure also indicates, by arrows, the 
threshold for laser operation and the values obtained when 
there is no seed. The influence of the seed laser is most pro- 




threshhold 
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Fig. 10. Pump powers /\j and Pn corresponding to the beginning (□) and 
ending (•) of self mode-locking versus power of the seed laser. Arrows 
mark the thresholds without external seed and the laser operation 
threshold. 
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trum are observed when mode-locking occurs. This con- 
trasts with the mode-locking behavior in a Ti:Saph FSF 
laser [8], There the spectrum in the mode-locking regime 
underwent only a small frequency shift while maintaining 
the same shape. The spectrum in Fig. 1 1 shows that two 
regimes (CW and mode-locking) seem to coexist in our 
laser, localized in different spectral areas. Mode-locking 
exists near the minimum of intracavity losses, whereas 
the spectrum of the CW regime is shifted to lower frequen- 
cies according to the sign of the AOM shift. Observations 
of the spectral changes occurring when the external seed 
laser is turned on (intermediate line on Fig. 1 1) confirm this 
interpretation. The laser power is redistributed in such a 
way that the CW-regime spectrum increases and mode- 
locking diminishes significantly. This behavior confirms 
the conclusion that the seed laser suppresses self mode- 
locking by promoting a redistribution of laser power from 
mode-locking to CW spectral regions. 

4,4, Comparison with theory 

Our earlier work [33] presented a theoretical description 
of the spectral characteristics of the FSF laser seeded both 
by spontaneous emission and a monochromatic external 
seed laser. That theory predicts that, apart from a small 
asymmetry, the FSF laser spectrum is accurately approxi- 
mated by a Gaussian, 



P v {Av) = P 0 exp 



( _ [(Jv - cop) 2 ] \ 



(2) 



where Av is the detuning of the observation frequency from 
the minimum of the filter losses. This Gaussian form de- 
scribes both the continuous spectrum of the laser seeded 
by spontaneous emission and the envelope of the comb 
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Fig. 11. Spectrum of the FSF laser when the broadband filter is inserted 
into the cavity: the thin line is valid for CW operation (regime (C)) without 
an external seed laser; the intermediate line is the self mode-locking 
(regime (M)) with a seed laser; the thick line is the self mode-locking 
(regime (M)) without a seed laser. The detuning Jv is the difference 
between the observation frequency and the frequency of maximum filter 
transmission. 



of discrete components for a FSF laser seeded by a mono- 
chromatic laser. Our previous work [33] presented simple 
formulas for the spectral width T and the shift a> 0 of the 
spectral maximum from the minimum of filter losses. For 
a FSF laser without external seed, the relevant formulas are 

y 



1/3 



(3) 
(4) 

(5) 
(6) 



oj 0 = ay, 
where 

and 

"(H) 

The parameters characterizing the spectrum depend 
mainly on the precisely known AOM frequency and on 
the filter width, parameterized as T f . As mentioned above, 
quantitative comparison of experimental results with the- 
ory prediction requires the use of well-characterized spec- 
tral intracavity filters. With such a filter emplaced, the 
parameter Ff is obtained by fitting the measured variation 
of filter transmission versus frequency to the theoretical 
curve. 

The parameter <x, determined by the small parameter 
Ps P <gi 1, characterizes spontaneous emission. In principle 
p sp can be estimated from characteristics of the fiber and 
the laser cavity. However, it is more convenient to derive 
the parameter a from Eqs. (3) and (4) as 



/coo\ 2 / 3 



(7) 



thereby allowing estimation from experimental data for the 
width and shift of the laser spectrum. 

Figs. 12 and 13 show measured values of the width T and 
frequency of spectral maximum co 0 versus pump power. 
Filled and empty points refer to laser operation with and 
without an external seed laser. Different combinations of 
intracavity etalon filters were used: (1) one etalon (FSR 
1000 GHz), (2) two etalons (FSR 1000 and 200 GHz) and 
(3) three etalons (1000, 200 and 50 GHz). The measurement 
errors are less than the symbol sizes. Table 1 gives the values 
of r f obtained, as described before, by fitting the measured 
variation of filter transmission versus frequency to the the- 
oretical Gaussian profile The table also collects calculated 
and measured parameters for laser operation without an 
external seed laser. Using the mean value of the parameter 
a — 1.98, we calculated the shift co 0 and width r. To take 
into account the apparatus width of the FPEs (5 GHz for 
cases 1 and 2, and 0.5 GHz for case 3), we evaluated a con- 
volution of a Gaussian spectrum with the apparatus func- 
tion and fitted it to a Gaussian with a width r ca i c . The 
straight lines in Figs. 12 and 13 mark the theoretical values 

Of r C alc and COo.calc- 

Table 1 and Figs. 12 and 13 show excellent quantitative 
agreement between theory and experiment when there is no 
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Fig. 12. The spectral width of the FSF laser r versus pump power P P for 
different combinations of intracavity etalon filters: one etalon (FSR 
1000 GHz, square points), two etalons (FSR 1000 and 200 GHz, round 
points) and three etalons (1000, 200 and 50 GHz, triangular points). Solid 
points: external seed is present; open points: no external seed. The straight 
lines show the calculated values of 



90 
80 
70 
60 

— • 50 
X 

O 40 

-3* 30 - 
20 
10 



(j O OQOOO < 



30 40 50 60 70 



80 90 100 110 120 130 140 150 
P p [mW] 



Fig. 13. The shift (o 0 of the spectral maximum of the FSF laser with 
respect to the minimum of the overall losses versus pump power P P for one 
etalon (FSR 1000 GHz, square points) and two etalons (FSR 1000 and 
200 GHz, round points). Solid points: external seed is present; open 
points: no external seed. The straight lines show the calculated ct> 0lC aic- 



Table 1 

Experimental and theoretical parameters 
Case 1 



FSR 


1000 GHz 


200 GHz 


50 GHz 


r f 


867 GHz 


120 GHz 


41 GHz 


y 


39.2 GHz 


10.5 GHz 


5.1 GHz 


a 


1.97 


1.99 




r caJc 


29.9 GHz 


7.64 GHz 


3.82 GHz 


F mcas 


29 GHz 


7.5 GHz 


3.6 GHz 




77.2 GHz 


20.8 GHz 


10.1 GHz 




80 GHz 


21 GHz 





external seeding. When an external seed is present, there is 
competition between spontaneous emission and seeding 
that leads to a slight, but still observable, dependence of 
f and co 0 upon pump power. 



5. Summary and conclusions 

This paper presents detailed properties of a FSF laser 
constructed in a ring configuration, using an Yb 3+ -doped 
fiber as the gain medium and an acousto-optic modulator 
as the frequency shifter. 

The results demonstrate the absence of longitudinal 
modes for the FSF laser when it operates without any 
external seed. When a monochromatic seed is present, the 
output spectrum exhibits the predicted discrete frequencies 
built by successive AOM increments upon the seed fre- 
quency. The envelope of these frequency components fol- 
lows the predicted analytic form, rising to a maximum at 
the predicted frequency. The peak frequency and the width 
of the spectrum are essentially independent of pump 
power, as theory predicts. 

Examples of the spectra and output power characteris- 
tics with different intracavity selective filters have been pre- 
sented here. These demonstrate that, by increasing the 
power of the seed laser, it is possible to increase the thresh- 
old power where mode-locking begins. Such an increase 
enlarges the CW-regime of the output radiation - the "use- 
ful" regime for application to distance measurements. 

By using an intracavity filter with known loss character- 
istics, we have been able to calculate important laser param- 
eters. Using only a single parameter, obtained by fitting 
observations, we are able to predict all other parameters 
in excellent agreement with experiment.The results are in 
perfect agreement with our earlier theory [33]. This detailed 
agreement with experiment gives us confidence that, within 
its expected domain of applicability, the theory provides a 
satisfactory basis for predicting FSF-laser properties. 
Recent work provides a detailed description of an applica- 
tion of this FSF laser to distance measurements [27]. 

Acknowledgements 

We acknowledge support by the Stiftung Rehinland- 
Pfalz fur Innovation. L.P.Y. acknowledges support by 
the Deutsche Forschungsgemeinschaft (436-UKR-l 13/ 
16). B.W.S. and L.P.Y. acknowledge support from the 
Max Planck Forschungspreis 2003. 

References 

[1] A.E. Siegman, Lasers, Univ. Sci. Books, Mill Valley, CA, 1986; 

P. Milonni, J.H. Eberly, Lasers, Wiley, New York, 1988. 
[2] F.V. Kowalski, S.J. Shattil, P.D. Hale, Appl. Phys. Lett. 53 (1988) 

734; 

K.. Nakamura, K. Kasahara, M. Sato, H. I to, Opt. Commun. 121 
(1995) 137; 

K. Nakamura, F. Abe, K. Kasahara, T. Hara, M. Sato, H. Ito, IEEE 
J. Quant. Electron. 33 (1997) 103; 

I. CM. Littler, S. Balle, K. Bergmann. Opt. Commun. 88 (1992) 514; 

I. CM. Littler, S. Balle, K. Bergmann, J. Opt. Soc. Am. B8(I991) 1412. 
[3] L.C Foster, M.D. Ewy, C.B. Crumly, Appl. Phys. Lett. 6 (1965) 6. 
[4] W. Streifer, J.R. Whinnery, Appl. Phys. Lett. 17 (1970) 335. 
[5] D.J. Taylor, S.E. Harris, S.T.K. Nieh, T.W. Hansch, Appl. Phys. 

Lett. 19 (1971) 269. 



V. V. Ogurtsov el ai I Optics Communications 266 (2006) 627-637 



637 



[6] A. P. Willis, AJ. Ferguson, D.M. Kane. Opt. Commun. 116(1995)87. 
[7] K.A. Shore, D.M. Kane, IEEE J. Quant. Electron. 35 (1999) 1053. 
[8] G. Bonnet, S. Balle, Th. Kraft, K. Bergmann, Opt. Commun. 123 
(1996)790. 

(9] J. Geng, S. Wada, Y. Urata, H. Tashiro, Opt. Lett. 24 (1999) 676. 
[10] J. Geng, S. Wada, N. Saito, H. Tashiro, Opt. Lett. 24 (1999) 1635. 
[11] K. Kasahara, K. Nakamura, M. Sato, H. Ito, IEEE J. Quant. 

Electron. 34(1998) 190. 
[12] K. Nakamura, T. Miyahara, M. Yoshida, T. Hara, H. Ito, IEEE 

Photon. Techn. Lett. 10 (1998) 1772. 
[13] K. Nakamura, T. Hara, M. Yoshida, T. Miyahara, H. Ito, IEEE J. 

Quant. Electron. 36 (2000) 305. 
[14] M.J. Lim, CI. Sukenik, T.H. Stievater, P.H. Bucksbaum, R.S. Conti, 

Opt. Commun. 147 (1998) 99. 
[15] M. Cashen, V. Bretin, H. Metcalf, J. Opt. Soc. Am. B 17 (2000) 530. 
[16] I. CM. Littler, H.M. Keller, U. Gaubatz, K. Bergmannn, Z. Phys. D 

18 (1991) 307. 

[17] H. Sabert, E. Brinkmeyer, Electron. Lett. 29 (1993) 2122. 

[18] H. Sabert, E. Brinkmeyer, J. Lightwave Technoi. 12(1994) 1360. 

[19] O.G. Okhotnikov, Electron. Lett. 34 (1998) 1493. 

[20] J.M. Sousa, O.G. Okhotnikov, Opt. Commun. 183 (2000) 227. 

[21] S.H. Yun, D.J. Richardson, D O. Culverhouse, B.Y. Kim, IEEE J. 

Sel. Top. Quant. Electron. 3 (1997) 1087. 
[22] S.K. Kim, M.J. Chu, J.H. Lee, Opt. Commun. 190 (2001) 291. 
[23] I.R. Perry, R.L. Wang, J.R.M. Barr, Opt. Commun. 109 (1994) 

187. 

[24] M. Yoshida, K. Nakamura, H. Ito, IEEE Photon. Techn. Lett. 13 
(2001) 227. 



[25] M. Yoshida, T. Miyamoto, N. Zou, K. Nakamura, H. Ito, Opt. Exp. 
9(2001) 207. 

[26] M. Yoshida, K. Nakamura, H. Ito, Rev. Laser Eng. 27 (1999) 
490. 

[27] V.V. Ogurtsov, L.P. Yatsenko. V.M. Khodakovskyy, B.W. Shore, G. 

Bonnet, K. Bergmann, Opt. Commun., in press. doi:I0.10!6/ 

j.optcom.2006.04.070. 
[28] M. Stellpflug, G. Bonnet, B.W. Shore, K. Bergmann, Opt. Exp. 1 1 

(2003) 2060. 

[29] J.R.M. Barr, G.Y. Liang, M.W. Phillips, Opt. Utt. 18 (1993) 1010. 
[30] F. Kowalski, J.A. Squier, J.T. Pinckney, Appi. Phys. Lett. 50 (1987) 
711; 

F.V. Kowalski. P.D. Hale, S.J. Shattil, Opt. Lett. 13 (1988) 622; 

P.D. Hale, F.V. Kowalski, IEEE J. Quant. Electron. 26 (1990) 1845; 

P.I. Richter, T.W. Hansen, Opt. Commun. 85 (1991) 414; 

M.W. Phillips, G.Y. Liang, J.R.M. Barr, Opt. Commun. 100 (1993) 

473; 

J. Martin, Y. Zhao, S. Balle, K. Bergmann, M.P. Fewell, Opt. 

Commun. 112 (1994) 109. 
[31] S. Balle, I. CM. Littler, K. Bergmann, F.V. Kowalski, Opt. Commun. 

102 (1993) 166. 
[32] C.C. Cutler, IEEE J. Quant. Electron. 28 (1992) 60; 

I.C.M. Littler, J.H. Eschner, Opt. Commun. 87 (1992) 44; 

S. Balle, K. Bergmann, Opt. Commun. 116 (1995) 136; 

K. Kasahara, K. Nakamura, M. Sato, H. Ito, Opt. Rev. 4 (1997) 180. 
[33] L. Yatsenko, B.W. Shore, K. Bergmann, Opt. Commun. 236 (2004) 

183. 

[34] M.G. Littman, H.J. Metcalf, Appl. Opt. 17 (1978) 2224. 




Available online at www.sciencedirect.com 



*o" ScienceDirect r^.^jMirAi 



Communications 

ELSEVIER Optics Communications 281 (2008) 1679-1685 — 

www.elsevier.com/locate/optcom 



An all-fiber frequency-shifted feedback laser for optical ranging; 

signal variation with distance 

V.V. Ogurtsov a , V.M. Khodakovskyy a , L.P. Yatsenko 3 , B.W. Shore b '*, 

G. Bonnet b , K. Bergmann b 

a Institute of Physics, Ukrainian Academy of Sciences, Prospect Nauki 46, Kiev-39, 03650, Ukraine 
b Technische Universitdt Kaiserslautern. 67653 Kaiserslautern, Germany 

Received 14 November 2007; accepted 17 November 2007 



Abstract 

A recent paper [L.P. Yatsenko et al., Opt. Commun. 242 (2004) 581] provided a first-principles prediction for the optical ranging sig- 
nals obtained when using a frequency- shifted feedback (FSF) laser system, seeded by a phase-modulated laser. Such a system has many 
useful advantages over other alternative FSF laser techniques. We report here experimental verification of that theory, specifically the 
variation of the amplitude modulation signal with both distance and modulation index of the seed laser. We describe the operation of an 
all-fiber FSF laser that uses an Er 3+ -doped active fiber as the gain medium. To improve the signal and minimize the noise we seed the 
FSF laser with a phase-modulated (PM) laser; the measurement of distance derives from a measurement of amplitude modulation within 
a narrow frequency interval. We demonstrate that the resulting system is capable of fast and precise measurements. With the bandwidth 
limitations of our current system we achieved an accuracy better than 0.1 mm. Although measurements based on interferometry offer the 
potential for much greater accuracy under carefully controlled conditions, the present method does not suffer from the presence of a 
material-dependent phase shift at the surface of the measured object. 
<g> 2007 Elsevier B.V. All rights reserved. 

PACS: 42.55.-f; 42.60.Da; 42.55.Ah 

Keywords: Lasers; Laser ranging; Frequency-shifted feedback 



1. Introduction 

A frequency-shifted feedback (FSF) laser is a special 
type of laser with a frequency shifter inside the cavity (it 
can be an acousto-optic modulator (AOM), a moving dif- 
fractive grating, etc.) [1-9]. The presence of such an ele- 
ment that shifts radiation frequency by a fixed value 
during each cavity round trip provides FSF lasers with sev- 
eral unusual properties not found with ordinary lasers. 
Most notably, the spectrum of a FSF laser has no well- 
defined longitudinal modes because the frequency shifter 
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destroys the constructive interference inside the cavity. This 
makes a FSF laser a potential source of wideband contin- 
uous radiation, a property that has already found many 
applications [10-14]. Another important feature of a FSF 
laser is that, because of periodic frequency shifting on each 
round trip, the output radiation can be regarded as fre- 
quency swept (chirped) [6,15]. After passing through a 
Michelson interferometer any chirped laser will produce 
an amplitude modulation at a frequency proportional to 
the difference in length of the two interferometer arms. This 
makes the FSF laser a promising tool for optical ranging of 
distances up to several kilometers [16,17]. 

As describe earlier [24], optical ranging with a FSF laser 
can be greatly improved by using an external phase-modu- 
lated monochromatic laser to seed the FSF laser, thereby 



1680 



V. V. Ogurtsov el all Optics Communications 281 (2008) 1679-1685 



improving the signal detection and lessening the noise. This 
method was demonstrated experimentally [18] with an 
accuracy better than 10 um. 

Ranging devices based on all-fiber FSF lasers hold 
promise for practical application because of their robust 
and compact construction. The use of an Er 3+ -doped active 
fiber as a gain medium offers particular advantages because 
their central wavelength, 1.55 um, is used in telecommuni- 
cation and therefore there are many readily available inex- 
pensive fiber elements and devices for that wavelength. 
Another advantage of erbium lasers, very important for 
any industrial application, is their safety for eyes. 

Our technique, like those based on interferometry, has at 
first glance an intrinsic ambiguity regarding the deduction 
of a length from a time or frequency: the measurements only 
give a result modulo an increment. Interferometric methods 
commonly use the technique of an effective wavelength, 
based on the use of two wavelengths or a phase modulation, 
to resolve the ambiguity. We have noted appropriate 
approaches to overcome this problem. 

The present paper describes the experimental character- 
istics of an all-fiber FSF laser with Er 3+ -doped active fiber 
as the gain medium. It is seeded by a phase-modulated 
(PM) laser for precise measurements of distances. Using 
this laser we are able to confirm the recent first-principles 
derivation [24] of the dependence of signal strength upon 
distance; the predicted Bessel function behavior is clearly 
seen. 

2. Er 34 " -doped fiber FSF laser 

As shown in Fig. 1 , the all-fiber FSF laser consists of a 
fiber cavity (optical length 19.2 m), into which is inserted 
an Er 3+ -doped active fiber (5 m long) pumped by a diode 
laser of wavelength 980 nm) through a wavelength division 
multiplexer (WDM). The frequency-shifting element is an 
acousto-optic modulator (AOM) that upshifts the intracav- 
ity radiation by 80 MHz with each round trip. Output and 
seed insertion occurs through two 10% couplers. We use 
two couplers instead of one to avoid placing 90% of the 




pump 
laser 



-Off) 



active fiber 



output j 




Fig. I. Layout of the all-fiber FSF laser. Light from a pump laser enters 
the fiber circuit through a WDM (top left) and provides gain in an Er 3+ - 
doped active fiber (shown as loops). Within the circuit the FSF light 
travels counterclockwise, starting from an external seed laser (lower right) 
entering the circuit through a coupler that transmits 10% of this light. An 
AOM induces a frequency upshift of 80 MHz prior to amplification in the 
active fiber. A second coupler diverts 10% of the light into output for 
observation. 
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Fig. 2. Power characteristics of the all-fiber FSF laser with and without 
the seed laser: output power P om , in mW, as a function of pump power 
^pump. also in mW. The insert shows the threshold region. The output 
exhibits mode-locking, with consequent formation of pulse trains, when 
the pump power is sufficiently large. The arrows show the beginning of the 
mode-locking regime with and without external seed: the presence of the 
seed increases the power at which mode-locking starts. The useful regime 
for ranging is to the left of these arrows. 



seed laser power into the output radiation [19]. All fiber 
components are joined by a fusion splice. 

The output characteristics of our FSF laser are typical 
(see Fig. 2): the output power increases linearly with the 
pump power above a threshold of about 1 5 mW pump. A 
common feature of FSF lasers [5,19] is that for sufficiently 
high pump power a self mode-locking regime usually exists. 
Such pulse-train output is not suitable for optical ranging. 
In our laser, in the absence of the seed laser, mode-locking 
appears at quite low pump power (about 60 mW). Injection 
of an external CW seed greatly increases the threshold 
pump power for self mode-locking from about 60 mW to 
about 170 mW, thereby extending the CW- regime that is 
suitable for distance measurements [19]. 

The spectrum of this laser with and without external 
monochromatic seed is a slightly asymmetric Gaussian 
(Fig. 3) with FWHM 150-160 GHz. The FSF laser spec- 
trum without a seed laser is continuous whereas with a seed 
it consists of equidistant components (unresolved in Fig. 3) 
separated by the AOM frequency [15,19,20]. 

3. Ranging with a FSF laser seeded by a phase-modulated 
external laser 

The output radiation from a FSF laser can be described 
by a so-called "moving comb" model [15,16,20-23]. In this 
model a FSF laser, characterized by a round-trip time r r 
and an AOM frequency A, has an output spectrum compris- 
ing a very large number (103 or more) of equidistant fre- 
quencies (comb teeth), separated by the cavity axial mode 
interval A c = 2n/r r , and moving with the rate y = A/r r , 

(o n (t) = aw + yt - nA c . (1) 
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Fig. 3. Spectrum of the all-fiber FSF laser with (thick black Une) and 
without (thin red line) the external seed laser. The arrow marks the 
frequency of the seed laser. The presence of the seed introduces some 
distortion of the spectrum, but it remains reasonably well-described by a 
Gaussian. (For interpretation of the references to color in this figure 
legend, the reader is referred to the web version of this article.) 



Here, as in previous work, the constant co max is the fre- 
quency at which the time varying component has maxi- 
mum amplitude. Fig. 4aa portrays, at a fixed time /, this 
frequency comb. Such a field, when delayed by time incre- 
ment T, will again appear as a frequency comb, but offset 
(at time /) by a frequency increment yT (see Fig. 4b). 

For application to ranging these two fields are those of 
the two arms of a Michelson interferometer whose arms 
differ in length by L. The resulting time delay T = 2L/c 
produces a set of difference frequencies, from combinations 
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Fig. 4. (a) A portion of the sequence of equidistant discrete FSF-output 
frequencies (solid vertical lines) at a fixed time t, separated by the cavity- 
mode spacing A c . (b) The output frequencies (dashed vertical lines) after 
the field undergoes a delay in time, as would occur with a reflected field. 
Shown is one possible value of yT that would produce this sequence, along 
with a single unshifted component (solid line); other possibilities differ by 
n d c , with integer n. (c) Arrows mark some of the frequency differences that 
are beat-frequency signals, observable as resonances. Those pointing to 
the right contribute to the "direct" resonances, those to the left contribute 
to the "reverse" resonances. The central rectangle encloses a frequency 
interval of size A c used for measurement. 



of different teeth of the reference and reflected fields. 
Fig. 4c shows some of these frequencies. Two classes occur, 
distinguished in the figure by arrows pointing to the right 
and arrows pointing to the left. These frequency differences 
are observed as beat signals at the frequencies 

<{q) =q*c + yT (2) 

and 

a)l(p) - pA c - yT (3) 

with integer q y p = 0, ±1,±2, . . . 

At each frequency a>^(n) the beating signal is a sum of 
many contributions from paired teeth of the reference 
and reflected fields (each with the same frequency difference 
(o%(n)). A very important feature of a FSF laser is that in 
the absence of fluctuations, e.g. those caused by spontane- 
ous emission, the phases of all these contributions are bal- 
anced in such a way that their sum is equal to zero. 
However, unavoidable fluctuations break this balance, 
and one can therefore observe a resonant increase in the 
spectral density of the RF spectrum of the Michelson inter- 
ferometer output intensity fluctuations near these beat fre- 
quencies. Such signals have been used for distance 
measurements in [16,17]. Recently we have proposed [24] 
and experimentally demonstrated with a Yb 3+ -doped 
fiber-ring FSF laser [18] another approach to distance mea- 
surements, based on the use of a phase-modulated laser as 
an external seed. When the output of such a laser passes 
through a Michelson interferometer, the phase modulation 
induces amplitude modulation. The variation of the magni- 
tude of the amplitude modulation as a function of the mod- 
ulation frequency Q reveals very sharp resonances. This 
occurs whenever Q coincides with the beat frequencies 
given by Eq. (2) or (3). 

From a measurement of the frequency associated with 
the maxima of the amplitude modulation signal we can 
deduce the delay T and hence the distance L. Because the 
monochromatic amplitude modulation of the FSF laser 
output can be measured with a very narrow-bandwidth fil- 
ter, the technique allows a significant decrease in the noise 
and a corresponding improvement in the signal-to-noise 
ratio. 

Fig. 5 shows the various resonance frequencies co^(q) 
and co^(p) for given delay time T, and hence for distance 
L = cT jl. The locus of values appear as lines marked by 
integers p and q. For best results the modulation frequency 
should be chosen within a region convenient for measure- 
ments. The most natural choice is a modulation frequency 
which is smaller than the free spectral range of the cavity, 
A c (or the cavity mode spacing) 



0 < Q < A c 



(4) 



shown hatched in Fig. 5. In this region, where L ^ knc/A, 
(k = 0, 1,2, , . .), there exist, for a given y7\ two resonance 
frequencies, Q = Q D and Q = Q Ry corresponding to specific 
values of integers p and q. With increasing distance L the 
"direct" resonance frequency Q D increases and the 
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Fig. 5. Possible frequencies of amplitude modulation resonances, Qo or 
Q R> as functions of yT = (2A/cz T )L. Thick lines, labelled by q y mark values 
for fl D , Eq. (5). Thin lines, labelled by p, mark values for f2 R , Eq. (6). The 
shaded portion marks the region 0 < Q < A c used for frequency 
modulation. 



"reverse" frequency Q R decreases. The relevant connection 
with distance comes from Eq. (2) or (3), namely 

Qu = q*c + (2A/cx r )L (5) 

and 

Q R =pA c -(2A/cx T )L. (6) 

As with any interferometry, there is an ambiguity in the 
present technique. Distances which differ by an integer 
number of intervals L A = nc/A give the same resonance fre- 
quencies in the interval 0 < Q < A c . This ambiguity is 
resolvable by measuring changes 6Q D or &Q R produced 
by changing the AOM frequency from A to A + 8 A 
[16,17]. The required theoretical relationship is obtained 
from Eqs. (5) and (6) for constant A Ci and with variable 
A. As long as the change dA of the AOM frequency is small 
enough that integers q and p do not change we obtain the 
distance L from the formula 

_ cr r 6Qo _ £tr 
T"oT~ 2~bA' 



(7) 



4. Variation of amplitude modulation signal with distance 

For practical applications it is important to know the 
degree of amplitude modulation at resonance. To derive 
it consider, following [24], the FSF laser seeded by a CW 
laser whose phase &{t) = w s t -h cp s (t) is modulated 
sinusoidally, 

<Ps(0 = <p 0 + fism(Qt). (8) 

Here co s is the seed laser frequency, /? is the modulation in- 
dex and (p 0 is the initial phase. To obtain the amplitude 
modulation signal we model the FSF-laser output as a dis- 
crete spectrum. In [15] we have shown that this model is 
equivalent to the moving comb model. In this model the 
FSF-laser output consists of a set of discrete optical com- 
ponents with the frequencies co s 4- nA (n = 0, 1,2, . . .). The 
phases of the components have a time-independent contri- 



bution that is quadratic in n and a time-delayed contribu- 
tion from the phase-modulated seed laser phase 
q> s (t - nx r ). The amplitudes of the components are propor- 
tional to a Gaussian function of n y 



On = On 



,exp [-(/i -n^) 2 /^], 



(9) 



where n w is the effective number of optical components 
within the Gaussian profile of the laser spectrum, n max iden- 
tifies the component closest to the maximum of the spectral 
profile and a max is the amplitude of the central component. 
The parameters n w and n max are large, n w » 1 and 
"max 3> 1 1 and can be evaluated from measurable laser 
characteristics (for details see Eqs. ( 1 3)— < 17) in [15]). 

When this output enters a Michelson interferometer 
whose arms differ in length by L, the interference term in 
the interferometer output intensity is a sum of independent 
contributions 3P\„l from each nth frequency component. 
Every contribution is proportional to a cosines of the phase 
difference between reference and reflected components 

</>„ = (co s + nA)T + <p s (t + T- nx r ) - cp s (t - nx x ) (10) 

and to the square of the component amplitude: 

SPgocfl^cos^,. (11) 

Because the phase (p s (t) of the seed laser is modulated, the 
difference phase <f>„ will also be modulated, with a time- 
dependent part 

b(j> n (f ) = 2P sin [QT/2] cos [Qt + QT/2 - nQr T ] , (12) 

as follows from Eq. (8) and the use of a trigonometric iden- 
tity. Thus the contribution of each individual component 
has an amplitude modulation with the period 2n/Q. Mak- 
ing use of a Jacobi-Anger expansion [25] 



exp[iZcos(0)] = i 7 ^/(Z)exp[i/0] 



(13) 



one can consider this periodic modulation as a sum of the 
harmonic components with the frequencies IQ (/ = 1 , 2, . . .) 
and with the amplitudes proportional to the Bessel func- 
tion Jt(Z) of integer order / and of argument 
Z = 2j9sin(^r/2). The amplitude S n (Q, T) of the harmonic 
of present interest, the first harmonic, varies with time de- 
lay T = 2L/c and with frequency Q in accordance with the 
formula, 



S„(Q, T) = s 0 J l (2p sin [QT/2])a n 2, 



(14) 



where (Z) is the Bessel function of order 1 and so is a nor- 
malization factor that incorporates the laser output power. 
For particular combinations of the modulation index /?, the 
modulation frequency Q and the time delay T, the ampli- 
tude modulation with frequency Q may actually be sup- 
pressed for each elementary interference term. This will 
occur when 



2p sin [or/2]=y M> 



(15) 
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Fig. 6. Density plot of the calculated amplitude-modulation signal 5 A as a function of distance L and the modulation index 0. The left frame corresponds 
to the "direct" signal, the right frame corresponds to the "reverse" one. The basic length scale of the figure is U = ncf A = 1.875 m for our AOM 
frequency of A = 80 MHz. 



where y' M = 0,3.8317,7.0156,... are the kih zeros of the 
Bessel function of order 1 . 

If the condition (15) is not fulfilled it is easy to see from 
(10) and (12) that there are two sets of elementary contribu- 
tions to the signal of the amplitude modulation with 
depending on n phases: a "direct" signal with phase 
= n(Qx r + AT) and a "reverse" signal having phase 
= n{Qx T -AT). 

The sum of these contributions, weighted by the 
smoothly varying squared amplitude a 2 n will be zero for 
all Q except for very narrow intervals near the frequencies 
Q p providing that 

Q p x T ± AT = 2np (16) 

with p = 0, ±1, ±2 . . . Note that the frequencies Q p deter- 
mined by this condition coincide with the beating frequen- 
cies of Eqs. (2) and (3) obtained from the moving comb 
model. 

Regarded as a function of modulation frequency fi, the 
amplitude modulation signal has a Gaussian profile: it has 
the width AQ = 2y/2/(x r n w ) and reaches a maximum 



S A (Q, T) = SoJ\ {2fi sin [QT/2]) 



(17) 



with 5*0 — n w s 0 y/n/2 i when the modulation frequency Q 
has the value determined by Eq. (16). Thus by measuring 
the modulation frequency for which the amplitude modula- 
tion of the interferometer output maximizes, and knowing 
the chirp rate y, we obtain the interferometer-arm differ- 
ence L. 

Fig. 6 shows a density plot of the amplitude modulation 
signal S Ay calculated from Eq. (17), as a function of the dis- 
tance L and the modulation index p for the "direct" and 
"reverse" signal. The basic scale length of the interferome- 
ter-arm difference L 0 for this plot is determined by the con- 
dition yT = AT/x r = A c . Using the relationships T = 2L/c 
and A c = 2n/x r we obtain L = nc/A. Thus, the scale length 
in Fig. 6 is determined by the AOM frequency A. Our 
AOM frequency is 80 MHz, leading to the distance scale 
of L Q — nc/A = 1.875 m, apparent in the figure. 



The Bessel function in the expression (17) for the signal 
produces a strong variation of the signal 5a with distance L 
and modulation index p. Fig. 6 illustrates this rather com- 
plicated behavior. The regions of minimal signal (dark 
blue 1 regions in Fig. 6) pose challenges for detection. How- 
ever, these can be circumvented by proper selection of the 
modulation index P at each distance L. 

The preceding analysis assumed a monochromatic seed. 
In reality the seed laser has a finite spectral width, set by 
the finite coherence length of the seed laser. This will 
decrease the signal (17) in a manner that can be approxi- 
mately described by a damping factor exp(-L/£ coh ), where 
Lcoh is the coherence length of the seed laser: 



S A = SqJ} (20 sin [QL/c]) exp {-L/L^). 



(18) 



The inclusion of a finite bandwidth of the seed laser does 
not alter the Bessel-function variation of the signal ampli- 
tude with distance L. 

5. Experimental results and comparison with theory 

Fig. 7 shows the experimental setup we used for the dis- 
tance measurements. It consists of an all-fiber FSF laser 
seeded by a phase-modulated external monochromatic diode 
laser. The output power of the FSF laser passes through a 
Michelson interferometer with one fixed and one movable 
arm. The distance difference between them can be changed 
from 1 5 cm to 5 m. After emerging from the Michelson inter- 
ferometer the laser light impinges on a photodetector with 
bandwidth 50 MHz connected to a lock-in amplifier. The 
latter is driven by the same generator that operates the 
electro-optic modulator. Thus the lock-in amplifier detects 
resonances of the amplitude modulation Sam that corre- 
spond to various distances L. We observed the spectrum of 
the FSF laser using a Fabry-Perot interferometer (FPI) with 
the free spectral range of 650 GHz, see Fig. 7. 



1 For interpretation of color in Fig. 6, the reader is referred to the web 
version of this article. 
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Fig. 7. Experimental setup for distance measurements. The upper 
rectangle encloses the FSF laser of Fig. 1. The seed laser is modulated 
by an EOM. The output coupler feeds 10% to the Michelson interferom- 
eter (rectangular box at bottom of the figure), one arm of which includes 
the length L to be measured. The output from a beam splitter (BS) in the 
interferometer is detected as an amplitude modulation signal 5a with a 
lock-in amplifier matched to the frequency of the EOM. The seed input 
coupler also provides output analyzed by a Fabry-Perot interfermometer 
(FPI). 

Fig. 8 shows a typical example of a narrow "direct" res- 
onance in the magnitude S A as a function of the modula- 
tion frequency Q. Because in this case the distance L is 
smaller than nc/A the frequency of the resonance maxi- 
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Fig. 8. Amplitude modulation signal S A versus modulation frequency Q. 

mum is directly proportional to L. For our laser system 
the proportionality constant is 8.316 kHz/mm. The 
FWHM of the resonance is about 7 kHz. Thus, the resolu- 
tion available with our system is better than 1 mm and the 
accuracy is better than 0.1 mm (assuming that the fre- 
quency position of the resonance maximum can be deter- 
mined within 10% of its width). 

As was mentioned above, the amplitude modulation sig- 
nal S A has a complicated dependance upon the distance L 
and the index of modulation /?, described by formula (18). 
Fig. 9a and b compare the experimental results with 
the theoretical calculation made for our laser system. The 
experimental results are in excellent agreement with the 
theory. 




4.0 4.5 5.0 5.5 



Fig. 9. Dependance of the amplitude modulation signal size upon the distance for a modulation index 0 = 1.35 (frame a) and /? = 2.7 (frame b). Dots are 
experimental results and lines are theoretical simulation. Black squares and thick lines refer to the "direct" maxima (D, see Fig. 5), red circles and thin lines 
refer to "reverse" maxima (R). 
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Fig. 9 shows that there are two minima of the "direct" 
signal for ($ = 2.7, at the distances 3.6 and 4.8 m. These 
do not exist for ft = 1.35. This means that for any distance 
L there is an optimal modulation index 0 such that it is 
possible to avoid signal minima. Thus the useful range of 
our method is limited only by the coherence length of the 
seed laser (this was 3.6 m in our work), which is set by 
the bandwidth of the seed, and the collection efficiency of 
the reflected light. For our setup this latter restriction is 
several tens of meters. 

6. Summary and conclusions 

This paper describes an ail-fiber FSF laser with Er 3+ - 
doped fiber as the gain medium and an acousto-optic mod- 
ulator as the frequency shifter. The spectral and power 
characteristics are presented for cases with and without 
an external seed laser; these are typical of FSF lasers. 

We have used our system, a FSF laser with an external 
PM seed laser, to measure distances with a Michelson inter- 
ferometer scheme. We display examples of the amplitude- 
modulation signal for distances ranging from 15 cm to 
5 m. Our results are in very good quantitative agreement 
with numerical simulations which take the finite bandwidth 
of the seed laser into account. We thereby verify the theory 
first presented in [24]. Notably, the signal exhibits pro- 
nounced minima at particular distances, set by the zeros 
of a Bessel function. These "blind" zones can be avoided 
by changing the modulation index of the PM seed laser 
and/or the AOM frequency. 

The all-fiber FSF laser seeded by an external PM laser is 
a powerful noise-immune tool for distance measurements 
.over a range of distances up to tens of meters. The band- 
width limitation of our present apparatus allows an accu- 
racy better than 0.1 mm and resolution better than 1 mm. 

Dual wavelength interferometric measurements offer the 
opportunity to determine lengths to within a small fraction 
of a wavelength, as numerous publications have attested 
[26-29]. Such high accuracy is only possible under carefully 
controlled conditions. However, measurements of dis- 
tances to arbitrary objects require detailed information 
about the surface, because there occurs a material-depen- 
dent phase shift that affects any measurement based on 
phases, such as interferometry [30]. Techniques based on 
the use of a FSF laser, being essentially measurements of 
time rather than phase, do not suffer from this problem. 
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ABSTRACT 



We present a detailed theoretical analysis, using correlation functions, of the coherence properties of the 
output from a frequency shifted feedback (FSF) laser seeded simultaneously by an external seed laser and 
by spontaneous emission (SE). We show that the output of a FSF laser is a cyclostationary process, for 
which the second-order correlation function is not stationary, but periodic. However, a period-averaged 
correlation function can be used to analyze the optical spectrum. From the fourth-order correlation func- 
tion of the output of a Michelson interferometer we obtain the essential characteristics of the radio-fre- 
quency (RF) spectrum, needed for describing the use of the FSF laser for optical-ranging metrology. We 
show that, even for a FSF laser seeded by SE, the RF spectrum comprises a sequence of doublets, whose 
separation gives directly a measure of the length difference between the interferometer arms. This dou- 
blet structure is a result of the correlation of interference terms of individual components of the cyclo- 
stationary stochastic process. It is not seen in the optical spectrum of the FSF laser but is observable in 
the RF spectrum. We analyze the competition between SE and continuous wave (CW) seeding to obtain 
an analytical expression for the ratio of power in the discrete CW signal to the background continuum 
spectrum from SE. We show that, unlike mode competition in conventional lasers, where there occurs 
exponential selectivity, here there is a balance between the two fields; the power in the fields is related 
linearly, rather than exponentially, to the control parameters. 

© 2008 Elsevier B.V. All rights reserved. 



1. Introduction 

The remarkable properties of frequency shifted feedback (FSF) 
lasers have drawn attention as broad-band light sources [1-4] 
and, more recently, as tools, for metrology [5-9]. In essence, an 
FSF laser comprises a closed optical path, of total optical length I 
(a linear cavity or a ring), in which light undergoes not only the 
usual gain of energy from excited states but also, during each cav- 
ity round trip, a discrete-frequency shift A. This differs from the ax- 
ial mode spacing of the cavity, the free-spectral range (FSR), which 
is 2nc/L. The operation of FSF lasers has been considered by several 
research groups [10-20] and it has been used for a number of prac- 
tical applications [21-26J. 

Two complementary viewpoints offer means for understanding 
the output characteristics, and uses, for an FSF laser. In one, the 
output field is regarded as a comb of chirped frequencies, sepa- 
rated in frequency at any instant by the free-spectral range 
27tc/L (the moving-comb model) [5,6.14-16.27]. The other approach 
considers the output growing from a discrete-frequency seed as a 
stationary set of frequencies starting from the seed frequency 
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and separated from each other by A [19,20,23,28.29]. Spontaneous 
emission (SE) provides a continuous distribution of such seed 
frequencies. 

For applications to metrology, specifically the determination of 
distances by optical ranging, it is important to understand the 
coherence characteristics of the output from an FSF laser. In the ab- 
sence of any external seed, this output originates from a distribu- 
tion of spontaneous-emission (SE) events within the gain 
medium. These are uncorrelated and intrinsically incoherent. We 
here examine, by following a sequence of cyclic passages within 
the FSF cavity, how this radiation acquires coherence properties, 
as evidenced in correlation functions. 

We discuss here the operation of a FSF laser whose field origi- 
nates either entirely from spontaneous emission or else from the 
joint action of spontaneous emission and a continuous-wave 
(CW) monochromatic laser seed. We pay particular attention to 
the RF spectrum of the output of a Michelson interferometer driven 
by a FSF laser. As we will note, this signal underlies metrology 
applications. 

Earlier we have presented the basic theory of the operation of a 
FSF laser [19]. That work treated the SE as a seed field of very large 
bandwidth. It also described the effect of an external CW laser seed 
field, one whose phase could be modulated. However, it did not 
discuss the competition between these two sources. Here we rem- 
edy these shortcomings by treating the two sources as coexisting. 
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All applications of FSF lasers to distance measurements (optical 
ranging) rely on measurements of the RF spectrum of the output of 
a Michelson interferometer with FSF-laser input. Hitherto there 
has been no quantitative ab initio theory of the properties of this 
spectrum. Such a description must originate with intensity correla- 
tion functions, as presented here. We will show how the fre- 
quency-shifting element produces a field which, when passed 
through a Michelson interferometer, produces spectral signals 
from which the interferometer-arm distance can be deduced. Spe- 
cifically, the information is contained in sets of frequency doublets. 
We present here a first-principles derivation of the "beat frequen- 
cies" that have been discussed earlier based on the moving-comb 
model of FSF operation. We show that such signals occur even 
when the FSF laser is seeded entirely by spontaneous emission, 
and that there is no need to invoke a moving comb model to ex- 
plain this structure. The interferences from such SE sources are 
not seen in the optical spectrum of the FSF laser; they are observa- 
ble only in the RF signals from the interferometer. 

Our earlier work |19] considered the FSF output field from a 
specific seed field. We used these results to model the output from 
an FSF laser seeded either solely by SE or exclusively by an exter- 
nally controlled laser field. As shown in [7), the latter option offers 
significant improvement in signal-to-noise ratio over SE seeding. 
The present paper treats FSF-laser operation when an external seed 
field is accompanied by the inevitable spontaneous emission that 
provides an unavoidable stochastic background. In normal laser 
operation the competition between modes leads to exponential 
growth of favored modes at the expense of those with less gain. 
We shall show that in the FSF laser the competition between SE 
and external seed does not act in this way. Instead, there occurs 
a balance between the two fields that is related linearly, rather 
than exponentially, to the control parameters. Although it would 
be desirable to suppress the amplified spontaneous emission 
(ASE) component of the FSF-laser output, thereby improving the 
signal-to-noise ratio, we shall show that this is not possible. This 
behavior differs, perhaps unexpectedly, from that of a normal laser, 
where suitable choice of gain characteristics can place all of the 
output into a favored mode. 

2. Basic model 

This article draws upon two recent papers [7,19] that present, in 
some detail, an idealized model for an FSF laser. We here first re- 
view the basic principles of that work and then extend it by exam- 
ining correlation functions, thereby allowing analysis of the 
coherence and spectral properties of the FSF laser. 

2.1. The FSF-laser layout 

Fig. 1 presents a schematic drawing of the essentials of a FSF la- 
ser in a ring configuration. Light propagates along a closed optical 
path that passes through a gain medium (G). a frequency shifter 



seed ^ 




output ^ 




AOM 





(usually accomplished by an acousto-optic modulator. (AOM)] 
and a filter (F). In C a combination of stimulated emission and 
spontaneous emission add energy. The AOM adds a fixed incre- 
ment A to each frequency component. A number of elements with- 
in the overall optical path have frequency dependent losses that 
we here represent as a single spectral filter. F . The output emerges 
through the zero diffraction order of the AOM. The AOM also serves 
to couple any external seed laser beam into the cavity. 

2.2. The field description 

The filtering elements (F). which remove energy, and the gain 
medium (C), which replenishes this energy, together provide an 
effective gain (the growth minus the loss) distributed over a range 
of frequencies. We take the center of this range to be co ft regarding 
this as a carrier frequency. We consider the propagation of the 
electric field along a one-dimensional path, coordinate z, within 
the ring. The round-trip time in the cavity is z r = L/c. 

We denote by E(t) the complex-valued electric field at position 
z = 0. taken to be the entrance surface of the AOM. We express this 
field in terms of a complex-valued envelope <f(r) and a carrier: 

E(t)-*(t)exp(-fa*t). (1) 

The intensity of the FSF laser at z = 0. is 

/(0=^l*(0l 2 - (2) 

The AOM diffracts into first order a fraction of the beam inten- 
sity, denoted & 2 . Typically & « 0.9. The frequency of this fraction 
is shifted, in passing through the AOM, by A. This portion continues 
to circulate within the cavity. The unshifted fraction (1 - g 2 ), usu- 
ally small, remains in zeroth order and forms the output field. 

We consider an arbitrary but fixed reference time T and a time 
window that endures for one round-trip time prior to T. Within 
this interval, T - x r < t < T, we introduce a sliding windowed Fou- 
rier transform (SWFT) of the envelope, writing 



(3) 



(4) 



g(t) = [ dwS(w, T) exp(-irot), for T - x r < t < T, 
where 

& (ro, T) = — I dts(t) exp(+rmr). 

Here and in the following we denote by m = co - a>f a frequency 
offset from the carrier co f . The SWFT field 6 (07, 7) for a fixed refer- 
ence time T provides the essential tool with which we express the 
operation of the FSF and the output field [ 1 9]. From it we construct 
the field envelope 6(t) for all time, and deduce the output charac- 
teristics of the FSF for any given experimental conditions. 

2.3. The field amplitude equations 

As in our earlier work (19], we derive the basic equation for the 
laser field by considering its history during the round trip prior to 
time T. It began this circuit, starting at time t = 7*- T r . with fre- 
quency co = coj + w - A. It then underwent a frequency shift A as 
it passed through the AOM. It changed amplitude in passing 
through the filter and gain medium. It also acquired a phase incre- 
ment (ojf + ru)z r as it propagated around the loop. To describe 
these changes it proves useful to introduce a two-dimensional 
time-frequency space X = (n?, T) in which to present the equations 
[19]. Successive circulations around the ring involve a change of 
the two-dimensional vector X by the increment X 0 = (A,T r ). The 
field change after one round trip is expressible as 



Fig. 1. Symbolic diagram of a ring cavity showing gain C. spectral filter F and 
frequency-shifter AOM as well as seed and output. 



S(X) = 6{X - X 0 ) exp [c(X)] + £(X) + e(X). 



(5) 
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Here the effect of gain and loss appears through the propagator 
argument G(X), described in detail below, and the terms {(X) and 
£(X) represent, in Fourier space, additions to the field during the 
time interval T r . The first of these, expresses the field incre- 
ment supplied by spontaneous emission. It is the SWFT of ^(t), 
the amplitude of the stochastic field created by spontaneous-emis- 
sion events, 

c(X) = z(mj)=±j^ dtt sp {t)exp(+iwt). (6) 

The second contribution to the field, £(X). describes the external 
seed field, if present. We take this to have a carrier frequency co 5 
offset from the central filter frequency <o f by w St 

co s = (Of + w s . (7) 

We are interested in seeding by monochromatic light; we write its 
SWFT as 



dt exp [-im s t + imt], 



(8) 



where e, is the amplitude of the seeding laser field within the 
cavity. 

We write the propagator argument C(X) as the difference be- 
tween gain g(a>) and losses /(a>), modified by the phase 
(ajf + w)r r acquired during the round trip. We approximate the ef- 
fect of losses by means of a filter function, taken to be a quadratic 
centered at a> = co f with characteristic width //, 



f(co) =/ m + 



1 ( aj-w f \ 2 

2\ r, J ' 



(9) 



The frequency-independent term f m describes fixed losses, such as 
those from mirror reflectivities. 

We will regard the spectral filter as having narrower bandwidth 
than the broader bandwidth of the gain medium; it fixes the range 
of frequencies we need to consider. In this approximation we ne- 
glect the frequency dependence of the gain and approximate it as 

go 



gM - g sat = 



l+7// M 



(10) 



Here g sat is the saturated frequency-independent gain, g 0 is the 
unsaturated gain, 7 is the laser intensity averaged over one round 
trip and I sat is the saturation intensity. The averaged intensity, and 
hence the saturated gain, must be evaluated from appropriate equa- 
tions describing the FSF laser, see Section 5.1 below. Combining 
these three contributions, we write the propagator argument as 

C(X) = G(w,T) = i(o)f + w)x r +g sat -f m -~(TD/r f ) 2 . (11) 

This depends on frequency co through the detuning tu = a> - cof. 

The basic parameters under the control of the experimenter are 
those of the cavity r r . A, those of the filter ojj, F f , / m , those of the 
gain g 0 , /jar. and those of the seed, £,, co s . The mean intensity /, 
and the concomitant saturated gain g sat must be determined such 
that the equations for the field, with predicted output intensity, 
are self consistent. 

2.4. The output intensity 

An important characteristic of the FSF-laser output is the inten- 
sity averaged over a round trip. 



T r Jt-t, 



(12) 



We define the spectral intensity l{w, T) by writing this average as a 
frequency distribution. 



J(T)= rdwl(w,T). 

J-ao 

Using Eq. (3) we write the intensity as 

'(0 = 8^7 J dWldW2 exp[-i(07i - ufi)t)] 



x &(m u T)6{m 2 JY 



(13) 



(14) 



From this integral we extract an expression for the spectral inten- 
sity of Eq. (13) as 



/(ro, T)=S- J™ duf<f{m, T)S(w + ar\ T)' 
1 - expl-itp'Trl 



lTO'T r 



exp(ia7T]. 



(15) 



The fields described here, and the corresponding output inten- 
sity, involve uncontrollable events, the stochastic processes that 
describe spontaneous emission. Observations inevitably introduce 
averages over stochastic realizations of the average intensity dur- 
ing a round trip. Thus we are interested in stochastic averages 
(■ ■ ) that give the spectrum 

J(w) = (l(wJ)), (16) 
and the mean intensity 

J = <7(T)>- (17) 

2.5. The generol solution for the output field 

The FSF laser based on the growth of spontaneous emission acts 
as a regenerative amplifier of spontaneous emission, as described 
by basic Eq. (5). The solution to this equation, as discussed earlier 
[19], is expressible as the contribution of two parts, 



^(X) = ^ p (X) + ^ ed (X) 1 

a field that grows from the spontaneous emission, 

= £(*) + £ £(* - nXo) exp[A„(X)], 
and a field that grows from the external seed. 
<WX) = e(X) + £ e(X - nX 0 ) exp[A„(X)]. 



Here 

A n (X) = £G(X-iXo). 

1=0 



(18) 



(19) 



(20) 



(21) 



These formulas express the field as the most recent added incre- 
ments from seed and spontaneous emission, together with the suc- 
cession of fields from earlier passes, each altered by the 
appropriate propagator exp[>L n (X)]. Our earlier work considered 
each of these separately. Here we consider the operation of the la- 
ser when both are present. 

3. Seeding by spontaneous emission. The optical spectrum 

Earlier work on the theory of the FSF laser emphasized the 
structure of the output fields. To gain a full understanding of the 
FSF laser it is necessary to evaluate various spectral properties. 
For the field emerging from the FSF cavity the significant observa- 
ble is the optical spectrum. When the FSF laser is used for 
metrology, the significant observable is the RF spectrum from a 
Michelson interferometer. To quantify either of these it is neces- 
sary to have a theoretical description of various correlation 
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functions. Specifically, we require the second-order correlation 
functions for the optical spectrum, considered in this section. We 
require fourth-order functions to describe the RF spectrum, dis- 
cussed in Section 4. 

We begin by considering the FSF laser seeded only by SE. This 
provides a simple foundation for more elaborate models of FSF-la- 
ser operation. 

3. J. Passive frequency-independent cavity 

To simplify the problem as much as possible, we first consider a 
passive-cavity model, one in which there is no gain and in which 
losses are independent of frequency. We assume that the cavity 
field originates with the intracavity spontaneous emission. The 
lack of frequency dependence allows usto construct the solution 
using the fact that after one round trip the field E(t) reproduces, 
with no distortion, the field at the earlier time t - T r . The results 
of this simplified approach clarify some of the apparent peculiari- 
ties of the FSF-laser field. 

For our idealized cavity the main losses occur from the small 
nonideality of the AOM, whose efficiency & — exp(-2/ m ) < 1 is 
taken here to be independent of frequency. We assume that the 
spontaneous emission source in the cavity creates the stochastic 
electric field with the amplitude f, p (t) and the carrier optical fre- 
quency co S p. We take f sp (r) to have mean value zero and the corre- 
lation behavior 



<6p(t)€i(r + T)) = ^F(T). 



(22) 



To avoid the unphysical occurrence of an infinite energy and 
infinitely broad spectrum we take the function F(t) to be 



F(t) = ?f expt-r^Tl + i{w 5p - OJ/lT), 



(23) 



rather than the usual delta function F(t) = <5(t). Here r sp and co^ 
are, respectively, the spectral width and the central frequency of 
the spontaneous-emission line, much larger than all others frequen- 
cies in the system. (For the moment we disregard gain; the assump- 
tion of very large r, p is equivalent to our earlier assumption of 
constant gain within the spectral range of interest.) 

Because we assume here that the cavity losses are independent 
of frequency we can write the amplitude S{t) of the electric field in 
the cavity as the sum of the immediate spontaneous-emission con- 
tribution £ jp (f) together with emission that was present at the ear- 
lier time t - r r along with similar emissions at a succession of 
earlier times, each diminished by a diffraction SI but with un- 
changed time dependance, 

S (r) = &p(t) + #&p(r - T r )e~ i&t + ^^(t - 2T r )e- fA e- T '»- fAf , 

+ «»{ v (t - 3T r )e- iA( '- 2T '>'-* (t - T '>'-'" + - * • (24) 
Thus the cavity field is an infinite sum 



*(t) = 6p(0 + W - nr r )e- 

n=1 



(25) 



where the phase of component n is <P„ = XX~o' At " 

We infer from Eq. (25) that the spectrum of the field is the sum 
of spontaneous-emission spectra weighted with and shifted by 
nA. We verify this by examining the second-order correlation func- 
tion for the electric field amplitude. 



G< 2) (r, r + t) = (<f (r)<r (r + t)) exp(fo>/T). 



(26) 



The term exp(ico/T) provides a shift to the optical frequency. 

It is customary to deal with stationary processes for which cor- 
relation functions depend on the time interval t but not on the 
time r. For such processes one can evaluate the power spectrum 



as the Fourier transform of the correlation function; that is the con- 
tent of the Wiener-Khinchine (WK) theorem. For the system con- 
sidered here this theorem does not apply. Instead, we use Eq. 
(25) and the correlation behavior of the spontaneous-emission 
source, Eq. (22), to obtain the result 

C< 2) (t,f + t) = tl ^ k F(x - [n - k)x r ) x exp[iS(r,n,fc) + ic^t] 

n.fc-0 

(27) 

where 

S(r, n,k) = -(n - k)At + kAx + <*>„ - <J> k . (28) 

This correlation function G <2, (r, f + t) depends explicitly on time 
t. This means that the electric field in the FSF laser is not a station- 
ary process. However, this time dependence is periodic, with peri- 
od 2 n/ A. Thus we have here an example of a stochastic process 
which is periodically correlated (or cyclostationary, periodically sta- 
tionary (30]). For such a process the spectrum can be defined as the 
output of a narrow-bandwidth optical filter. The output of such a 
filter, for such a process, is obtained from the Fourier transform 
of the period-averaged correlation function [30], When we average 
over time the terms of Eq. (27) with k^n cancel. leading to the 
result 



g< 2 >(T) = C (2) (t,t-i-T) = 51) ^ 2k F(T)e lkAT+, ^ T . (29) 

fc=0 

The laser spectrum is the Fourier transform of this time-aver- 
aged correlation function, 

/ 2) (o>) = ^^^~dT^ 2 )(T)exp(-^T). (30) 

This can be considered as an extension of the WK theorem to cyclo- 
stationary processes. 

Because the correlation function g {2) (x) is the sum of indepen- 
dent terms, so too is the spectrum the sum of separate spectra. 

#s P or>t(CO) - ^ jf~ dzFlT)***-*", (31) 

each offset by an increment A and weighted by a power of # 2 . 



f\oj) = ^^^spondCO - co 5p - kA). 

This result confirms our expectation based upon Eq. (25). 



(32) 



3.2. Frequency-dependent effective gain 

We next consider a more realistic model of the FSF laser by 
including gain and frequency-dependent losses. The field is given 
by Eq. (19). We assume that the width r sp of the spontaneous- 
emission spectrum is much larger than the width r f of the intra- 
cavity filter. Thus, the spontaneous-emission process f sp (r) can be 
considered here as a delta-correlated (Wiener-Levy) stochastic 
process: 

(33) 



When we use this property of £ jp (f) to calculate the field, using 
Eq. (19) we obtain the correlation function 

'n-1 n+M-1 

xexp 5Zfcar -Urn - \A)] + £ l&« ~f( w - 

. 1=0 /=0 



x exp 



-/A)]] 

+l(G7 + oj f )Sz 4- iMAf + ^ AM 2 T r - ~ AMr r j . 



(34) 
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We have here defined M and <5t by expressing r as r = Mx r + St 
with -T r /2 < <5t < T r /2. 

As in a passive FSF cavity, the correlation function G ,2, (f,f + r) 
given by Eq. (34) depends explicitly on time f. As discussed above, 
to obtain the spectrum we need the period-averaged correlation 
function. When we average over time the terms of Eq. (34) with 
M ^ 0 cancel. The resulting average vanishes unless -T r /2 < t < 
T r /2, when it gives the result 

s (2) w = i&i r dw s ex pi 5 " + * w + (35) 

The amplitude of component n here depends on the exponent 

Because the spectral width of J(vu) is much larger than 2n/x r 
this averaged correlation function is just the Fourier transform of 
the spectral density. 



^g (2) (T) = jT dmf\m) exp[+i(n7 + o*)t|, 



with 



(37) 



(38) 



This result was obtained earlier [19] using the definition (15) of the 
optical spectral density. 

4. Seeding by spontaneous emission. The radio-frequency 
spectrum 

As we have shown above, the FSF laser output is an example of a 
cyclostationary stochastic process |30]. This property has no sub- 
stantial consequences for the second-order correlation functions, 
such as those associated with the optical spectrum, other than 
the need to average over a period. However, the cyclostationarity 
has significant consequences in the fourth-order correlation func- 
tions. These determine, for example, the RF spectrum of the FSF 
output intensity. In this section we analyze such signals. 

4.1. Passive, frequency-independent cavity 

The output intensity of the idealized FSF cavity, seeded only by 
spontaneous emission, will obviously be noisy. The RF spectrum of 
the fluctuations of the output intensity is of special interest for 
applications to metrology, where it is observed the output of a 
Michelson interferometer. Observations of the RF spectrum were 
important in the use by Nakamura et al. of a FSF laser, not seeded 
from an external source, for ranging measurements [5,6). RF spec- 
tral measurements continue to be central to the use of externally 
seeded FSF lasers for ranging [9]. 

The RF spectrum derives from the second-order correlation 
function of the intensity, i.e. the fourth-order correlation function 
of the electric field of Eq. (25).The intensity is 

l(t) = ±4(t)r(t) 

= i^F £ a"** W - n ^W 5p (t - kz r ) x e i *- k >"- i *°+ i *> . (39) 

n,fc=0 

We define the correlation function for this intensity, based on 
deviation from the mean value J, as 

C (4) (t, T) = {[l(t) - HV(t + T) - J]) = </(t)/(t + T)) - P- (40) 

We express the averages as the multiple sums 



and 



P = (C/871) 2 



(/(r)/(t + T)) = (c/8ti) 2 

x (Zs P (t - n?rK p (t - kx r )Z sp {t + x- mJrHlp 

x (t + t - st r )^ x exp [i(n - k)At + i(m - s) A(r + i) 

- i<P n + i<P k - i<P m + 1*0,] (41 ) 

£ ^ k U 5P (t - nx r )^(t - fc Tr ))e l * (n - k)AC - | * B ^l - 

n,k=0 J 

(42) 

To average the product of four random functions we need to in- 
voke additional assumptions about the statistical behavior of the 
function £ sp (t). We will assume that this is a Gaussian random pro- 
cess [31]. This assumption allows the factorization 

<^p(ti)^ p (t 2 )^(t 3 )<j; p (t4)) = (^p(ti)^(t 2 ))<C sp (t3)^ p (t4)) 

+ <^(t,K; p (t4))<^(t 3 K; p (t2))^ 

(43) 

The width of the spontaneous-emission spectral distribution is 
much larger than all the relevant frequency parameters. Therefore 
we assume that the x dependence of the function F(t) in Eq. (23) is 
much more localized than are the variations with x in Eq. (41 ). We 
carry out the ensemble average of the particular products shown 
above and find that the fourth-order correlation function is inde- 
pendent of time: 



tf 4) (t f T) = G (4) (T) = Co ]T \F{X - MT r )| 2 ^ 2!M| 
with 



(44) 



(45) 



This correlation function does not depend on the frequency shift 
A. Therefore particular characteristics of the FSF behavior is not 
important; we obtain the same result by considering a cavity that 
has no frequency shifter. The periodic structure of the correlation 
functions occur because the field repeats itself after one round trip. 
The field at any time results builds on superpositions of many con- 
tributions, each shifted in time by x r and weighted by the reflection 
coefficient The interference terms involving different compo- 
nents do not contribute to the mean intensity but they do give 
important contributions to the correlation function. This correla- 
tion function is independent of time and so. by the WK theorem, 
the RF spectrum J l4) (Q) is its Fourier transform. 



J (4) (0)=^ |^dTC (4 >(T)exp(-iflT). 



(46) 



Thus the spectrum of output intensity fluctuations is the infinite 
sum 



J 



dTexp(-ifiT) |F(T-MT r )| 2 # 2|M| . 



(47) 



We write this as the product of a scaling factor G 0 and two fre- 
quency-dependent functions, 

p A) (Q) = G Q U p (Q)<?{Qy (48) 
Here the smooth envelope 

wot = & £ <*\m? ex P( -,«t) = ± 1+( y; rtp)2 m 



♦ 
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incorporates the broad but finite width of the spontaneous-emis- 
sion spectrum. For the RF spectrum it suffices to use the approxima- 
tion Lsp(fl) * Ljp(O) = r sp /Bn. The sum over M in the spectrum (47) 
produces the factor 



1 -i* 4 



(50) 



This is a periodic function of *2 with a period 27t/T r . For a high 
quality cavity, one with 1 - & «: 1 , this function describes a set 
of narrow lines, each centered near a value Q q = 2nq/x r for some 
positive integer q. The line profile ¥{Q - Q q ) of each individual 
component is a Lorentzian 

^> = (i-V+(«/rc)Y (51) 

where r c = 2(1 - #)/r r . 

Thus the RF spectrum of the passive cavity output comprises a 
sequence of equidistant spectral features, each with the same nar- 
row profile i?(x), and each centered at an integer multiple of the 
cavity axial-mode frequency-difference 2n/x r . The width of each 
peak is determined by the lifetime of a photon in the cavity. The 
amplitudes of the peaks differ very little. The origin of this spec- 
trum regularity is a periodic repetition of any initial fluctuation 
of the spontaneous emission source. Note, as remarked above, that 
the FSF behavior of the cavity has no effect on these properties of 
the RF spectrum. 

Frames (a) and (d) of Fig. 2 illustrate the properties of the cor- 
relation function C } (x) of Eq. (44) and the spectrum S(Q) of Eq. 
(47). For this figure the various times have been chosen to illus- 
trate the essential physics rather than to describe any actual real- 
ization. The individual peaks of frame (a) have widths set by the 
atomic spontaneous emission lifetime; in practice those widths 
would be much less than the spacing, and the figure would show 
only vertical lines. The dashed-line envelope has a decay time set 
by the lifetime of a photon in the cavity; in practice this decay 
would be much more gradual than is shown. The spectral proper- 
ties of frame (d) are similarly exaggerated. The widths of the indi- 



vidual peaks are set by the cavity decay time, while the dashed-line 
envelope expresses the spontaneous-emission lifetime. 

42. Michelson interferometer with input from a passive cavity 

Here we analyze the RF spectrum of the output of Michelson 
interferometer. For simplicity we consider the interferometer 
formed by an ideal beam splitter and two totally reflecting mirrors. 
We take the difference of interferometer arms to be the distance 
L - cT/2, where T is a total time delay. 

The electric field in the output arm is 



*MW.(r) = j[*(t)+*(t-T)]. 

The interferometer output intensity, 

can be written as a sum 

WO = i|/(t) + /(t - T) + / jm (t, t - T)] t 



(52) 



(53) 



(54) 



where /(t) is the output intensity given by Eq. (39). The interference 
term 



/»{f,t-T) = 2»[*(ry'(t-T)] 



(55) 



makes no contribution to the mean intensity but, as we will show, it 
has important effects on the RF spectrum of the Michelson-interfer- 
ometer output The correlation function for the intensity of Eq. (54), 



C£L(t,T) = (Wt)Wt + T))-]J 1 



(56) 



J MWi^» W - VMich\t>)'MichK*- T *,)/ - 'Mich 

has two distinct contributions, which we denote as G l *l h (t, t) = 
C3 4) (t, x) + G\%(t y x). The first of these is the sum of three individual 
correlation functions 



C<">(r, T) = CSV) = ^ [2C (4, (T) + C«>(T + T) + C'V - T)] , 



(57) 





1 2 3 4 S 



3 



e 




-6-5-4-3 -2 



Fig. 2. Left-hand frames (a), (b). (c) show the correlation function, right-hand frames (d), (e). (0 show the resulting spectrum. Top row, frames (a) and (b) show the second- 
order correlation function of Eq. (44) and the corresponding optical spectrum of Eq. (47). Middle row. frames (b) and (e) show the correlation function of Eq. (57) and the 
correspondingRF spectrum ./J 4 ' of Eq. (62). Bottom row. frames (c) and (d)show the correlation function GjjJ/ ofEq. (58) and the corresponding RF spectrum./^ of Eq. (64). The 
observable spectrum is the sum of contributions (e) and (0- 
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and is independent of time. The second part originates with the 
interference intensities. 



1 



(58) 



This interference term can be averaged in the same way as was 
done when obtaining C* 4) (t) in the previous subsection. The result is 

Cffitt. t) = Cffl(T) = Co £ F (* - Mr r ) 2 ^' M ' cos[MArj. (59) 

M=-oo 

The additional phase shift AT between consecutive components in 
the cavity leads to some modulation of the correlation function with 
a period 2nx r /6J. 

The RF spectrum of the output from the interferometer receives 
contributions from two parts, corresponding to the two parts of the 
correlation function, 

jSL(0)=^ ) (fi)+vS5(fi)- (so) 

The contribution } { ^(Q) comprises a sequence of equidistant spec- 
tral profiles, each with the same profile i?(ft) but with different 
amplitude: 

» = 3^ 

x f~ dz[2G iA) (x) + C (4) (t + T) + G (4) (t - T)j exp(-iflr) 

= £ ® 2m rdxexp(-iQx) > 

x [2F(T - Mr r ) 2 + F(T + T - Mx r ) 2 + F(x-T- Mx r ) 2 ] . 

(61) 

This function is expressible as the product of a scaling factor and 
three frequency-dependent functions (in contrast with the two 
functions of Eq. (48)], 



(62) 



The first factor, U P {Q). can be considered constant, in keeping 
with our assumption of large r sp . The second factor. &{Q). is that 
of Eq. (48): it describes a set of narrow peaks, centered at the dis- 
crete frequencies Q Q = 2nq/x r . In contrast to the previous section, 
here there occurs an additional frequency dependence: amplitudes 
of the narrow discrete components are weighted by the factor 
cos 2 (OT/2). evaluated at Q q . This revision occurs because the cor- 
relation function C 3 originates from an incoherent sum of three 
components delayed by time T. Fig. 2b illustrates the properties 
of the correlation function G^x) of Eq. (57) and of the spectrum 
y^ 4) (0) of Eq. (62). The latter exhibits an overall damped cosine 
modulation, as expected. 

The contribution f*](Q) to the RF spectrum from C{ 4) . discussed 
next, is qualitatively different. From the definition we obtain the 
expression 

-C(G) = oT- r dxG ° £ F(T-MT r ) 2 ^ |M| cos|MA7]exp(-i^T) 
= ^L $p (fl) £ ^ |M, cos[MAnexp(iMfiT r ). (63) 
This spectral function therefore comprises a set of doublets, 
jffiOT = £ f§^(0)JSf(O - 2nq/x T - yj) 

+ ^(Q)<e(Q - 2nq/x r + y e T)] . (64) 



Fig. 2f illustrates this spectral structure. The quantity y e = A/T r 
appearing here is recognizable as the chirp rate used for optical 
ranging. 

To summarize: the total RF spectrum of the Michelson-interfer- 
ometer output consists of narrow discrete components centered at 
the frequencies Q q - 2nq/x r . each weighted by a factor 
cos 2 (fl g T/2). Each of these components is accompanied by a dou- 
blet of components shifted from Q q by ±y c T and whose amplitudes 
are almost independent of q. The frequency separation of these 
doublets incorporates information about the difference in length 
cT/2 of the interferometer arms that are used in ranging applica- 
tions of FSF lasers [5-9|. Fig. 2 illustrates, on the third line of 
frames, these features. 

It is worth noting that the existence in the RF spectrum of com- 
ponents shifted by ±y c T is readily understood from a model of a 
FSF laser seeded by a phase fluctuating CW laser [7]. Here we see 
that these components appear in a model of the passive FSF cavity 
seeded by uncorrelated spontaneous emission as a result of inter- 
ference of cyclostationary stochastic processes. 

4.3. Interferometer input from a frequency-dependent cavity 

Here we consider a more realistic model of the FSF laser by 
including both gain and frequency-dependent losses. The field is 
given by Eq. (19). We will continue with assumptions made in Sec- 
tion 3.2: the width r sp of the spontaneous-emission spectrum is 
much larger than the width r f of the intracavity filter, so we ide- 
alize the spontaneous-emission process £ 5p (t) as a delta-correlated 
stochastic process; its correlation function is given by Eq. (33). In 
addition, we assume that the width r f is much larger than both 
the frequency shift A and the cavity axial-mode frequency-differ- 
ence 2n/x T . 

Using the solution (19) for the FSF laser field and the condition 
ry » A,27i/T r we obtain, after some complicated but straightfor- 
ward algebra, the following expression for the intensity correlation 
function C (4) (t,T) defined by Eq. (40): 

C (4) (t,T) = tf 4> (T) 

= Uj<5 J_l g exp(im5T) g exp[S(n7,n,M)| 



(65) 



Here the argument of the exponent is zero unless 0 < t - Mx r < x r 
when it is 



S(uj, n, M) = |> fl , -f(w - /A)] + "£'"(g $ct -f{w - /A)]. (66) 



n-1+M 

£ 

1=0 



To simplify this expression, we use the condition r f » A,27c/t r 
and replace the summation over n in (65) by integration. The expo- 
nent S(fi7,n,M) reaches a maximum at 



n = rw t + M/2, 



with td m = yjxol - M 2 A 2 /4. In the Gaussian approximation justified 

in our earlier paper [19] the parameter m Q = y j2{g sat - f m )rj is the 
frequency shift between the maximum of the Gaussian fit to the 
optical spectrum J {2) {co) = J 0 exp [-(to - w 0 - eo/) 2 /f 2 ] and the fil- 
ter frequency a) f . As shown in [19], and justified experimentally in 
[20], the shift w 0 and the width T of the FSF-laser spectrum are re- 
lated by the expression 

f 2 t37o = y 3 = (Ar 2 ). (67) 
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Using the Gaussian approximation we then write 
jr exp[S(tn t n t M)] ~Aexp\S(w,n max ,M)]. (68) 



The phase and amplitude occurring on the right hand side are 
-07 3 + 3tcrgtP + 2wl M 2 



3y> 



(69) 




(70) 



(71) 

is the effective number of round trips of a SE photon in the FSF-laser 
cavity. Integration over xo gives finally 



. 4M oJ 



(72) 

(73) 
(74) 



C< 4 »(t) = C 1 J F,(T-MT r ) 2 exp 
where 

Q -"5b U^J exp [ V 3 . 
f((T ) = _exp^j. 

The correlation function (72) is very similar to the correlation 
function of the output intensity of the passive FSF cavity, Eq. 
(44), but there are two differences: 

1. When gains and losses are present the width of the narrow 
function F,(t) is determined by the reciprocal of f , the width 
of the optical spectrum. In the passive cavity this width is fixed 
by the spontaneous-emission correlation time. 

2. The successive maxima of the correlation-function decrease in 
keeping with a Gaussian law exp(-M 2 /Mj) with M 0 T r equal 
to the effective lifetime of a photon in the active FSF cavity. 

It can be shown that the same is true for the correlation func- 
tion of the Michelson-interferometer output excited by the active 
FSF laser: the contribution G^ 4) (t) is given by Eq. (57) with C 4) (t) 
defined by Eq. (59) and the interference term 



G! 4 r(*) = C, £ F,(T-Mt r ) 2 exp 



cos [MAT]. 



(75) 



Thus we conclude that the RF spectrum of the FSF laser output 
intensity and the Michelson-interferometer output have the same 
behavior as in the passive cavity. The only difference is that the 
amplitudes of narrow discrete components centered at the fre- 
quencies Q q = 2nq/Xr and fl^= 2nq/x r ± yj decrease according 
to the Gaussian law exp(-^/r 2 ) and the shape se{Q - Q q ) of each 
narrow component is not Lorentzian but a Gaussian: 

(76) 



jfi?(x) oc exp [-x 2 /(M 0 T r ) 2 ] . 



5. FSF laser with seeding by both spontaneous emission and 
continuous-wave radiation 

5.1. The intensity equation 

We consider here an FSF laser seeded simultaneously by spon- 
taneous emission and a monochromatic seed laser. Because there 
is no nonlinearity in our model (only saturation of gain by the total 
intensity) the laser field is the sum of two fields described by Eqs. 



(19) and (20), each created by one of the sources. These sources are 
statistically independent and so the averaged intensity is also a 
sum of spontaneous and discrete spectrum intensities, 

7 = 7sp+7„ed. (77) 

To evaluate the contribution of the external seed, / Sfrd . we con- 
sider a monochromatic seed, meaning that its bandwidth is much 
smaller than either the axial mode spacing 27t/T r or the frequency 
shift A. For such a seed the field ^^(t) is (19) 



n=0 



(78) 



where 

* n = -T r n[oj s + (n + 1)A/2]. (79) 

We write the constant real-valued amplitudes as 
a n = exp[S w «, (g sa „ n)], where 

Sseed(gsat, n ) 

6w 2 + 6m s nA + 6ro s A + 2n 2 A 2 + 3A 2 n + A 2 



gsat —fm 



(80) 

The contribution l iee(t of the discrete spectrum to the mean 
intensity 7 reads 



- oo 

heed = g^s«d(t)| 2 = IseedYlM 2 = 1 seed F seed (gsat ), 



(81) 



where the intensity of the seed laser radiation injected into the FSF 
laser cavity is 



i--£i*r" 

and 

oo 

Fseedigsat) = £ exp[2S seed (g Mt , n)]. 

n=0 



(82) 



(83) 



The mean intensity of the spontaneous-emission spectrum. 

2 



has been evaluated earlier [19]. It can be written as 

fjp = UpF sp(g; 0 t)' 

Here 



(84) 
(85) 

*-tsf (86) 

is the intensity of spontaneous emission within the spectral interval 



r s and 

Fsp(g 5Q t)=^-f2 ex PP*(&*.n)]. 



with 

Ssp(gsat>n) = n 



2(g«-/m)- 



(n 2 -l) 



- ln(n)/2. 



(87) 



(88) 



Because the stochastic averaging removes any interference con- 
tribution the sum of the two terms (81) and (85) gives the total 
intensity 7. It is this intensity that saturates the gain 
gsat = go/0 +7//wf). The determination of 7 = hat{f£ - 1) requires 
the solution of the following transcendental equation for g sat : 



'sot hat \Ssat J 



(89) 
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Our interest is with FSF lasers in which the output spectrum is 
both smooth and broad. Such conditions are possible if J c . 
'o lsot(*1 - !) and if tne spectral filter is sufficiently broad. 
r f » A. For these conditions the general operation of the FSF laser 
is readily understood from Fig. 3. The figure shows the cavity gain 
as a constant, and the loss as varying around the central frequency. 
The frequency shifter displaces any seed frequency toward the 
right. By contrast to the operation of a conventional laser, there ex- 
ists a finite band of frequencies a> f - ro 0 < co < co f + w 0t symmet- 
rically centered about the central frequency of the filter co/, 
where the frequency-dependent effective gain is positive. 
Ssat -/m - Jr* > °- Anv fr et l uenc y component within this band will 
undergo groWth as successive round trips increase the frequency 
by J. However, once the frequency exceeds co f + w 0 the losses 
dominate, and the frequency component will diminish as further 
frequency shifts occur. Thus we expect the saturated gain g sat to 
be slightly greater than the minimum loss,/ m . 

To proceed we make two approximations. First, because we ex- 
pect the saturated gain g sat to be slightly greater than f m we 
approximate it as g sat = f m on the right hand side of Eq. (89). Sec- 
ond, as in our earlier work [19], we assume that the seed frequency 
is far from the central frequency of the laser spectrum, so that the 
laser output spectrum is very close to a Gaussian form. 

We express Sspig^^n) and S s «d(g sar ,n) as a Taylor series, cen- 
tered around the value n = n„ where dSj{n)/dt = 0 (i ~ sp or 
i = seed), and we retain only terms through second order in 
(n - tii). Within this Gaussian approximation the number of dis- 
crete components in the interval ru 0 is large, m 0 /A :» 1 . and so 
we replace the summation over n in Eqs. (83) and (87) by an inte- 
gration. We thereby obtain the results 



F S e e d y/nMo exp 



'2ml + U}3 s ~ 3rc s m 0 



3y3 



and 



(90) 



(91) 



Eq. (89) now reads 



Here the small parameter 0. defined as 

-ft /q/> 
p "2/ Sflf tt-l)y' 



(93) 
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Fig. 3. Qualitative portrait of the magnitude of the FSF laser output-field spectrat- 
density as a function of frequency, together with frequency dependence of the gain 
(constant) and loss (quadratic). Two choices of seed frequency are marked. 



describes the spontaneous-emission source. The parameter 

(94 > 

expressing the relative importance of the seed laser, is the ratio of 
the intensity of the seed laser to the spontaneous-emission inten- 
sity in the spectral interval y/ny/2. Using experimental data |20) 
we estimate that for an Yb 3+ -fiber laser the parameter e is equal 
to 1 when the power of the seed laser injected into the FSF cavity 
is about 50 u,W. 

5.2. Competition between spontaneous emission and continuous-wave 
seeding 

In this subsection we discuss solutions to Eq. (92). We consider 
two specific choices for the seed frequency. In the first, marked a in 
Fig. 3, the seed frequency is at the minimum of the loss curve, 
ct s = 0. In the second, marked b in the figure, the seed is near the 
frequency where growth first becomes possible. 

5.2. J. Resonant seed, co s = a) f 

When the seed frequency coincides with the minimum of 
losses, tu, = 0, Eq. (92) takes the form of a quadratic equation 



Y 2 + eY = 



ft 



involving the variable 
From the solution, 

y 2/fi 

e+s/e 2 +4/p } 
we obtain the intensity of the continuous spectrum l sp as 
4/g 



r7/«, t (>?-l) 



and the intensity of the discrete spectrum as 
2e 



heed — 



Wf-l). 



(95) 



(96) 



(97) 



(98) 



(99) 



We see that I sraf > I sp for e > 1 / y/2fi » 1 . This is to be expected: 
when the seed frequency is resonant with the loss minimum. 
cd s = then the seed laser must compete with spontaneous 
emission that has been amplified throughout the frequency inter- 
val from (Of - rz7 0 to co/. 

5.2.2. Seed at gain edge. cv 5 = a>f - ta 0 

A better choice for the seed laser is (o t - co f - w 0 . Eq. (92) then 
reads 



(100) 



The solution to this equation for w 0 with 0 «: 1 gives the following 
expressions for the intensity l sp of the continuous spectrum 

1 



/sat(?7-l) 

and the intensity 7 5eed of the discrete spectrum 

Mi? - 1) 



'seed = 



1 + tsfo 



(101) 



(102) 
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where the parameter a is defined as 

-(3 ta i) • (103) 

This parameter is almost independent of laser parameters. In 
previously reported experimental work |20] the measured value 
of a was (7^2. The discrete spectrum dominates in the case of a 
detuned seed laser starting with much smaller intensities: 
> hp for e > l/v/ff- 1. 

We see that even for optimal detuning of the seed laser there 
exists a large background of continuous radiation originating from 
spontaneous emission. For a strong seed laser, e » 1. the ratio of 
spontaneous to seed intensity is 

U- (104) 
Led erf 

This means that, roughly speaking, the ratio of continuous to dis- 
crete spectra is equal to ratio of seed laser intensity to the intensity 
of spontaneous emission in the spectral interval xn 0 . 

6. Summary and conclusions 

We provide here the first complete theory of the spontaneous- 
emission seeded FSF laser based on correlation functions. We pres- 
ent the second-order correlation function, from which one obtains 
the optical spectrum of the radiation emerging from the FSF cavity. 
Although the correlation function is not stationary, it is periodic, 
and so a period-averaged correlation function leads to a simple 
expression for the spectrum. 

We also present the fourth-order correlation function analysis 
with which one can evaluate the RF spectrum of the laser output 
intensity and of the output of a Michelson interferometer, as is 
used in optical-ranging metrology. We show that even in a FSF la- 
ser seeded solely by spontaneous emission this spectrum com- 
prises a set of doublets, whose spacing gives directly a measure 
of the length difference between the interferometer arms. The exis- 
tence of this structure in the RF spectrum of a the Michelson-inter- 
ferometer output results from correlation of interference terms of 
individual components of a cyclostationary stochastic processes, 
and is not seen in the optical spectrum of the FSF laser. 

We discuss the operation of a FSF laser whose field originates 
from the joint action of spontaneous emission and a continuous- 
wave monochromatic laser seed. We show that in the FSF laser 
the competition between SE and external seed does not act in 
the usual way, of exponential growth of modes with greater gain. 
Instead, there occurs a balance between the two fields that is re- 
lated linearly, rather than exponentially, to the control parameters. 
Although it would be desirable to suppress the spontaneous-emis- 
sion component of the FSF laser output, thereby improving the sig- 
nal-to-noise ratio for optical-ranging technique based on phase 



modulation of the CW seed laser, we show that this is not possible. 
This behavior differs from that of a normal laser. 

We have not here discussed the temporal properties of the laser 
output. The statistical approach we have used in this paper, based 
on an analysis of the correlation functions and the spectral charac- 
teristics, is necessary because we consider a FSF laser seeded by 
spontaneous emission - a sequence of events that is completely 
chaotic. A time series of the laser output, just like the time series 
of the seeding spontaneous emission, is chaotic and conveys no 
information. 
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1. Introduction 

Many optical ranging technologies make use of interference be- 
tween two ray paths, one of known distance and the other to a tar- 
get object, to measure distances with subwavelength accuracy. 
Lasers have served well as the radiation sources for optical ranging. 
In particular, frequency shifted feedback (FSF) lasers have been 
proposed and demonstrated 11-7]; their properties offer potential 
advantages over other laser-based techniques for ranging [8-13]. 

In recent papers |1 0,1 1,13] we have proposed and demon- 
strated the use of a FSF laser whose seed frequency is phase 
modulated. This easily implemented modification improves the 
signal-to-noise ratio by many orders of magnitude and thus over- 
comes some recognized problems. Although the mathematics has 
been presented in detail [10], and experimental demonstrations 
have been published [1 1,13], the mathematical formalism has pro- 
vided no intuitive picture of the role of phase modulation. 

We here offer a very simple presentation of the principles in- 
volved in the phase-modulated FSF laser and its use for optical rang- 
ing. To explain the technology we first review the basic operation of 
a Michelson interferometer, because this forms the heart of the 
ranging measurement. We start by pointing out. as is well known, 
how classical interferometry with monochromatic light has an 
intrinsic ambiguity in measuring lengths: these are determined 
only within a multiple of the optical wavelength. One technique 
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for overcoming this limitation is the use of phase-modulated light 
as interferometer input, but this technique has limitations that we 
mention. 

Having presented this background, we discuss how a phase- 
modulated FSF laser, being the superposition of a very large num- 
ber (typically more than 10 4 ) of mutually coherent components of 
a radiation field, overcomes the shortcomings of a conventional la- 
ser source when applied to ranging. It relies on a measurement of 
frequency, which is inherently more accurate than measurements 
of response amplitudes. As will be explained, the FSF laser with 
phase modulation greatly improves the signal-to-noise ratio. 

7.7. Basic Michelson interferometer 

The schematic diagram of Fig. 1 shows the essentials of a 
Michelson interferometer: from a source of collimated radiation 
two paths, whose lengths differ by L, recombine at a detector that 
produces a signal proportional to the power, i.e. proportional to the 
square of the local electric field. The two contributions to the de- 
tected field are derived from a common source field but differ by 
the delay time T = 2L/c Thus the detector signal has the form 

S(tJ) = Jf{\E(t)\ 2 + \E(t - T)\ 2 + 2Re[£(t)£(t - 7*)*]}. (1) 

Here, the factor Jf expresses the efficiency of the detectionprocess 
as related to the square of the electric field. Note that when the 
material of the target differs from that of the reference there occurs 
upon reflection a difference in phase that acts as an increment of 
length [14]; we here ignore this phase difference. 
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Reference 



Source 



path 2 i 



Unknown 



Detector 

Fig. 1. Basic schematic layout of a Michelson interferometer, showing light source, 
power detector, and two arms. For clarity the overlapping beams are shown here as 
slightly offset spatially. 



The simplest implementation of such a Michelson interferome- 
ter occurs when the source is monochromatic, of wavelength 
X = 2nc/Q). We write this field, at the position of the detector, in 
the complex form 



£(t) = £exp[-ia>t). 



(2) 



Here, <f is the amplitude of the source field and co its frequency. It is 
this field that, after passing through the two paths of the interfer- 
ometer, appears at the detector to produce a signal. 

With such a field the detector signal is independent of time t 
and is expressible as 



S(t,r) = #(l +cos(coT)], 



(3) 



where coT = 4nL/X. The constant = 2\6\ 2 ^ determines the size of 
the signal for a given electric field. Measurements of this signal, for 
varying optical frequency a>, provide values of cos(o/T) and with 
that a measure of the ratio L/X. Thus the technique is capable of pro- 
viding values of L with an accuracy better than the optical wave- 
length A, but with an ambiguity (of A) arising from the periodicity 
of the cosine. When used for measurements of distance a classical 
interferometers cannot give values of absolute distance, only dis- 
tances relative to some integer multiple of the wavelength. There 
are many methods capable of resolving this ambiguity [15-17]. 
One of them uses a phase-modulated light source. The following 
paragraphs review the operation of that technique, as a preliminary 
to the discussion of phase modulation in FSF lasers. 

1.2. Absolute distance measurements using a phase-modulated light 
source 

One of the suggestions for making absolute measurements of 
distance using an interferometer came from Webb et al. [18], 
who suggested the use of a phase-modulated laser as interferome- 
ter input. We here follow their approach and consider a light 
source upon which we impose a specified phase variation (p{t), 

E(t) = <?exp[-icor - i<p(t)]. (4) 

The detector signal is then expressible as 

S{t, T) = ^{1 + cos[a>7 + (p(t) - (p(t - 7")]}. 

A particularly useful phase variation is a sinusoidal modulation, 
achievable using an electro-optical modulator driven at a radio fre- 
quency G, 



<p(t) = (p 0 + ps\n(Qt+d). (6) 

Here, p is the modulation index, <p 0 is the initial optical phase and d 
is the initial phase of the modulation signal. 

For such harmonic modulation of the phase <p(r) the phase dif- 
ference <p(t) - (p{t - T) is also sinusoidal 

(p{t) - <p(t - T) = 0{sin(flt + 1?) - sin(«(t - T) + tf]} 

= 2^sin(OT/2)cos(fl(t - 7/2) + 0} 

= zsin{Qt + <t>), 



(7) 
(8) 

(9) 

(10) 
(11) 



with the amplitude and phase 
z = 20sin(GT/2), <f> = - (QT + 7t)/2. 
The resulting interferometer signal is 
S(tJ) = <<?{! + cos[cuT + zsin(f?t + <f>)]} 
= <€{\ + cos{coT)cos[zsin(Qt + 4>)] 
- s\n{coT) sin[zsin(Gt + </>)]}• 
By using the Jacobi-Anger identities 

sin(zsinx) = 2 £j 2k _ 1 (z)sin[(2fe - 

oo 

cos(zsinx) =; 0 (z) +2]T; 2k (z)sin(2fc;r;)> 
we rewrite the signal as [18] 

S(r, T) = <e\ 1 + cos(ojT)J 0 (z) + 2 cos(coT) £j 2k (z) s'm[2kQt + 2k<f>\] 

-2 sin(coT) £j 2k _ t (z) sin[(2k - 1 ) Qt + (2k - 1 )<j>\ \ . (12) 

The signal of Eq. (12) is a superposition of the harmonics of the 
modulation frequency Q, each weighted by an appropriate Bessel 
function. Viewed in the frequency domain these appear as discrete 
components. We can use a narrow-band filter to select any individ- 
ual harmonic, using phase-locked detection or other means. In this 
way we can restrict consideration to a single value of k. 

For example, consider a filter that selects just the first-har- 
monic, k =.1. The signal responsible for this harmonic is 

S } (t, T) = -2<$J X (z) sin(o)T) s'm(Qt + <f>) 

= <ej } (z)[cos(«t + (f> + a>T) - cos(Qt + 0 - a>T)]. (13) 

From the measured amplitude of this signal we obtain the value of a 
Bessel function (z). This gives information about z, from which we 
obtain information about the absolute value of the delay T. Because 
the modulation frequency Q is radiofrequency (RF), it has a wave- 
length that is orders of magnitude longer than the optical wave- 
length of the light source. It is this long wavelength, rather than 
the optical wavelength, that determines the ambiguity with which 
we measure the distance L 

This method of ranging has two drawbacks. First, we have to ad- 
just the interferometer to maximize sin(a>T) in Eq. (13). Second, 
the method relies on a quantitative measurement of light intensity. 
From the magnitude of the signal we obtain the argument of the 
Bessel function z and, from Eq. (8), the time delay T. It is difficult 
to calibrate measurements of intensity, and hence the method 
has intrinsically low accuracy. These drawbacks have prevented 
practical application of the method for ranging using Michelson 
interferometer and phase-modulated light. 

( 5 ) 2. Optica] ranging using a FSF laser 

The development of frequency shifted feedback (FSF) lasers 
[1-13] offers a means of overcoming the previous limitations of 
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RF-modulated interferometry and of obtaining very high accuracy. 
The basic principles underlying a FSF laser are readily understood: 
it comprises an optical cavity (or closed loop) in which there occurs 
a frequency-shifting element, typically an acousto-optic modulator 
(AOM). and material that provides gain to a range of frequencies. A 
small portion of the cavity radiation emerges as output; see Fig. 2. 

From a single seed field of frequency co, and phase <p s (t) a suc- 
cession of round trips within the bounded laser cavity (or fiber 
loop) produce a set of frequencies, spaced equidistantly by the fre- 
quency shift A induced by the AOM. Typically the gain bandwidth 
is sufficient to amplify a large number (>10 4 ) of discrete frequency 
components, each with definite frequency and phase. These form a 
frequency comb spaced by the AOM frequency A. 

As we explain below, instead of one laser field an FSF laser pro- 
vides a large number of field components, each of which gives a 
similar contribution to the first-harmonic signal, but with different 
phases. One might expect that the averaging of so many signals 
would eliminate any interference effects. On the contrary, we will 
show that, when the modulation frequency matches a delay- 
dependent value, the signal is greatly enhanced; from a measure- 
ment of this frequency one obtains a very accurate value of the de- 
lay T and hence of the distance. 

2.1. The FSF laser 




Fig. 2, Symbolic diagram of a ring cavity showing gain C. spectra! filter F and 
frequency-shifter AOM. as well as seed and output (from [71). 




Gain 



Fig. 3. Frequency components of FSF laser, spaced by A. as a function of frequency 
u>. Gain and filter loss are shown, as is the seed frequency co,. 



The output field from a FSF laser is the sum of many 
components 

E(0 = £V n exp[-i(fl>s + nA)t - i<P n - i(p s (t - nt)), (14) 

n 

Component n has a fixed frequency that has been shifted from w s by 
nA. It has the phase of the seed laser taken at an early time 
(p s (t - nx) and it has undergone an additional phase shift, 

■4>„ = nx[G>, + (n + 1W2). (15) 

where t is the round-trip time within the FSF cavity. Gain from the 
amplifying medium and loss from the filter alter the original seed 
amplitude S 0 to the value <f„. As earlier work has shown [5], the dis- 
tribution of amplitudes is well described by a Gaussian centered at 
a peak value rw and of width n w . In this approximation the ampli- 
tude S n is 



& n - #o exp 



(n - rw) 



n 2 w 



(16) 



Fig. 3 illustrates the distribution of field components, showing the 
varying amplitudes, spaced by the AOM frequency A. The figure 
shows a representative distribution of gain and of loss from the fil- 
ter. The peak of the Gaussian, parametrized by n max . occurs close to 
the frequency at which gain equals loss. 

2.2. Coherence properties of the FSF laser output 

It is well known that the output of a multimode laser can. when 
the modes are phase-locked, appear as a periodic succession of 
very short pulses J 19]. The frequencies contributing to this output 
are separated by the free spectral range (FSR) of the laser cavity. 
27T/T for a round-trip time of t. and their phases vary linearly with 
frequency. A FSF laser can also operate in a regime in which the 
output is a train of short pulses [20]. The field is then a set of 
mode-locked frequencies spaced by the FSR, 2n/x. 

However, for use in optical ranging, the operating regime of the 
FSF laser is such that no mode-locking occurs. The frequencies of 
the output field, are spaced by the AOM frequency A, and their sta- 
tic phases d> n vary quadratically with component number n. The 
temporal output is periodic with the period Inj A and appears as 
a train of frequency-chirped pulses [5]. 



2.3. Optical ranging using a FSF seed having sinusoidally modulated 
phase 

The extension of the discussion of a single phase-modulated la- 
ser to a FSF laser is straightforward. Let us assume that the phase of 
the FSF seed laser is modulated as 

<p = 0sin(Gr + tf). (17) 

When such phase modulation of the seed is present the phase of a 
specific FSF component n reads 

cp n = 0sin(i?(r - m) + tf). (18) 
Therefore, we have 

cp n = psm(Qt + !?„), tf n =tf-nflT, (19) 

. meaning that each FSF component has. in addition to the static 
phase 0 n . a modulated phase <p n that originates with the seed-laser 
phase. The initial phase of that modulation varies linearly with 
component number n. 

Eq. (12) has given the expression for the signal for a single mod- 
ulated laser. With the FSF laser we have many such modulated 
components. The detector signal sums the results from each of 
them, appropriately phased. 

For example, consider filtering that selects just the modulation 
frequency Q (the first of the harmonics). Then the signal, general- 
izing Eq. (13), is the sum of all the individual components of the 
FSF field. 

S, (t, T) = 2JTU (z) ]T \S n \ 2 {zos(Qt + <f> n + [o) s + nA]T) 

n 

- cos(flr + 4> n - [o>5 + nd]T)}, (20) 

where 

0 n = i9 n -(OT + 7c)/2, 4 = tJ-nflr. (21) 

We rewrite this as two separate sums 

St (t, T) = 7.JTU W H l^nl 2 cos l^ " n ( QT ~ T ^ + *J 

n 

- 2jfU (2) Y, ^«i 2 cos i fit - n ( Qx + TA) + (22) 
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Fig. 4. Schematic distribution of vector phases, (a) a uniform distribution of vectors 
resulting from small constant phase increments between components of the 
radiation field, (b) A narrow distribution of phases, as occurs when the resonance 
condition is nearly but not exactly fulfilled, (c) All vectors have phases that obey the 
resonance condition. 



where 

tfr ± = ±co s T - (QT + n)/2 + tf. (23) 

Each of these sums is responsible for a distinct resonance (an 
enhancement of the signal within a narrow frequency range) as 
noted in the following paragraphs. 

2.4. Vector picture of interferometer output signal 

Because the terms of the summation have amplitude and phase, 
they can be regarded as vectors. When the modulation frequency Q 
obeys the resonance condition 



Qx = ±TA + m2n, m = 0,1,2 



(24) 



then all vectors in one of the sets are aligned in the same direction. 
The resulting vector has. as magnitude, the arithmetic sum of the 
amplitude of each of the many components. The resulting first-har- 
monic signal is 

S,(t,T) = ±2^(2) £ |^„| 2 . (25) 

n 

Fig. 4c illustrates this situation. By contrast with the results of a sin- 
gle laser, here we have the sum of many components (typically 
>10 4 ), and so the signal is greatly enhanced by comparison with a 
single laser. 

However, when the modulation frequency is nonresonant. so 
that the argument {Ox ± TA) is nonzero, the vectors point, with 
angular increment n(Qx±TA), into different directions. Their 
summed amplitude will therefore typically be small, or zero. 
Fig. 4a illustrates this situation. 

To evaluate the interferometer-arm delay T, from which we de- 
rive the ranging information, we must evaluate the modulation fre- 
quency Q for which the resonance condition of Eq. (24) holds. This 
is recognizable as the modulation frequency which maximizes the 
interferometer output signal at a frequency 

(26) 



, . , A „ 271 

The components will be in phase, and the signal will be large, only 
within a small interval 

c ^ 27t 7" 

6Q A. 

n w x 



(27) 



of frequencies near such a resonance frequency, set by the width n w 
of the component distribution, see Eq. (16). Fig. 4b illustrates this 
situation. 

Our optical ranging procedure, using a phase-modulated seed 
for a FSF laser, is based upon identifying the frequency where the 
resonance occurs. We do not need to measure relative intensity 
and interpret this through a Bessel function. It is this change, from 
measuring relative intensity to measuring a frequency (facilitated 



by strong signal enhancement at resonance), that gives our method 
such an advantage. 

Our procedure also dramatically enhances the signal observed 
with an unseeded FSF laser [8,9]. The phases of a FSF laser seeded 
by a phase-modulated input laser have well defined phases. By 
contrast, when the FSF grows from spontaneous emission the FSF 
field comprises a continuous distribution of frequency combs. 
Within any comb the phases are fixed by geometry, and are not 
adjustable, but the phases between different combs (started by 
spontaneous emission) are random. Our method uses the coher- 
ence of the various components to enhance the signal compared 
with growth from noise. 

2.5. Observing higher harmonics 

We remark that the resonances exist in the amplitudes of all 
harmonics of the output of Michelson interferometer excited by 
FSF laser with phase-modulated seed laser. For example, the signal 
of the second harmonic is 

S 2 (tJ) = 2jV] 2 (z) l^nl 2 sin[2f2t - n(2Qx - TA) + ^] 

n 

- 2JT} 2 {z) £ |^n| 2 sin[2Gt - n(2Qx + TA) + (28) 



where 

4, ± = ±co s T - (2QT + 7i)/2 + 2tf . 
The resonance condition is 
1 , 



(29) 



(30) 



Qx =-[TA + 27im\. 

In a similar way one can show that for the kth harmonic the reso- 
nance condition is 



Qx=l{TA + 2nm]. 



(31) 



Each harmonic signal varies in proportion to a Bessel function 
whose argument z, defined in Eq. (8), varies with distance L. When 
a distance is such that the Bessel function J, (z) for order 1 is small, 
thereby limiting the signal strength, the Bessel function J k (z) of a 
higher order will then give a larger signal. 



3. Summary and outlook 

Phase modulation interferometry has been earlier suggested 
and is well understood [18]. The physics of FSF laser operation is 
equally well understood [5]. The concept of ranging with a free- 
running FSF laser has been described and demonstrated some 
years ago by Japanese researchers [8.9]. However, the free-running 
FSF laser is not a practical device for industrial applications. 

The use of a phase-modulated FSF seed overcomes the draw- 
backs noted for a free-running FSF laser. As has been demon- 
strated, this technique provides dramatic enhancement of the 
interferometer output signal by orders of magnitude if. and only 
if, the modulation frequency Q obeys the resonance condition of 
Eq. (24). The accuracy of the ranging technique stems from the 
accuracy with which frequency measurements can be made, in 
contrast to measurements of relative intensity. 

Although the literature contains rigorous mathematical treat- 
ments of FSF lasers [10] the present discussion ofTers a simple 
straightforward picture of the physics behind the observation of 
the dramatic signal enhancement resulting from the phase modu- 
lation of the FSF seed. 
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